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PROGRAMA DE POSGRADO EN CIENCIAS DE LA SALUD 
ESPECIALIZACIÓN EN INMUNOHEMATOLOGÍA DIAGNÓSTICA 



Objetivo General 

Formar recursos humanos con capacidad para prestar servicios 
de diagnóstico hematológico e inmunológico, capaces de proponer 
alternativas para la investigación de problemas de diagnóstico 
así como para transmitir sus conocimientos mediante actividades 
docentes y de difusión de información científica. 



MAESTRIA EN CIENCIAS DE LA SALUD 



Objetivo General 

Formar recursos humanos con capacidad para realizar 
investigaciones de alto nivel en el campo de la salud y para transmitir 
sus conocimientos mediante actividades docentes y de difusión de 
resultados, con el interés de que su formación tenga un impacto 
decisivo en el perfil epidemiológico poblacional. 

DOCTORADO EN CIENCIAS 
(QUÍMICO BIOLÓGICAS Y DE LA SALUD) 

Objetivo General 

Formar personal altamente calificado en las áreas de química, 
biología y ciencias de la salud, con una visión interdisciplinaria, que 
contribuya a solventar la problemática de salud. 



PROGRAMA DE MAESTRIA EN CIENCIAS Y TECNOLOGIA 
DE LOS ALIMENTOS (CIENCIA Y TECNOLOGÍA DE GRANOS; 
CIENCIA Y TECNOLOGÍA DE PRODUCTOS MARINOS) 



Objetivo General 

Formación de personal calificado y competente de alto nivel 
que sea capaz de resolver problemas alimentarios en los sectores de 
investigación, educación, de los sectores público y privado, vinculados 
con el manejo, almacenamiento, conservación y procesamiento de 
granos o productos marinos. 



PROGRAMA DE DOCTORADO EN CIENCIAS DE LOS 
ALIMENTOS (CIENCIA Y TECNOLOGÍA DE GRANOS; 
CIENCIA Y TECNOLOGÍA DE PRODUCTOS MARINOS) 



Objetivo General 

Formación de personal de alto nivel académico que sea capaz de 
generar conocimiento, a través del planteamiento, diseño y desarrollo 
de proyectos de investigación relevantes al desarrollo de la ciencia de 
los alimentos y a la resolución de la problemática de la producción y 
procesamiento de alimentos 



POSGRADO EN BIOCIENCIAS (MAESTRIA) 



Objetivo general 

Formar recursos humanos con una visión amplia de su área de 
especialidad y con capacidad de identificar problemas y procesos 
básicos y aplicados en las diferentes disciplinas de las ciencias 
biológicas. Tendrá la capacidad de desarrollarse en áreas emergentes 
de las biociencias, combinando diferentes aproximaciones, 
conceptualizaciones y niveles jerárquicos de la biología y áreas de 
especialización del posgrado 



POSGRADO EN BIOCIENCIAS (DOCTORADO) 



Objetivo General 

Formar investigadores de alto nivel con una amplia capacidad 
de realizar investigación independiente, la cual contribuya 
significativamente a la generación, innovación y adaptación 
de conocimientos científicos y/o tecnológicos en su área de 
especialización. 

Posgrados no incluidos en el Padrón Nacional de Posgrados de 
Calidad de CONACYT 



MAESTRIA EN CIENCIAS AGROPECUARIAS 



Objetivos Generales 

1. Formar profesores investigadores con una alta preparación 
académica y capacidad innovadora, capaces de generar y transmitir 
el conocimiento científico; aptos para crear y desarrollar proyectos 
y programas de investigación científica y tecnológica, tendientes a 
mejorar la productividad agropecuaria nacional. 

2. - Preparar personal altamente capacitado para ejercer labores 
docentes de licenciatura y posgrado en las Instituciones de Educación 
Superior dedicadas a la generación de conocimientos y tecnología en 
las áreas de fisiología vegetal, genética y mejoramiento de plantas, 
manejo de agostaderos, praderas irrigadas, producción de carne y 
leche, nutrición y reproducción animal. 




http://www.investigacion.uson.mx/programas-posgrado.htm 
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En este primer semestre de 201 3, la H. Junta Universitaria 
convocó y desarrolló el proceso de designación de 
Rector para el periodo 2013-2017, el resultado fue 
designar al Dr. Heriberto Grijalva Monteverde para su segundo 
periodo como Rector de nuestra máxima casa de estudios, 
quien tomó posesión y rindió protesta el 1 6 de junio. 

Luego de protestar y recibir el nombramiento de parte del 
presidente en turno de la Junta Universitaria, Alfonso Rangel 
Guerra, el rector aseguró que encabezará una administración 
que busque los consensos para garantizar un clima de 
estabilidad laboral y su trabajo impulse la construcción de 
una sociedad con mejor calidad de vida e inclusión social. 

Los responsables de edición de la revista Epistemus 
hacemos llegar nuestra felicitación al Dr. Heriberto Grijalva 
Monteverde por lo que esperamos será un fructífero segundo 
periodo como Rector de la Universidad de Sonora. 

Por otra parte, también en este semestre se celebraron dos 
eventos de especial relevancia académica, el Departamento 
de Física de la División de Ciencias Exactas y Naturales convocó 
y realizó a nombre de la Universidad de Sonora y en conjunto 
con los Subsistemas de Educación Media Superior del Estado 
de Sonora, dependientes de la Secretaría de Educación y 
Cultura del Gobierno de Sonora, el Primer Congreso Estatal 
de Ciencias Exactas y Naturales y Encuentro de Enseñanza de 
las Ciencias, que se realizó del 30 de mayo al 1 de junio de 
2013, con el objetivo de compartir experiencias y estrategias 
sobre la enseñanza y el aprendizaje de las Ciencias Exactas y 
Naturales en los niveles básico, medio superior y superior, así 
como plantear, analizar y generar acciones de colaboración 
con los Sistemas de Educación Básica y Media Superior en 
este rubro. Además, promover la divulgación de trabajos 
de investigación en ciencias desarrollados en el estado de 
Sonora en las áreas de Física, Matemáticas y Ciencias Químico 
Biológicas. 

También en el presente semestre el Departamento 
de Investigación y Posgrado en Alimentos realizó el VIII 
Congreso del Noroeste y IV Nacional de Ciencias Alimentarias 
y Biotecnología, en San Carlos Nuevo Guaymas, Sonora, 
del 27de febrero al 1 de marzo de 2013, que consistió en 
varios cursos de actualización y presentación de trabajos de 
investigación científica en las áreasde las Ciencias Alimentarias 
y la Biotecnología, contando con la participación de cientos 
de asistentes de todo el país. 

Asimismo, es necesario informar que el Ing. Rafael 



Pacheco Rodríguez ha decidido concluir con su valiosa 
colaboración que abarca desde la creación y mantenimiento 
de la difusión y divulgación científica que han significado 
los 13 números en los que participó, que aunado al presente 
consigue Epistemus, un caluroso agradecimiento y el deseo 
de que continúe cosechando más logros. 

Contenido 

En este número 14 de Epirtemus se abordan diversos 
temas a través de sus secciones de Investigación, Desde 
la Academia, Políticas de Ciencia y Tecnología, y Ciencia, 
Tecnología y Sociedad. 

En la sección de Investigación se incluye Modelo 
no lineal con potencial seno-gordon para un cristal 
fotónico unidimensional. También Búsqueda de discos 
protoplanetarios en cúmulos abiertos jóvenes; Materiales 
ultrablandos: diagrama de fase de una suspensión coloidal 
de polímeros estrella; Modelo predictivo de riesgo de 
morosidad para créditos bancarios usando datos simulados; 
y Una aproximación al comportamiento de los valores de la 
vivienda: caso Hermosillo. 

En la sección Desde la Academia se presenta El modelo 
atómico de Bohr: una aplicación; Armonización de modelos 
de calidad del software; Diagnóstico de los instrumentos 
utilizados para la medición de la presión arterial en los 
pacientes atendidos en un centro médico; Diseño de 
alternativas de conservación de carreteras en el sur del 
estado de Sonora; Método alterno para determinar el costo 
de producción en líneas de ensamble; y Elementos para la 
presentación de las formas diferenciales. 

En la sección Políticas de Ciencia y Tecnología se 
incluye el trabajo intitulado La geometalurgia: una 
excelente herramienta multidisciplinaria para la vinculación 
universidad-industria. Este tema toma auge ya que la 
geometalurgia involucra a especialistas en el campo de la 
geología, la microscopía, la minería y la metalurgia extractiva 
para optimizar los procesos de explotación y beneficio de los 
minerales. 

En la sección Ciencia, Tecnología y Sociedad están 
publicándose los artículos: De meteoros y meteoritos; Usos 
y ocurrencia de los principales metales que se producen 
en Sonora y una reseña periodística sobre los Aniversarios 
astronómicos en Cananea, Sonora. 
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BASES 

^ PARA PARTICIPAR CON ARTÍCULOS O COLABORACIONES EN EPISTEMUS 



La educación, la ciencia y tecnología son consideradas como 
pilares fundamentales sobre los que se sustenta el desarrollo de 
un país; por eso es importante fortalecer el enlace entre los que 
generan el conocimiento y los beneficiados de ello: la sociedad. 
Uno de los aspectos que distingue a la Universidad de Sonora 
es la generación de conocimiento y su impacto en la sociedad. 
Es así como las Divisiones de Ingeniería, Ciencias Exactas y 
Naturales y Ciencias Biológicas y de la Salud crearon el proyecto 
editorial Epistemus el cual, constituye un medio de información y 
comunicación para dar a conocer las investigaciones realizadas o 
en proceso, las actividades académicas, las reflexiones en torno a 
la ciencia y tecnología, la cultura científica y la sociedad, proyectos 
de vinculación y extensión de la cultura y los servicios. 



OBJETIVO 



Divulgar el conocimiento que se crea, genera y enseña en la 
Universidad de Sonora con la finalidad de fortalecer la educación, 
la cultura científica, tecnológica y de la salud de los diversos 
sectores de la sociedad. 



ÁREAS GENERALES DE CONOCIMIENTO 



• Ingenierías: materiales, metalurgia, civil, minas, industrial, 
ambiental, hidráulica, sistemas de información, mecatrónica, 
alimentos, energía, agua, entre otras. 

• Ciencias exactas y naturales: geología, física, matemáticas, 
electrónica y ciencias de la computación. 

• Ciencias biológicas y de la salud: investigación en alimentos, 
desarrollo regional, acuacultura, salud, biología, agricultura, entre 
otras. 



DIRIGIDA A 



Los Sectores relacionados con la educación y la investigación, 
empresarios, dependencias gubernamentales, estudiantes de 
nivel medio superior, superior y posgrado, y sociedad en general. 



CONTENIDO DE LA REVISTA 



Artículos originales e inéditos de proyectos de investigación, 
reseñas, ensayos, información sobre ciencia y tecnología, eventos 
relevantes, convocatorias, noticias sobre educación y cultura. 



ENFOQUE DE LOS ARTÍCULOS 



• Los artículos reflejarán lo más trascendente de la producción 
académica, los proyectos de ciencia, tecnología y sociedad que 
se realizan en la universidad y en particular en las tres divisiones 
académicas. 

• Los artículos deberán de hacer mención de la trascendencia de lo 
expuesto, su impacto en la solución de problemáticas específicas 
de la sociedad, del sector industrial, de la educación, de la cultura, 
entre otras. 

• Se incluirán artículos que integren y reflexionen en torno a la 
ciencia, la tecnología y la sociedad, que aporten elementos precisos 
que permitan profundizar en el análisis y proponer esquemas 
de colaboración entre los que producen el conocimiento y los 
beneficiarios o usuarios potenciales. 

• El lenguaje escrito debe ser de buen nivel, con el rigor 
científico pero de divulgación, comprensible para un público no 
especializado de nivel bachillerato, empresarios y profesionistas 
de otras especialidades. 
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ARBITRAJE 

La revista es arbitrada en varios niveles. Comprende dos aspectos: 

• Arbitraje académico: En todos los casos, los artículos serán 
arbitrados por pares académicos internos y/o externos a la 
Universidad de Sonora . 

•Arbitraje de divulgación: Se evaluará el lenguaje y el significado 
de lo expuesto para que no pierda la idea original al tratarlos como 
artículos de divulgación. Participan en el arbitraje de académicos, 
expertos en divulgación científica. 

• Se les sugerirán adecuaciones y modificaciones para que sean 
tomadas en cuenta por los autores. En todos los casos se definirán 
los tiempos para las correcciones. 

CARACTERÍSTICAS DE LOS ARTÍCULOS 

Recepción 

El artículo deberá ser original e inédito y se hará llegar de manera 
electrónica al correo: revista.epistemus@correom.uson.mx 
El comité editorial hará una primera revisión del artículo, si 
cumple con el formato de la revista, se enviará al menos a 2 
árbitros académicos. Los arbitrajes se enviarán a los autores, 
quienes reenviarán los artículos corregidos. El comité editorial 
decidirá qué artículos se publicarán en cada número. Los autores 
de los artículos aceptados, deberán firmar una carta de cesión de 
derechos de autor y enviarla por e-mail al comité editorial. 
Título 

Utilizar letra tipo Arial de 12 puntos, negrita espacio sencillo, 
centrada. 

Autores, grado académico, afiliación, correo electrónico 

Incluir los datos curriculares de los autores, especificando su grado 
académico, institución para la cual labora y correo electrónico. 
Respecto a los autores, el primero será el autor y el resto serán 
coautores, siendo un máximo de cinco personas. 
Resumen 

Los artículos deberán de contener un resumen con letra cursiva, 
con un máximo de 200 palabras del total establecido. 
Contenido 

Utilizar letra tipo Arial de 12 puntos, normal, a doble espacio, 
justificada. La extensión de los artículos será de un máximo de 
3,300 palabras incluyendo el texto, figuras y citas bibliográficas, 
en Word con márgenes de 2.5 cm por lado. Los ensayos, reseñas 
e información de eventos académicos serán de un máximo 
de 1,500 palabras. Las ecuaciones deberán estar numeradas 
consecutivamente entre paréntesis. No incluir notas en pie de 
página. 
Secciones 

El autor señalará en qué sección de la revista desea se publique su 
artículo. Cuando la colaboración sea en la modalidad de ensayo, 
hacer uso de subtítulos que faciliten la lectura del texto. 
Tablas, fotos, figuras 

Incluir fotos y gráficos de buena calidad en jpg o tiff de máxima 
resolución a 1,200 pixeles. 
Referencias bibliográficas 

El formato de las referencias será de acuerdo a la American 
Psychological Association (APA) y deberán aparecer citadas en el 
texto con un número consecutivo entre corchetes. 
FECHAS DE RECEPCIÓN DE ARTÍCULOS 

• Fecha límite de recepción de artículos: 20 de Septiembre de 201 3. 

• Impresión y distribución: Diciembre de 201 3. 

Mayores informes al correo: 

revista.epistemus@correom.uson.mx 




MODELO NO LINEAL CON POTENCIAL SENO-GORDON 
PARA UN CRISTAL FOTÓNICO UNIDIMENSIONAL 

MARGARITA FRANCO ORTÍZ*, ADALBERTO CORELLA MADUEÑO, JUAN ADRIÁN REYES 
CERVANTES, ARNULFO CASTELLANOS MORENO Y RODRIGO ARTURO ROSAS BURGOS 

Se presenta un modelo para un cristal fotónico unidimensional formado por una sucesión deplanosconstituidos 
por material no lineal, inmersos en un medio lineal. Se resuelve la ecuación de Schródinger introduciendo en el 
hamiltoniano modelo del cristal un potencial no lineal tipo Seno-Gordon. Se obtiene analíticamente una versión 
generalizada, dependiente de la amplitud de la onda incidente, de la ecuación trascendente que caracteriza al 
modelo del cristal de Kronig-Penney, de donde se determina la estructura de las bandas de energía del cristal. El 
modelo puede aplicarse a una versión finita de un cristal fotónico para un número limitado de capas alternadas 
de material lineal y no lineal para el cual se han calculado la reflectancia como una función de la intensidad de 
la onda electromagnética, el índice de la banda y el número de períodos. Es posible construir un sistema con 
estas características alternando capas muy delgadas de material de materia blanda no lineal con capas sólidas 
más gruesas, con el cual puede diseñarse un dispositivo para controlar la propagación de luz para intervalos de 
longitudes de onda específicos e intensidades de la luz de la misma señal que se propaga. 
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INTRODUCCIÓN 

Los cristales fotónicos (CFs) son materiales compuestos, 
espacialmente periódicos, que pueden exhibir brechas 
en la banda fotónica (BF) para la propagación de la luz, los 
cuales han llegado a un estado maduro de desarrollo con sus 
propiedades ópticas bien conocidas y con muchas posibles 
aplicaciones (1). La mayor parte de las investigaciones en 
esta área se refieren a CFs cuyas características están fijas, 
es decir, una vez que se han fabricado no hay posibilidad 
de alterar su respuesta óptica. Sin embargo, una tendencia 
reciente se refiere a CFs activos que por medio de algún 
agente externo pueden cambiarse sus propiedades ópticas 
de forma continua y reversiblemente. Esto podría dar 
lugar a guías de ondas ópticas sintonizables, interruptores, 
limitadores y polarizadores, redes ópticas reconfigurables, 
e interconexiones electro-ópticas en la microelectrónica. 

Podemos clasificar a los CFs sintonizables de acuerdo 
a dos grandes categorías: para una de ellas, un agente 
externo causa cambios estructurales sin alterar las 
constantes dieléctricas de los materiales constituyentes; 
en la otra categoría la configuración del CF permanece 
igual y es una propiedad material del CF la que es afectada 
por el agente externo. Se ha realizado la sintonización 
estructural por medio de tensión mecánica aplicada a 
un ópalo polimérico de esferas a escala nanométrica; 
aplicando un campo eléctrico a un cristal fotónico sobre 
un sustrato piezoeléctrico; mediante la incorporación de 
diodos en una estructura 2D de alambres; y aplicando un 
campo magnético a partículas magnéticas distribuidas 
periódicamente (2). Por otro lado, la sintonización a 
través de la alteración de alguna propiedad del material 
involucra la incorporación de algún material ferroeléctrico 
o ferromagnético en el CF para ser sintonizado por un 
campo eléctrico o magnético externo, respectivamente 



(3). En particular, los cristales líquidos (CLs) son materiales 
electro-ópticos bien establecidos que pueden sintonizarse 
por medio de la presión, calor y el campo eléctrico o 
magnético aplicado. Su incorporación dentro de un 
cristal fotónico es de particular interés por la posibilidad 
de sintonizar selectivamente brechas BF, como se revisó 
recientemente (4). 

La ecuación Sine Gordon aparece en 1970 como una 
forma de incluir efectos no lineales en las teorías físicas y su 
uso se extiende desde la teoría del campo relativista hasta 
el estudio de la dinámica de las líneas de transmisión. Su 
interés se centra en la existencia de soluciones llamadas 
solitones, cuyo nombre obedece a que son ondas que se 
propagan sin dispersarse. Existen varias formas de obtener 
la ecuación y entre ellas se pueden mencionar desde 
modelos de cadenas de osciladores hasta la postulación de 
una densidad lagrangiana adecuada, tal que por medio de 
la formulación lagrangiana es posible obtener la ecuación 
correspondiente. Ésta es una labor altamente técnica que 
se ubica al margen del propósito de este trabajo. Nuestro 
objetivo se centra en modelar una estructura periódica 
capaz de intercalar dos materiales: uno cuya respuesta 
a la perturbación ondulatoria es lineal y otra de carácter 
no lineal situada en sitios separados una distancia a. Este 
es el papel que juega la función seno y constituye una 
generalización de otro modelo que hemos estudiado 
anteriormente. 

Los efectos no lineales son muy amplios, pero es 
posible presentar un ejemplo que muestre la diferencia 
entre la conducta lineal y la que no lo es. Si enviamos 
un pulso luminoso a través de un material lineal, la luz 
se propagará de forma que el pulso se va haciendo cada 
vez más ancho. En términos técnicos, se dice que es un 
síntoma de conducta dispersiva. En cambio, en un medio 
no lineal pueden aparecer varios efectos en el pulso que 
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se propaga, Puede ser que el pulso se haga más angosto, o 
bien, que simplemente mantenga su forma. 

ESTRUCTURA DE BANDAS EN EL MODELO 
DE KRONIG-PENNEY LINEAL 

El modelodeKronig-Penney(KP)es un modelo quedescribe 
los estados de energía de un electrón perteneciente a un 
cristal, para ello se supone que la estructura cristalina 
configura un potencial periódico de cambios abruptos que, 
si bien es hipotético, es de gran utilidad para los cálculos. 

El modelo KP es una aproximación al comportamiento 
real de los electrones en el cristal y toma en cuenta la 
periodicidad de la estructura cristalina de los sólidos; 
el potencial periódico real debido a los iones en la red 
cristalina se simplifica del espacio tridimensional (3D) al 
unidimensional (ID) considerando una sucesión periódica 
de pozos y barreras rectangulares unidimensionales, 
donde el centro de cada pozo corresponde a la posición 
que ocupa cada uno de los iones en la red (Figura 1). 



k =— y <i --¡¿(Vo-e) o) 

La función de onda dentro del pozo n + 1: 

Vn+i = A n+1 sen kx n+1 + B n+1 coskx n+1 # < x n+1 < b (4) 

La función de onda y su derivada deben ser continuas 
en el punto A (x n = b; a — x n = c)' 

A n senkb + B n coskb = A' n senhqc + B' n coshqc 
k[A n coskb — B n senkb] = —q[A' n coshqc + B' n sen hqc] (5) 

y también en el punto B ( x n = a; x n + 1 = ): 



2m 




Figura 1 . Modelo de cristal unidimensional (1 D) infinito de 
Kronig-Penney. 

La posición de un punto x en el pozo n es x n , medida 
desde la pared izquierda de este pozo: x n = x - na. 
El punto B tiene la coordenada x n = a, que equivale a 
x n+i = ; a representa la constante de red del cristal, es 
la distancia entre una pareja de iones sucesivos en la red. 

La ecuación de Schródinger correspondiente para este 
sistema está dada por: 

h 2 d 2 ip 

Los potenciales en los pozos son cero v(x) = 0, y en las 
barreras es diferente de cero, es decir, v(x) = V . 

La función de onda electrónica dentro del pozo n es: 

ip n = A n sen kx n + B n eos kx n / < x n < b (i ) 

y dentro de la barrera n: 

xjj n = A' n senhq(a — x n ) + B' n coshq(a — x n ) t b < x n < a (2) 
donde: 



—qA n — kA n+1 



(6) 



Sustituyendo las ecuaciones 6 en las ecuaciones 5 se 
llega a: 

A n senkb + B n coskb = — A n+1 senhqc + B n+1 cosh qc 

A n eos kb — B n sen kb = A n+1 cosh qc — (j^ B n+1 senh qc 

Para simplificar la solución del sistema de 
ecuacionesypasamosallímitel/o -> ^simultáneamente 
las barreras se hacen más y más estrechas, c -> 0, de 
tal manera que el producto cV permanece constante, 
obtenemos así una sucesión periódica de potenciales 
delta repulsivos tipo peine de Dirac (la función peine de 
Dirac S T es una distribución periódica de deltas de Dirac 
espaciadas T y se representa por S T = £-oo S (t — nT) , a 
esta función también se le conoce como tren de impulsos 
unitarios), para los cuales la ecuación de Schródinger 
adopta la forma: 

- — + — S^-.oo S(x - na)ip = Ei/j (8) 
En este límite las ecuaciones 7 se reducen a: 



A n senka + B n eos ka = B n+1 
A n eos ka — B n sen ka = A n+1 — B n+1 



(9) 



donde A es un parámetro adimensional real y positivo 
tal que C V = — . 

U 2771 

Se propone una solución de la forma: 

A n = cJ n y B n = c 2 f n (10) 

Sustituyendo las ecuaciones 10 en las ecuaciones 9 se 
obtiene: 

c ± sen ka + c 2 (cos ka — f) = 
^(cos ka — f) — c 2 {sen ka — = 

Para que existan soluciones con c lt c 2 ^ , se 
requiere que el determinante del sistema se anule. Esta 
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condición conduce a: 

f 2 -2fcosxa + l = (12) 

donde: 

eos x a — cos ka + ^ sen ka (13) 

Las soluciones de la ecuación 1 2 son: 
f = e ±ixa ( 14 ) 

las cuales sustituidas en las ecuaciones 1 1 dan: 
/ — cos ka f — cos ka 



Cl ° 2 sen ka ^ senka ' ^ 2 
Para electrones desplazándose hacia la derecha se 



= A- 



■; A = c 2 



tiene: 

= A f -eos ka QÍxna g Ae i X na 
senka 



(15) 



La sustitución de estos valores en la ecuación 1 da 
como resultado una función de onda de Bloch: 

x/j n = Ae i x x u n (16) 

con: 

^ ^ e l X( a ~ x n) S enkx n +e~ l X x nsen(a-x n ) 
71 senka 

Las ondas de Bloch obtenidas son un caso particular del 
Teorema de Floquet, el cual establece que las ecuaciones 
diferenciales lineales con coeficientes periódicos tienen 
soluciones de la forma e ±ikx u(x) donde u(x) es una 
función periódica de x. 

La energía de los electrones queda determinada por el 
valor de X. Los eingenvalores de la energía se obtienen de 
la ecuación 1 3 para k real (Figura 2). Se debe cumplir que: 



— 1 < cos ka H sen ka < 1 (18) 

ka 
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Figura 2. Gráfica de la ecuación 18, mostrando bandas 
permitidas (sombreadas) separadas por brechas 
prohibidas. 



ESTRUCTURA DE BANDAS EN EL MODELO DE 
KRONIG-PENNEY NO LINEAL 

Ahora se extenderá el estudio para resolver la ecuación 
de Schródinger considerando el caso de un potencial 
periódico no lineal tipo seno-Gordon; podemos escribir la 
ecuación de propagación unidimensional como: 

d 2 ip 
dx 2 



+ k 2 ip- [^Sn=-oo5(x- na) Sen xp^i/j = (19) 



en donde n 2 es el coeficiente de no linealidad y cp es la 
amplitud de la onda de propagación. Para puntos x ^ na 
, en donde n es un número entero, la solución general de 
esta ecuación es de la forma: 

xp = Acoskx + Bsenkx (20) 

Después de una integración estándar de la ecuación 
1 9 alrededor de los puntos x = na , se pueden establecer 
las condiciones de frontera: 

e- if ip(a) = V(0) 



(21) 



= n 2 Sen ib(0) 

dx dx z 

en donde el parámetro de Bloch es f , el cual 
caracteriza el momento en cada celda unitaria. Después de 
aplicar estas condiciones a la ecuación 20 y de eliminar la 
constante B, se obtiene: 



cosf = coska + p 



senA sen ka 
A ka 



(22) 



donde ü = 

p 2 



La ecuación 22 claramente representa una 
generalización de la ecuación 13 que corresponde al 
modelo de cristal lineal de Kronig-Penney. En este modelo 
no lineal, la estructura de bandas depende también de la 
amplitud de la onda incidente. 

RESULTADOS 

En la figura 3 se muestran las bandas permitidas y 
prohibidas que se han encontrado al resolver la ecuación 
22. El eje vertical corresponde al lado derecho de esta 
ecuación, mientras que el eje horizontal representa a ka. 




Figura 3. Gráfica de la ecuación 22: (a) Con A = 1,p = 
0.1,0.2,0.3,0.4,0.5,1,2,3,4,5; (b) p = 2,A = 

0.1,0.2,0.3,0.4,0.5,1,2,3,4,5. 
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La figura 4 muestra la estructura de bandas del CF no 
lineal como una función de las variables independientes 
z — kay A. La figura 5 exhibe la misma información, pero 
en una gráfica de superficies de nivel. 




CONCLUSIONES 

Para amplitudes de la función de onda A constante, 
y coeficiente de no linealidad p pequeño, sólo se tiene 
una banda prohibida en la región de ka cercana a cero, 
conforme p crece, las bandas prohibidas se ensanchan, 
mientras que las bandas permitidas se vuelven cada vez 
más estrechas, hasta quedar localizadas en los puntos nn. 

Para coeficientes de no linealidad p constante, y 
amplitud de la función de onda A variable, la anchura 
de las bandas prohibidas decrece conforme A aumenta, 
mientras que las bandas permitidas se ensanchan, pero 
cuando A alcanza un valor específico inicia un fenómeno 
de oscilación de la anchura tanto de las bandas prohibidas 
como permitidas, tal que a medida que A aumenta las 
bandas prohibidas se tornan cada vez más angostas y las 
bandas permitidas más anchas. Para A suficientemente 
grande desaparecen las bandas prohibidas. 



Figura 4. Gráfica de la ecuación 22: (a) Con p = 100, en 
función de z = ka y A. (b) Gráfica del valor absoluto de la 
ecuación 22 con p = 20, en función de z = ka y A. 



A 



v 




aa 



« z ■ ka 



Figura 5. Superficies de nivel con p = 20, en función de z 
= ka y A. Las regiones oscuras corresponden a las bandas 
permitidas. 
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BÚSQUEDA DE DISCOS PROTOPLANETARIOS EN 
CÚMULOS ABIERTOS JÓVENES 
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LORENZO OLGUÍN RUÍZ Y JESÚS HIGUERA VALDENEBRO 

La formación de planetas está íntimamente ligada a la aparición de discos circunestelares: estructuras aplanadas 
de gas y polvo que se encuentran alrededor de algunas estrellas. Para estudiar el proceso de formación de los discos, su 
frecuencia de aparición y la evolución secular que eventualmente llevará a la formación de planetas, es necesario contar 
con una base observacional amplia que nos permita relacionar las propiedades de los discos con diversas características 
de su estrella central y del medio ambiente donde éstas se forman. En este trabajo presentamos una breve descripción 
del proceso de formación de los discos circunestelares, su posterior evolución y los esfuerzos que hemos realizado un 
grupo de investigadores universitarios en colaboración con colegas nacionales, para obtener resultados observacionales, 
algunos de los cuales se presentan en este artículo, con el fin de contribuir con más y mejores datos al entendimiento de la 
formación planetaria. Hasta ahora, los resultados obtenidos por nuestro grupo son un ingrediente más que ayuda a definir 
y esclarecer el proceso de formación planetaria, particularmente contribuyendo a la identificación de estrellas candidatos 
a poseer un disco circunestelar en cúmulos estelares jóvenes con la finalidad de muestrear el porcentaje de estrellas con 
disco en cúmulos, además del estudio de candidatos individuales. 
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INTRODUCCIÓN 

Desde la antigüedad, el hombre se ha maravillado con la 
observación del cielo nocturno. Su curiosidad lo llevó a 
preguntarse y tratar de comprender la naturaleza de los 
objetos que brillaban en el firmamento y luego, de forma 
natural, le surgió la pregunta. ¿Entre esos objetos celestes 
existen mundos similares al nuestro? Aún en la actualidad 
tenemos sólo una respuesta parcial a esas inquietudes. 
Sin embargo, cada vez comprendemos mejor lo que son 
las estrellas y sabemos - como resultado de modernas 
técnicas de observación astronómica - que alrededor de 
muchas de ellas orbitan objetos de tipo planetario. Al día 
de hoy registramos la existencia de más de 861 planetas 
extrasolares orbitando alrededor de 677 estrellas (128 
tienen más de un planeta) e incluso sabemos que algunos 
de ellos comparten algunas características con nuestro 
propio planeta. Aunque muchos de estos sistemas los 
conocemos únicamente de manera indirecta - gracias al 
efecto que producen en el espectro, el brillo o la posición 
de su estrella-, la presencia de cuerpos planetarios 
extrasolares se ha revelado como un fenómeno común en 
el Universo. 

A medida que ha avanzado el conocimiento en la 
astrofísica y las herramientas con las que le es posible 
trabajar, las observaciones se han vuelto más precisas 
y detalladas. Así, hoy en día no sólo podemos observar 
características de los planetas de nuestro propio sistema 
solar sino de otros sistemas estelares. También es posible 
estudiar los discos que se forman alrededor de las estrellas 
durante su nacimiento, llamados discos circunestelores, 
y que se cree son las estructuras fundamentales para 
la formación planetaria. Para el estudio de estos discos 
primordiales, además de las herramientas observacionales, 
tenemos un mayor avance en las técnicas computacionales 
que nos permiten crear modelos cada vez más precisos 
sobre cómo se da el proceso de formación de sistemas 
similares a nuestro Sistema Solar. Este conocimiento nos 
acerca cada vez más a la respuesta sobre si algunos de esos 
mundos extrasolares son similares al nuestro e, incluso, a 
conocer si pudieran llegar a sustentar alguna clase de vida. 

Las primeras detecciones de discos circunestelares 
que mostraban la presencia de un planeta, se hicieron en 
luz visible con el conocido caso de la estrella Beto Pictoris 
(1) (Figura 1). Desde entonces y hasta la actualidad, con 
técnicas de observación cada vez mejores, disponemos de 
imágenes de discos circunestelares jóvenes que van desde 
los formados principalmente por gas con una componente 
minoritaria de polvo - como los discos mostrados en 
la figura 2 pertenecientes a la región de Orion - hasta 
imágenes de discos más evolucionados, compuestos 
principalmente por material sólido y apenas rastros de gas 
(discos "de residuos"), como el caso de la estrella Fomalhaut 
(Figura 3). Estudiar la evolución de estos discos nos llevará 
a establecer una imagen completa del proceso físico de 
formación de los sistemas planetarios que actualmente 
son comunes alrededor de las estrellas. 




Figura 1. La imagen muestra un disco de polvo y gas 
alrededor de la estrella Beto Pictoris a unos 50 años luz de 
la Tierra. El disco se muestra de canto en colores falsos y se 
extiende por más de 400 UA (1 UA - unidad astronómica 
- corresponde a 149 597 871 km). Esta histórica imagen 
fue la primera confirmación en el óptico de la presencia 
de discos de gas alrededor de estrellas jóvenes. Crédito: 
Carnegie Institution, procesada por la University of Arizona 
y el Jet Propulsión Laboratory (JPL), cortesía de NASA. 
Imagen del telescopio de 2.5 m del Observatorio Las 
Campanas en Chile. 




Figura 2. Imágenes del Hubble Spoce Telescope (HST) de 
discos protoplanetarios (proplyds) en Orion. Las imágenes 
fueron obtenidas con un filtro de luz Ha. El disco aparece 
en absorción para la mayoría de las líneas ópticas, p. ej. HST 
294-606 contiene un disco visto de lado de 
aproximadamente 400 UA de radio. En todos los casos 
puede observarse claramente el frente de ionización que 
corresponde a la zona más brillante de cada imagen. 
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Figura 3. Disco de residuos que rodea al sistema Fomalhaut. 
Estos discos son la etapa final de los discos protoplanetarios 
donde el disco ya se ha vuelto ópticamente delgado 
(menos opaco) y contiene anillos de material muy 
evidentes, además de la presencia de cuerpos con masa 
planetaria y muy poco gas. Ésta es una de las etapas finales 
de la evolución de los discos protoplanetarios y la 
culminación de la formación de un sistema planetario. 
Imagen: NASA/ESA/P. Kalas (Univ. of California, Berkeley) 
etal. 

FORMACIÓN ESTELAR Y LA APARICIÓN DE LOS 
DISCOS CIRCUNESTELARES 

El estudio de discos circunestelares comienza con la 
formación misma de las estrellas. Actualmente, la teoría de 
formación estelar más aceptada nos dice que las estrellas 
se forman a partir del colapso gravitacional que ocurre 
en enormes nubes de gas molecular que se encuentran 
en el plano de nuestra Galaxia (2). Dentro de estas 
extensas nubes se forman núcleos con mayor densidad 
que ocasionalmente son afectados por una perturbación 
externa (colisiones entre nubes, explosiones de supernovas, 
turbulencia del medio interestelar, etcétera) y responden a 
ésta colapsando sobre sí mismos por acción de su propia 
gravedad. Durante el colapso, los núcleos se fragmentan 
dando lugar a una gran cantidad de protoestrellas, es decir, 
estrellas en proceso de formación. Debido a que la nube 
tenía un movimiento de rotación inicial, la conservación 
del momento angular causa que a medida que los núcleos 
colapsan para formar protoestrellas, se forme alrededor 
un disco aplanado de material, perpendicular al eje de 
rotación de la protoestrella y que también se encuentre en 
rotación. Finalmente, la protoestrella termina con un disco 
en rotación a su alrededor llamado disco circunestelar. 
De acuerdo a las teorías más aceptadas actualmente, 
este disco atraviesa por diversas etapas de evolución que 
ocasionalmente conducen a la formación de planetas. Así 
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pues, es la rotación de la nube (a través de la conservación 
del momento angular) el fenómeno físico que causa 
la formación de discos de gas y polvo alrededor de las 
estrellas y será en gran medida la rotación del mismo la 
que determine su posterior evolución. Es importante 
resaltar que la formación de un sistema planetario surge 
como subproducto de la formación de una estrella y 
aunque ocurre con gran frecuencia, no todas las estrellas 
tienen planetas. 




Los discos primordiales están embebidos en la nube 
de polvo y gas molecular en colapso, por lo que es difícil 
observarlos en el intervalo de luz óptica. Sin embargo, 
pueden ser detectados en longitudes de onda infrarroja y 
milimétrica. Es posible estudiar observacionalmente esta 
primera fase en estrellas muy jóvenes de baja masa tipo 
T-Taur¡, al realizar un análisis minucioso de la distribución 
espectral de energía del sistema. Particularmente, el exceso 
de emisión en el infrarrojo cercano y medio es característico 
de la presencia de un disco, pues ésta es producida por 
el polvo presente en él. Con base en esta emisión es 
posible determinar la morfología del disco y determinar 
la temperatura y cantidad de polvo en él, así como sus 
características químicas (3). Además, es posible observar 
la emisión infrarroja de la envolvente más externa todavía 
en colapso esférico. En la evolución posterior, el polvo del 
disco sufre coalescencia (4), formando primero cuerpos de 
tamaño de metros y posteriormente alcanzando tamaños 
de unos pocos kilómetros, llamados planetesimales, que 
son cuerpos del orden de 10 km o más. Eventualmente, 
las colisiones entre planetesimales permitirá la formación 
de cuerpos mayores (5). En esta etapa, los discos son 
llamados discos protoplonetorios. También puede suceder 
que las colisiones sean tan violentas que el efecto sea la 
destrucción de los planetesimales, dejando una gran 
cantidad de polvo y cuerpos menores (2). 

En el punto final de la evolución, después de la 
formación de un sistema planetario, la distribución 
espectral de energía de la protoestrella se vería como se 
muestra en la figura 4 (objeto Clase III). Ahí se muestra el 
espectro de un cuerpo negro que representa la emisión 
fotosférica de la protoestrella central, ya muy cerca de su 
entrada a la secuencia principal, y por lo tanto a punto 
de convertirse en una verdadera estrella. La presencia 



del disco en luz infrarroja ha desaparecido, posiblemente 
porque el disco ha dado lugar a un sistema planetario. 
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Figura 4. Diagrama esquematizado del proceso de 
formación de sistemas planetarios. Posteriormente a 
la formación inicial del disco, éste evoluciona como 
resultado del transporte de momento angular, así como 
por su interacción con la radiación de la protoestrella 
central. El polvo presente en el disco forma agregados 
cada vez mayores, dando lugar a planetesimales y 
posteriormente, planetas. Las fases intermedias de este 
proceso de formación de planetas aún se encuentran en 
constante análisis y discusión, al igual que las fases finales 
de disipación total del disco (6). 



El proceso de formación planetaria dura del orden de 
pocos millones de años como puede verse en la figura 
5, donde se muestra el porcentaje de estrellas con disco 
en cúmulos estelares jóvenes como función de la edad. 
Puede observarse que a mayor edad del cúmulo estelar, 
la frecuencia de discos entre sus miembros disminuye. De 
esta forma, cúmulos con edades entre 5 y 10 millones de 
años tienen un porcentaje prácticamente cero, es decir, 
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que los discos han desaparecido. Esto puede deberse 
a la evaporación causada por los vientos estelares, a la 
migración dinámica, a la formación de planetas o a una 
combinación de todos ellos. 
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Figura 5. Diagrama de la frecuencia de discos en cada 
cúmulo como función de la edad del cúmulo, ésta se 
presenta en millones de años (Ma). Puede observarse que 
los cúmulos más jóvenes muestran un alto porcentaje de 
estrellas con disco y a medida que el cúmulo es más viejo, 
este porcentaje decrece hasta llegar a valores cercanos a O 
alrededor de los 1 Ma de edad (7). 



DISCOS PROTOPLANETARIOS 
Y CÚMULOS ESTELARES 

Debido a que la formación y evolución de los discos 
protoplanetarios ocurre simultáneamente en las primeras 
etapas de la formación estelar, es evidente pensar 
que el estudio de los sistemas estelares más jóvenes 
nos dará información sobre la ocurrencia, frecuencia, 
comportamiento y desarrollo de éstos. En nuestra Galaxia, 
la formación estelar ocurre en las regiones del plano 
galáctico donde se encuentran las nubes moleculares, 
las cuales dan lugar a agrupaciones estelares de cientos o 
miles de estrellas, llamadas cúmulos abiertos o galácticos. 
Para la mayoría de estos cúmulos es posible calcular su 
edad, lo que nos permite seleccionar sólo aquellos que 
sean realmente jóvenes, astronómicamente hablando, 
con edades menores a 10 millones de años, para realizar 
en ellos una búsqueda de discos protoplanetarios. En 
cúmulos con edades mayores a ese límite, es probable que 
los discos ya hayan evolucionado a un sistema planetario y 
que ya se haya dado la dispersión casi total del gas y polvo 
que dejaron como residuo. 

Ya que los cúmulos abiertos cuentan con una 
gran cantidad de estrellas, es posible hacer un estudio 



estadístico de la frecuencia de los discos y la ocurrencia de 
éstos en cada una de sus diferentes etapas. Además, debido 
a que es posible determinar su edad, podemos obtener 
una buena estimación de los tiempos de cada una de las 
fases de evolución de los discos. Por ejemplo, si el tiempo 
de vida de un disco protoplanetario es de 5 a 7 millones 
de años y observamos cúmulos entre 7 y 10 millones 
de años, será muy poco probable encontrar emisión de 
algún tipo de disco en sus estrellas. Estudios de este tipo 
nos permiten conocer cuál es la frecuencia observada de 
discos protoplanetarios, estudiar sus etapas de evolución y 
hacer predicciones de la cantidad de sistemas planetarios 
que se encuentran en las estrellas. 

Actualmente, nuestro grupo de trabajo realiza el 
estudio de cúmulos estelares jóvenes, con edades menores 
a 10 millones de años, buscando la emisión proveniente 
de estrellas con disco. Como primer paso, es conveniente 
detectar cuántas y cuáles estrellas de los cúmulos en 
estudio son candidatas a poseer un disco protoplanetario. 
Para esto se construyen diagramas fotométricos, en los 
cuales es posible detectar estrellas que presentan un 
exceso de emisión infrarroja. 

En la figura 6 se presenta el diagrama color-color 
del cúmulo abierto NGC 2414, realizado con datos del 
catálogo público 2MASS (Two Micron All-Sky Survey). 
Este catálogo contiene información de la emisión en el 
infrarrojo cercano (alrededor de 2 micrómetros) de todas 
las estrellas y otros objetos celestes, hasta cierto límite de 
brillo. En nuestro estudio de NGC 241 4, los datos de 2MASS 
junto con resultados de modelos teóricos de la emisión de 
radiación de los discos, nos ha permitido detectar aquellas 
estrellas que son fuertes candidatas a poseer un disco 
protoplanetario. Una vez detectadas, es posible construir 
la distribución espectral de energía de cada candidata 
con base en las observaciones infrarrojas en conjunto 
con observaciones ópticas obtenidas en el Observatorio 
Astronómico Nacional en San Pedro Mártir, B.C. (OAN-SPM, 
IA-UNAM), con el fin de confirmar la presencia de un disco. 
De confirmarse, el disco puede ser catalogado dentro de 
alguna de las fases de evolución que ya se han discutido. 
Esto nos dará una ¡dea de la frecuencia de ocurrencia de 
discos como función de la edad del cúmulo, así como 
del estado evolutivo en que se encuentran. Asimismo, es 
posible estudiar la correlación entre las características de 
la estrella y el disco que la rodea. Todo con la finalidad de 
comprender tanto las propiedades físicas de los discos 
protoplanetarios (tamaño, masa, densidad, temperatura), 
como su evolución y los procesos físicos que dan lugar a su 
transformación en sistemas planetarios. 

A pesar de lo mucho que se ha avanzado hasta ahora 
en el conocimiento de las características y la formación 
de sistemas planetarios, éste sigue siendo un tema de 
investigación en constante desarrollo. Debido a la gran 
diversidad de propiedades y configuraciones que muestran 
los planetas extrasolares encontrados hasta el día de hoy, 
ninguno de los mecanismos propuestos para su formación 
puede catalogarse como universal. Sin embargo, 
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continuar con el desarrollo de estos estudios nos llevará 
a un mejor entendimiento de la evolución de los discos 
protoplanetarios y la subsecuente formación de planetas. 
De esta manera, nos acercamos más a la respuesta acerca 
de nuestro propio origen y nuestra posterior evolución en 
el contexto astronómico. 
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Figura 6. Diagrama color-color de fotometría de 2MASS del 
cúmulo abierto NGC 2414. Los filtros J (1.24 jim), H (1.65 
|jm) y K (2.16 |jm) permiten el paso de la luz de un objeto 
en bandas del infrarrojo y el diagrama de los colores 
derivados, (J-H) vs (H-K), nos permite encontrar aquellos 
objetos con exceso de emisión en estas bandas. En esta 
figura, los círculos indican las estrellas que muestran un 
exceso en infrarrojo y que fueron seleccionadas como 
objetos candidatos a poseer un disco de acreción (8). 
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MATERIALES ULTRABLANDOS: DIAGRAMA DE FASE DE 
UNA SUSPENSIÓN COLOIDAL DE POLÍMEROS ESTRELLA 

GREGORIO ALEJANDRO VEJAR DELGADO*, 
HERIBERTO ACUÑA CAMPA Y LAURA L YEOMANS REYNA 

La determinación experimental y predicción teórica de diagrama de fases para las sustancias es un 
aspecto de la mayor importancia en la ciencia de materiales. La cristalización y fusión a pesar de ser 
fenómenos cotidianos y de estar dentro de la agenda de estudio de la mecánica estadística desde sus 
inicios, aún no provee una teoría que las explique de primeros principios. La dificultad principal se debe 
a que estas transiciones se presentan en sistemas concentrados convirtiéndose en un problema colectivo 
de muchos cuerpos. Es en este contexto que la prescripción de criterios fenomenológicos que permitan 
la localización de las líneas de transición es altamente valorada por la comunidad científica. Motivados 
por ello, en esta comunicación se presentan resultados obtenidos con el criterio de Lówen, mediante 
simulaciones computacionales, para una suspensión coloidal de partículas ultra suaves conformada por 
polímeros estrella. 
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INTRODUCCIÓN 

A pesar de que la transición sólido-líquido es no universal, 
existen algunos criterios fenomenológicos basados en 
propiedades de una de las fases en coexistencia. Las 
ventajas de estos criterios consisten en que permiten 
estimar la línea de coexistencia sólido-líquido sin necesidad 
de realizar cálculos de energías libres. Además, predicen 
valores cuasi-universales de propiedades estáticas y 
dinámicas. Por valores cuasi-universales se entiende 
que los valores de la propiedad estática o dinámica que 
caracteriza la transición no es propiamente una constante; 
sin embargo, es un número cercanamente fijo en alrededor 
deun10%(1). 

Se debe distinguir entre dos clases de universalidades. 
Primero, se puede variar la temperatura de un sistema 
manteniendo fijo el potencial par entre las partículas y 
estudiar las propiedades a lo largo de la línea de transición. 
Haciendo uso de escalamiento de propiedades de 
potenciales de interacción de potencia inversa, por ejemplo 
l/r n , se observa de inmediato que cualquier cantidad 
adimensional adecuadamente escalada es universal a lo 
largo de la línea de fusión y cristalización. Esta primera 
universalidad es trivial para potenciales de potencia 
inversa, sin embargo es menos trivial para potenciales más 
complicados. El segundo tipo de universalidad, que es más 
general, se manifiesta si una propiedad tiene el mismo 
valor en la transición sólido-fluido para diferentes sistemas, 
es decir, para sistemas cuyas partículas interaccionan con 
diferentes tipos de potenciales pares. 

A continuación se describen brevemente tres de los 
criterios fenomenológicos más utilizados en la descripción 
de las transiciones sólido-fluido: criterio de Lindemonn, 
criterio de Hansen-Verlet y finalmente el criterio de Lówen. 

Criterio de Lindemann: Este criterio fue la primera regla 
fenomenológica para describir la transición de fusión 
(sólido-fluido), propuesta por Lindemann en 1910 en su 
artículo The colculotion of molecular vibrotion frequencies 
(2). El criterio establece que el cociente L de la raíz cuadrada 
del desplazamiento cuadrático medio y la distancia 
promedio entre partículas adquiere un valor universal en 
la línea de transición: 



a 

La pregunta sobre el carácter de la universalidad de 
este criterio, particularmente, sobre si es del segundo 
tipo, es decir, válido independientemente del potencial 
de interacción entre las partículas del sistema, fue 
intensamente discutido durante las décadas de 1980-1990. 
Mediante cálculos con simulaciones computacionales 
sabemos ahora que el valor de i depende ligeramente del 
tipo de potencial de interacción y su valor aproximado 
es de 0.15 (L=0.185 para potenciales muy suaves; 1-0.17 
para potenciales de partículas cargadas tipo Yukawa; 
1=0.129 para esferas duras). El criterio de Lindemann se 
puede entender intuitivamente como el rompimiento del 



ordenamiento de la celda formada por partículas vecinas 
si su desplazamiento excede un valor crítico, es decir, si 
éste adquiere un valor un orden de magnitud menor que 
el correspondiente a la distancia media de sus partículas 
vecinas. 




Criterio de Hansen-Verlet Este criterio constituye una 
regla fenomenológica para la transición de cristalización, 
misma que fue formulada en 1969 por Hansen y Verlet en 
su artículo Phase Transitions of the Lennard-Jones System 
(3). Mediante simulaciones por computadora encuentran 
que para un líquido de Lennard-Jones el valor máximo 
de la amplitud del factor de estructura adquiere un valor 
constante a lo largo de la línea de transición convirtiéndose 
con ello en un criterio de universalidad del primer tipo: 

S(/W)~2.85 (2) 

Sin embargo, mediante experimentos de dispersión 
y cálculos de simulaciones moleculares para otros tipos 
de sistemas cuyos potenciales de interacción no son de 
Lennard-Jones, se ha confirmado que la transición de 
cristalización se presenta siempre cuando el valor máximo 
de S(k) es cercano al valor de 3. De esta forma se considera 
que el criterio de Hansen y Verlet también posee la 
universalidad del segundo tipo, es decir, tiene un sentido 
más general. 

Criterio de Lówen-Palberg-Simon: En 1993, Lówen, 
Palberg y Simón proponen un criterio estrictamente 
dinámico para describir la transición de cristalización para 
sistemas coloidales. En su artículo Dynamical Criterion for 
Freezing of Colloidal Liquids (4) aportan una prescripción 
sencilla que establece que el cociente del coeficiente de 
difusión de tiempos largos D L y el coeficiente de difusión 
de tiempos cortos D s adquiere un valor constante a lo largo 
de la línea de transición de cristalización: 

D L 

-~0.1 ,3, 
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Este criterio fue planteado como producto del 
estudio de la cristalización de fluidos coloidales cargados 
mediante simulaciones de dinámica Browniana para un 
sistema de partículas que interaccionan con un potencial 
de interacción modelo tipo Yukawa, así mismo, presentan 
resultados para un fluido de esferas duras y para plasmas 
de una componente (OCP). En el mismo artículo reportan 
resultados experimentales de dispersión de Rayleigh 
de suspensiones coloidales estabilizadas con carga que 
confirman esta regla fenomenológica. De lo anterior 
se desprende que el criterio dinámico es universal en 
el sentido más general (segundo tipo) por cumplirse 
para diferentes tipos de potenciales de interacción. Más 
aun, en 1996 el propio Lówen reporta que el criterio 
dinámico propuesto y verificado para el caso de sistemas 
tridimensionales, también se satisface para describir 
transiciones de cristalización en sistemas de partículas en 
espacio bidimensional, que interaccionan mediante un 
potencial que depende del inverso de una potencia de la 
separación entre las partículas. 

En los últimos años se ha explorado el alcance de estos 
criterios fenomenológicos en condiciones diferentes, 
por ejemplo, Hoffmann y Lówen estudian mediante 
simulaciones de no-equilibrio a un sistema coloidal de 
partículas cargadas en presencia de un campo eléctrico 
oscilatorio, obteniendo que en estado estacionario los 
criterios fenomenológicos de Lindemann, Hansen-Verlet y 
de Lówen-Palberg-Simon también se satisfacen en sistemas 
ligeramente polidispersos (5). Por otra parte, en años 
recientes y producto de la interacción entre académicos 
a través de las Redes Temáticas de Colaboración, como 
la Red Física de la Materia Blanda (PROMEP) y la Red 
de la Materia Condensada Blanda (CONACYT), se ha 
podido avanzar en el estudio de la equivalencia dinámica 
entre fluidos atómicos y coloidales. Una consecuencia 
importante de dicha equivalencia ha sido la extensión del 
criterio dinámico de Lówen, planteado originalmente para 
sistemas coloidales, al caso de la transición de cristalización 
en sistemas atómicos. Este resultado fue corroborado 
mediante la comparación de simulaciones de dinámica 
molecular y dinámica Browniana para sistemas de esferas 
duras y para otros sistemas de partículas con potenciales 
repulsivos suaves (Lennard-Jones truncado y Yukawa) que 
son dinámicamente equivalentes a un fluido de esferas 
duras (6). 

SISTEMA MODELO DE POLÍMEROS ESTRELLA 

Los polímeros estrella son macromoléculas compuestas 
por f cadenas poliméricas ancladas químicamente a un 
centro común. Al número de brazos f conectados se le 
llama funcionalidad de la estrella y cuando la longitud 
característica de las cadenas es suficientemente grande, la 
longitud del tamaño del núcleo central es irrelevante. En la 
figura 1 se ¡lustran configuraciones de polímeros estrella 
para dos funcionalidades diferentes (7). 





Figura 1 . Polímeros estrella obtenidos con simulación para 
(a)f=70y(b)f=50. 

La síntesis de los primeros polímeros estrella se hizo 
en la década de 1990 y como lo señala Likos, desde el 
punto de vista experimental la síntesis de estructuras 
regulares (en funcionalidad y longitud de las cadenas) 
hizo posible explorar la física de sistemas monodispersos 
ultrasuaves. Desde el punto de vista de aplicación 
tecnológica, los polímeros estrella han sido importantes 
como modificadores de la viscosidad en la industria de 
hidrocarburos, así como también como materiales de 
recubrimiento de superficies, entre otras. Los polímeros 
estrella constituyen un sistema importante dentro del 
campo de la materia condensada blanda, toda vez que en 
su estudio confluyen dos grandes campos de ella: la física 
de polímeros y la física coloidal (8). 

Para funcionalidades pequeñas (f=1y2) los polímeros 
estrella se asemejan a una estructura lineal, mientras que 
para valores grandes de la funcionalidad adquieren una 
forma esférica y de consistencia más rígida, de forma tal 
que en el límite de funcionalidad extrema se aproximan al 
comportamiento de un sistema de esferas duras. 

El modelo de potencial de interacción para polímeros 
estrella tiene una forma tipo Yukawa para r ><j y una 
forma logarítmica para r < a, mismas que se unen 
adecuadamente en r = o pidiendo continuidad en 
sus primeras derivadas. La longitud característica a se 
identifica como el diámetro de corona dentro del modelo 
de Daoud-Cotton. La dependencia logarítmica del modelo 
fue derivada por Witten y Pincus en 1986, mientras que el 
decaimiento tipo Yukawa se incluyó de forma empírica 
para ajustar a los resultados experimentales obtenidos 
con dispersión de neutrones de bajo ángulo (SANS). De 
esta forma, la expresión matemática del potencial de 
interacción efectivo por pares se expresa como: 



- in + 



1 + 



si r < g 



1+ f 



exp 



2 V 



(4) 



si r > o 



donde k B es la constante de Boltzmann y T la temperatura 
absoluta. En la figura 2 se ¡lustra el comportamiento del 
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modelo de potencial para 18 < / < 256. Como se observa, 
conforme se incrementa el valor de la funcionalidad se 
presenta más abrupta la transición para r < a, consistente 
con el hecho de que a funcionalidades suficientemente 
grandes el potencial de interacción debe corresponder al 
potencial de esferas duras. 




I r. 

Figura 2. Potencial de interacción entre polímeros estrella. 
La curva del extremo izquierdo corresponde a f = 32 
mientras que la del extremo derecho a f=256. 

CONSTRUCCIÓN DEL DIAGRAMA DE FASE 

En 1999 Watzlawek y colaboradores presentan el diagrama 
de fases (Figura 3) de un sistema monodisperso de 
polímeros estrella, obtenido mediante el método del 
cálculo de energías libres e integración termodinámica a 
partir de cálculos de Montecarlo (9). 




7 ¿ j 

Figura 3. Diagrama de fases obtenido por Watzlawek, Likos 
y Lówen. 



Como se observa en la figura 3, para fracciones en 
volumen ( rj = np*/6 , con p* = per 3 la densidad reducida) 
menores a 0.7 y funcionalidades mayores a f=34 se 
predicen estructuras cristalinas fec y bec (cúbico centrado 
en la cara y cúbica centrada en el cuerpo, respectivamente); 
mientras que para fracciones en volumen mayores a 0.7 y 



funcionalidades mayores a 44, establecen la existencia de 
estructuras cristalinas más complejas tipo bco y diamante. 
Finalmente, se predicen fases fluidas tanto para valores 
pequeños de la fracción en volumen (situación normal 
en sistemas de partículas duras) y también para valores 
grandes, siempre y cuando en este segundo caso las 
funcionalidades de los polímeros estrellas asuman valores 
adecuadamente pequeños de la funcionalidad. Con lo 
anterior el diagrama de fases presenta diferentes tipos de 
re-entrancias fluido-sólido y también entre fases cristalinas 
diferentes. Por otra parte, en 2007, Laurati y colaboradores 
(10), obtienen resultados experimentales de SANS para un 
sistema micelar PEP-PEO de polímeros estrella, en la región 
de funcionalidades intermedias y fracciones en volumen 
rj < 0.5. 

Como se señaló previamente, en el contexto del 
estudio de las propiedades estructurales y dinámicas 
de sistemas coloidales de partículas ultrasuaves, nos 
planteamos explorar la construcción de un diagrama 
de fases para el sistema de polímeros estrella a partir de 
la aplicación del criterio fenomenología) de Lówen y 
mediante simulaciones de dinámica Browniana. Para ello 
se elaboró e implemento un código en lenguaje FORTRAN 
haciendo uso del algoritmo de Ermak y McCammon, que 
permite calcular la posición como función del tiempo de 
N partículas de una suspensión coloidal en el régimen 
difusivo a partir de las Ecuaciones de Langevin. 

Conocidas las configuraciones del sistema (posiciones 
de todas las partículas a tiempos diferentes) es posible 
calcular las propiedades estadísticas dinámicas de 
la suspensión, como son el coeficiente de difusión 
dependiente del tiempo D(t) y el desplazamiento 
cuadrático medio W(t). Al igual que con simulaciones de 
Montecarlo, también es posible calcular propiedades 
estructurales como la función de distribución radial g(r). 
En la figura 4 se ¡lustra la concordancia de los resultados 
obtenidos con estos dos métodos de simulación, 
mediante la comparación de las funciones de distribución 
radial obtenidas para un sistema de polímeros estrellas de 
baja funcionalidad (partículas muy blandas) y un amplio 
intervalo de valores de la fracción en volumen <t>. 
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Figura 4. Función de distribución radial de suspensiones 
de polímeros estrella obtenidos mediante simulaciones de 
Montecarlo y dinámica Browniana. 
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Para implementar el criterio de Lówen fue necesario 
calcular coeficientes de difusión como función del tiempo, 
de forma tal que fuese posible identificar si en un tiempo 
suficientemente largo, el cociente D L /D adquiría un valor 
aproximado a 0.1. Es decir, buscamos la isodifusiva de 0.1 
para una suspensión coloidal monodispersa de polímeros 
estrella. Ello implicó realizar simulaciones sistemáticas: 
(i) para valores fijos de la funcionalidad aumentando la 
fracción en volumen desde las regiones fluidas de baja 
concentración; (¡i) para las regiones de baja funcionalidad 
y alta concentración, fue necesario realizar simulaciones 
sistemáticas para valores fijos de fracción en volumen 
aumentando el valor de la funcionalidad desde las regiones 
fluidas de sistemas ultra blandos (funcionalidad pequeña). 
En la figura 5 se ¡lustra una gráfica de exploración de la 
transición de cristalización. 




Figura 5. Coeficiente de difusión dependiente del tiempo. 
Aplicación del Criterio de Lówen en un sistema de fracción 
en volumen fija y funcionalidad variable. 

En la figura 6 se muestran los resultados obtenidos 
haciendo uso del criterio de Lówen (círculos rojos) 
con las simulaciones de dinámica Browniana descritas 



previamente; así mismo y de forma comparativa, se 
incluyen los resultados obtenidos mediante los cálculos 
de energía libre e integración termodinámica de la figura 
3 (círculos azules). Con símbolos (cruces rojas para fluido 
y negras para cristal) se han colocado los resultados 
experimentales para el sistema de polímeros estrella de 
micelas PEP-PEO obtenidos por Laurati y colaboradores 
con la técnica de SANS. En la figura también se incluyen 
resultados teóricos obtenidos del cálculo de energías 
libres a partir de propiedades estructurales obtenidas con 
teorías de líquidos (Roger-Young) y para la región sólida 
mediante teoría de perturbaciones. 
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Figura 6. Diagrama de fase de un sistema monodisperso de 
polímeros estrella. 

Como se concluye de la figura comparativa, en el 
régimen de alta funcionalidad y baja fracción en volumen, 
los resultados obtenidos con el criterio de Lówen coinciden 
con los reportados haciendo uso de cálculos de energía 
libre con Montecarlo y los resultados experimentales 
reportados en la literatura. Esto nos brinda confianza en 
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la metodología seguida para la exploración de la línea 
de transición de cristalización mediante el cálculo de 
coeficientes de difusión dependientes del tiempo con 
simulaciones de dinámica Browniana. En la región de 
partículas extremadamente suaves (funcionalidad baja) 
y fracciones en volumen intermedias y altas (presencia 
de traslapes entre los polímeros estrella), los resultados 
obtenidos con el criterio de Lówen aparentan estar 
comprendidos entre los resultados reportados con cálculos 
de energía libre con Montecarlo y teóricos. 

CONCLUSIONES 

En este artículo se presenta una descripción de los 
criterios fenomenológicos más valiosos en el estudio 
de las transiciones de fase sólido-líquido en el campo 
de la materia condensada. Dos de ellos corresponden a 
criterios de carácter estructural (Lindemann y Hansen- 
Verlet) y el tercero (Lowen-Palberg-Simon) de carácter 
estrictamente dinámico. Mediante simulaciones de 
Montecarlo y dinámica Browniana se estudian las 
propiedades estructurales y dinámicas de una suspensión 
coloidal monodispersa de partículas ultrablandas de 
polímeros estrella para diferentes condiciones de suavidad 
(funcionalidad) y concentración (o fracción en volumen). 
Como aplicación relevante, se usa el criterio dinámico de 
Lówen para construir el diagrama de fase de la suspensión, 
concretamente para obtener la línea de transición de 
cristalización. 

Por los resultados obtenidos, se puede dar cuenta 
del valor de la herramienta de la simulación molecular 
para entender, extender y comparar con resultados 
teóricos y experimentales de sistemas fisicoquímicos de 
estructura compleja. En particular para el sistema formado 
por polímeros estrella es relevante la comparación de 
los resultados mostrados en el diagrama de fase, ya que 
manifiesta, por una parte, la validación del modelo de 
interacción involucrado en la simulación, así como la 
utilidad de los criterios fenomenológicos, de Lówen en 
este caso, para estudiar comportamientos en líneas de 
transición fase. 

Para el futuro se abre la perspectiva de estudiar estos 
sistema ultrasuaves (o ultrablandos), entre otros, para el 
caso de mezclas. Particularmente nos interesa explorar a 
profundidad el caso de sistemas polidispersos en tamaño 
y funcionalidad. Estudios preliminares nos indican que 
para polidispersidad en tamaño de aproximadamente 
12% el diagrama de fases no se modifica sustancialmente 
para fracciones en volumen menores a 0.7, sin embargo 
para valores mayores, el sistema se convierte en uno 
altamente interactuante que requiere de una exploración 
más detallada y una revisión técnica más minuciosa, sobre 
las condiciones de continuidad del potencial demandadas 
en las simulaciones de dinámica Browniana. Consideramos 
que la incorporación del factor de polidispersidad (tanto 
en funcionalidad como en tamaño) para responder a 
preguntas tales como ¿Cuál es la modificación que se 



induce por polidispersidad en el diagrama de fases? 
o ¿Es posible frustrar por completo una transición de 
equilibrio (como la de cristalización) para obtener estados 
metaestables como los vidrios? es de la mayor relevancia, 
toda vez que en las suspensiones coloidales reales, la 
polidispersidad es una propiedad inherente a los sistemas 
y es todo un reto el control experimenal de esta propiedad. 
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MODELO PREDICTIVO DE RIESGO DE MOROSIDAD PARA 
CRÉDITOS BANCARIOS USANDO DATOS SIMULADOS 
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Un problema grave en las instituciones de crédito son los préstamos no recuperados, por lo que al 
momento de la solicitud, es deseable un modelo matemático para predecir el riesgo de tener un cliente 
potencialmente moroso, ya que su capacidad económica no garantiza por si sola el pago. Este trabajo 
presenta, a través de datos simulados, una técnica estadística conocida como regresión logística para 
predecir la probabilidad de no recuperar un crédito, con base al perfil socioeconómico del cliente, 
considerando el número de dependientes económicos, número de impagos previos, su empleo, destino 
del crédito y salario mensual. El modelo predictor logra una tasa aciertos de 95% para clientes morosos y 
una tasa global del 70% de clasificaciones correctas. Los resultados sirven para apreciar los indicadores 
estadísticos conocidos, como la razón de ventajas por categorías. Por ejemplo, por cada crédito previo no 
pagado del cliente, se estima que el riesgo es 3.7 veces mayor y por cada $ 1,000.00 de incremento en el 
salario se estima que el riesgo (multiplicativo) disminuye por un factor del 80%. 
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INTRODUCCIÓN 

La regresión logística ha sido aceptada durante las 
últimas décadas como modelo de predicción-clasificación 
desde sus orígenes en investigaciones epidemiológicas y 
clínicas (1). Hoy en día se emplea comúnmente, aunque 
no limitada a, investigaciones biomédicas, negocios y 
finanzas, criminología, ingeniería, salud pública, política, 
bilogía de la vida salvaje y psicología. Las aplicaciones van 
desde créditos no recuperados (2), votación en elecciones 
políticas (5), asignación de una beca, hasta predicción 
de cáncer de próstata (4), riesgo de contagio de VIH (4), 
etcétera. 

La regresión logística es una técnica de estadística 
multivariada basada en los principios de regresión lineal, 
pero a diferencia de ésta, la regresión logística utiliza una 
variable dependiente categórica binaria o dicotómica (en 
vez de una variable cuantitativa continua) cuyo resultado 
o valor de respuesta es o 1; por ejemplo, en créditos 
bancarios se asigna 1 si los datos observados provienen de 
un cliente moroso y en caso contrario. Debido a que la 
morosidad es una variable dicotómica que depende como 
variable de respuesta de otros indicadores tales como el 
salario y las cargas familiares, entre otros, es factible aplicar 
la regresión logística como modelo clasificador. Además, 
el modelo puede utilizarse para estimar la probabilidad de 
que un cliente no devuelva el crédito; a dicha probabilidad 
se le conoce como riesgo de morosidad. El status de 
morosidad se categorizará dependiendo si la probabilidad 
estimada, evaluada con los datos del cliente, es mayor 
que un punto de corte, el cual debe prefijarse por cada 
institución particular con base en su experiencia, por 
ejemplo: Pe = 0.30, o bien usar una técnica formal para 
seleccionarlo (3). 

EL MODELO DE REGRESIÓN LOGÍSTICA 

El modelo de regresión logística surge cuando queremos 
estimar la probabilidad de un evento dicotómico de 
Si (1) o No (0), en función de un conjunto de variables 
predictoras comúnmente llamados factores de riesgo, que 
pueden ser discretas o continuas, categóricas (nominales 
u ordinales), cualitativas o cuantitativas. En ocasiones la 
variable de respuesta Y, o variable dependiente, tiene más 
de dos categorías, para la que existen otras técnicas no 
tratadas aquí (1) y (3). Estos modelos son especialmente 
útiles cuando de manera natural o controlado 
experimentalmente, la probabilidad de ocurrencia del 
evento de interés (etiquetado con un 1) sigue una forma-S 
o sigmoidal, por ejemplo, en casos reportados de estudios 
clínicos, (1) y (3), el riesgo p de padecer hipertensión 
típicamente se comporta en la forma-S conforme aumenta 
la edad x(/.e., hasta cierta edad el riesgo es mínimo, luego 
conforme avanza la edad el riesgo crece paulatinamente 
hasta estabilizarse en un valor cercano a uno después de 
ciertos años); así mismo ocurre típicamente o de manera 
aproximada con el riesgo de no pagar un crédito conforme 
aumenta la deuda previa del solicitante (forma-S creciente) 



o al aumentar su ingreso (forma-S decreciente). Cuando el 
comportamiento no corresponde de manera aproximada 
a una forma-S, se recomienda utilizar otras técnicas 
distintas a la regresión logística, por ejemplo, árboles de 
clasificación o los modelos de riesgo aditivo (4). 

Para adentrarnos en el conocimiento de la técnica 
de regresión logística, utilizaremos como referencia los 
conceptos de la técnica de regresión lineal, ya que la 
regresión logística utiliza una metodología muy similar, 
excepto que no son necesarios supuestos de normalidad 
en la variable de respuesta ni en los residuos (lo cual es 
una de sus ventajas). Aunque ambas técnicas se basan en 
un ajuste lineal, la regresión logística utiliza primero una 
transformación de la variable de respuesta en unidades 
logarítmicas para deducir el modelo de regresión lineal 
final. El objetivo es estimar la probabilidad de ocurrencia 
p(x) de un evento, dado un conjunto de covariables 
predictoras x = (x v x 2 ,-~ ,x k ), ajustadas o controladas a 
ciertos valores particulares. El modelo matemático de 
forma-S comúnmente utilizado es la conocida función 
logística (para más detalles ver las citas (1 )-(4)): 

e f{x) e Po+PiX 1 +-+P k x k i 

PW = 1 + = 1~T e Po+PiXi+-+p k x k = 1 + e -p -ií =1 p iXi 

Una transformación sencilla en términos de lo que 
se conoce como razón de ventajas, expresada como el 
cociente de la probabilidad p(x) de ser moroso contra la 
probabilidad de ser cumplido 1- p(x), nos lleva al modelo 
de regresión en unidad logarítmica o función logit (del 
inglés log-unit): 

logit => In [ - ^ ] = f(x) = p + p lXl + - + P k x k 

El cual es ahora un modelo lineal para el que se deben 
estimar los parámetros jSs mediante la técnica de máxima 
verosimilitud, (1 )-(5), en vezde utilizar la técnica de mínimos 
cuadrados de la regresión lineal. Aunque, cabe mencionar 
que en el caso particular de la regresión lineal ambas 
técnicas producen las mismas estimaciones (1). El término 
entre corchetes en la ecuación anterior, se conoce como 
la razón de ventaja entre la probabilidad de ocurrencia de 
un éxito con respecto a la de fracaso, dado un valor de x 
particular. Esta razón de ventajas (odds de la literatura en 
inglés) es muy importante ya que nos ayuda a entender la 
relación entre clientes sujetos o no a un factor de riesgo, 
como puede ser tener o no dependientes económicos. Por 
ejemplo, en el caso simple de una variable dependiente 
predictora dicotómica,x= tiene dependientes económicos, 
denotando con x = 1 para "si" y x = para "no"; calculamos 
la razón de ventajas con ambos valores de x, al dividir 
ambos odds obtenemos lo que se conoce como el odds 
ratio OR (que representa la proporción en la que la razón 
de ventaja de estar expuesto a un factor supera a la de no 
estar expuesto). Por ejemplo, si el 0R - °^adX = =^) = 2 ' s '9 n 'f' ca 
que un cliente con cargas familiares (x = 1) tiene 2 veces 
más riesgo de no pagar el crédito que un cliente que no las 
tiene (x= 0). La interpretación también es válida en estudios 
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clínicos o epidemiológicos, donde, por ejemplo, se puede 
decir que el riesgo de contraer el VIH es dos veces mayor en 
personas drogadictas que las que no lo son. Lo interesante 
del modelo logístico estriba en que las estimaciones de 
los coeficientes del modelo de regresión (las jSs) sirven a 
su vez para estimar las odds ratios precisamente con la 
transformación logaritmo inversa, tal que ORi = e& (esta 
igualdad es fácil verificarla directamente de la definición 
anterior del OR). Así también para obtener los intervalos 
de confianza de las ORs, se calculan primero los intervalos 
de confianza de los parámetros (3s usando el estadístico 
Wald, el cual tiene una distribución aproximadamente 
normal bajo ciertas condiciones, no discutidas aquí ya que 
utilizaremos los resultaos directos del paquete SPSS, ver (1 ) 
y (3); después, mediante exponenciación de los límites de 
los intervalos obtenidos calculamos los intervalos para las 
ORs. 




MODELO DE RIESGO DE MOROSIDAD 
PARA CRÉDITOS BANCARIOS 

Dada la dificultad de vinculación con el sector empresarial 
para utilizar datos reales, presentamos un caso de estudio 
simulado sobre créditos bancarios, el cual intenta ampliar 
los resultados y análisis reportados para regresión 
logística en (2), ya que los autores comparan la capacidad 
predictora de un modelo empírico logístico frente a 
uno de análisis discriminante (http://www3.unileon.es/ 
pecvnia/pecvnia01/01_1 75_199.pdf). Nuestro objetivo es 
determinar que otros factores, además de los reportados en 
(2), caracterizan la morosidad en las entidades financieras, 
así como incluir variables adicionales de interés práctico 
como el ingreso, el cual sorprende no se contempla como 
variable significativa en (2). 

Nuestra propuesta incluye las cuatro variables 
reportadas con significancia estadística, a saber, la variable 
dependiente dicotómica Y, estatus de la morosidad (Y 
= 1 para moroso, Y = O, no moroso); y como variables 
independientes o predictoras: Destino del crédito (variable 
categórica: O = "traspaso de negocio"; 1 = "compra de 
automóvil"; 2 = "otros" el cual incluye compra de vivienda, 
remodelación, pago de deudas, compra de negocio o 



activos para actividad profesional y otros fines); Nuevo 
residente en lo localidad (1 = si, O = no); Número de impagos 
anteriores (O, 1, o 2) como variable discreta de razón. 
A diferencia de los autores de (2) y buscando tener un 
modelo más realista, incluimos tres variables adicionales 
consideradas de interés práctico, pero que fueron excluidas 
por los autores en (2) basados en la poca significancia 
estadística. Es bien sabido que la capacidad de pago está 
íntimamente relacionada con el Ingresos del cliente, el 
número de dependientes económicos (variable discreta con 
valores 0, 1 , 2, 3, 4, 5 o más) y su tipo de empleo, por lo que 
son agregadas al modelo. 

A continuación utilizaremos la técnica de regresión 
logística progresiva del paquete SPSS (5) para analizar 
tanto la significancia estadística mediante la prueba Wald 
como la prueba de de Hosmer-Lemeshow, para analizar 
qué tan bien se ajustan el modelo a los datos. Una vez 
que se obtiene la estimación de parámetros, se puede 
usar el modelo para medir su capacidad predictora, la 
interpretación de las ORs e intervalos de confianza. Para 
más detalles de justificación y uso de variables ver (2). 

Para seleccionar la muestra hipotética de clientes, es 
importante considerar los tipos de entidades financieras 
como bancos, cajas de ahorros, cooperativas de crédito. 
Se recomienda un muestreo por conglomerados en dos 
etapas donde las unidades primarias sean las entidades 
y las unidades últimas los clientes (nuestro estudio 
presupone tal diseño de muestreo). Además, utilizaremos 
1 00 casos en vez de los 72 reportados en (2); es importante 
notar que para la variable empleo, con nueve categorías 
el tamaño de muestra resulta muy pequeño para confiar 
estimar sus ORs, por lo que se debe tener cuidado con los 
resultados en casos reales, verificando que los resultados 
sean consistentes antes de asumirlos como válido, ya 
que para un usuario con poca experiencia podría pasar 
desapercibido (4). Asimismo, se debe utilizar un tamaño de 
muestra suficiente para garantizar que todas las categorías 
tienen una frecuencia de casos observados distinta de cero, 
para evitar divisiones con denominadores nulos al calcular 
los odds, en particular la categoría de referencia. Si este 
no fuera el caso, se recomienda fundir varias categorías 
adyacentes o similares en algún sentido práctico. Por 
cierto, esto es una de las posibles desventajas al aplicar 
regresión logística (4). 

GENERACIÓN DE LA BASE DE 
DATOS MEDIANTE SIMULACIÓN 

La base de datos utilizada (ver apéndice) fue simulada 
mediante el paquete SPSS. Cabe mencionar que dicha 
base sirve como conjunto de entrenamiento y prueba 
para generar el modelo, por lo que la variable morosidad 
es generada desde un principio para luego contrastarla 
con el modelo clasificador-pronosticador, una vez que se 
estiman las probabilidades o riesgo de morosidad. En este 
sentido, se da por sentado que la base de datos captura una 
supuesta relación intrínseca entre la variable de respuesta 
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y los valores de las covariables, los cuales usaremos para 
generar el modelo inmerso. 

Las variables fueron generadas de acuerdo a los 
siguientes criterios: 

y = Morosidad: Se genera una muestra de números 
aleatorios con distribución Bernoulli de parámetro p = 
0.26 que corresponde a la probabilidad de morosidad 
reportada en (2). 

x 1 = Destino: Números aleatorios con distribución 
binomial de parámetros n = 2 y p = 0.4. Donde = Traspaso 
de negocio, 1 = compra de automóvil 2 = otras inversiones 
(casa, local, activos, pago de deudas y otros gastos). Es 
posible utilizar un esquema discreto distinto a la binomial 
más realista, sin embargo, no hay pérdida de generalidad 
usando la binomial. 

x 2 = Nuevo Residente: Variable 0-1; muestra de 
números aleatorios con distribución Bernoulli con p = 0.35, 
suponiendo de manera empírica que el 35% son nuevos 
residentes. 

x 3 = Número Impagos: Variable aleatoria (v. a.) 
discreta con valores enteros 0, 1 y 2; más de tres impagos 
suponemos rechazada automáticamente la solicitud de 
crédito. Se genera mediante una binomial (n = 2, éxito = 
número de impagos en su dos operaciones anteriores, con 
probabilidad general de impago p = 0.08). Dicho valor del 
8% sólo es de referencia para generar la base, no confundir 
con la probabilidad de riesgo que buscamos predecir. 

x, = Salario: v. a. continua con distribución normal (en 
miles de pesos) con media 15 y desviación estándar 3. 
"Promedios" más o menos realistas de ingresos mensuales 
actuales para un trabajador de clase media-alta. 

x 5 = Empleo: v. a. con 9 categorías, generadas usando 
Excel con una distribución multinomial de acuerdo a la 
siguiente tabla de probabilidades propuestas. 



Tabla 1. Probabilidades por categorías para la variable 
Empleo 



Número de 
categoría 


Probabilidad 


Categoría de empleo 


1 


0.20 


Cuenta propia 


2 


0.20 


Empelado ejecutivo 


3 


0.24 


Empleado indefinido 


4 


0.14 


Empleado eventual 


5 


0.02 


Cuenta propia y 
empleado ejecutivo 


6 


0.01 


Cuenta propia y 
empleado indefinido 


7 


0.05 


cuenta propia y 
empleado eventual 


8 


0.04 


Jubilado o prejubilado 


9 


0.10 


Desempleado 



x 6 = Cargas familiares: variable discreta generada como 
una v.a. binomial con n = 5 p = 0.20, que representa la 
probabilidad de tener un hijo. Muchos casos resultaron 
con 0, 1 ó 2 hijos; lo que concuerda para clientes jóvenes 
solteros (29%, carga = 0); en su mayoría casados 47% 
(carga = 1), con un hijo (carga = 2) y no más de dos hijos 
(carga = 3). 

Se debe notar que cargas familiares se genera como 
variable discreta pero simbólicamente se emplea en el 
modelo como variable continua (en escala de razón) para 
analizar el aumento del riesgo por cada dependiente 
económico adicional (incremento unitario), en vez de 
comparar categorías por cantidad de hijos; o con cargas o 
sin cargas (dicotómica). Así también el número de impagos 
se emplea internamente como variable continua. Además, 
la variable destino del crédito es de varias categorías y sólo 
fueron significativas "traspaso de negocio" o "compra de 
automóvil". Es importante mencionar que los autores en 

(2) tomaron las categorías como variables independientes, 
lo cual no resulta congruente al estudiar el modelo. Las 
categorías de remodelación o compra de vivienda, compra 
de local o activos para actividad profesional, pagar deudas 
y otros gastos no fueron variables significativas y se 
agrupan en la categoría de "otros". 

DISCUSIÓN DE RESULTADOS 

La base de datos del apéndice fue procesada con SPSS 
(5); las variables con más de dos categorías (empleo y 
destino) fueron codificadas con "celda de referencia" 

(3) , usando 1 para dicha categoría y cero para el resto; la 
categoría de referencia fue la última y surge cuando todas 
las demás categorías toman el valor cero. Lo primero es 
una evaluación general del modelo mediante la prueba 
-2Log de la Verosimilitud, la cual se utiliza para decidir si 
es el modelo es confiable estadísticamente hablando 
para predecir la variable de respuesta. El SPSS utiliza 
como hipótesis nula que el modelo es significativo y, 
como hipótesis alternativa lo contrario, asumiendo que el 
estadístico -2Log tiene una distribución chi cuadrada con 
k -2 grados de libertad. El estadístico -2Log resultó 82,880 
con 14 grados de libertad y un nivel de significancia de 
0.246, por lo tanto, no rechazamos la hipótesis nula de 
que el modelo logístico en general es significativo como 
modelo predictor. Analicemos ahora el ajuste de los datos. 

La tabla 2 contiene la bondad de ajuste, la cual resulta 
muy significativa 55.9%. Lo anterior se aprecia en la tabla 3 
que muestra la diferencia entre las frecuencias observadas 
y esperadas muy aproximadas, agrupando los valores 
de probabilidad estimada en 10 clases percentiles. Para 
descripción de la prueba de bondad ver (1) y (4). 



Tabla 2. Prueba de bondad de ajuste Hosmer y Lemeshow. 





gi 


Sig. 


6,792 


8 


0.559 
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La tabla 4 muestra la clasificación de los casos observados 
contra los casos pronosticados, dichos resultados se pueden 
obtener comparando en la tabla del apéndice la primera 
columna (clase observada), contra la última columna que 
representa la probabilidad p de ser clasificado como moroso 
si dicho valor es mayor que el punto de corte utilizado, Pe = 
0.15. El modelo tuvo una tasa total de pronósticos correctos 
del 70%, es decir, los pronósticos tanto de morosos como no 
morosos clasificados de manera correcta con respecto a su 
clase observada. 

Tabla 3. Prueba de Hosmer y Lemeshow para frecuencias 
observadas y esperadas. 





Si el crédito fue 
pagado o no = No 
moroso 


Si el crédito fue 
pagado o no = Moroso 


Total 


Observado 


Esperado 


Observado 


Esperado 




Paso 


1 


10 


9.819 





0.181 


10 




2 


10 


9.556 





0.444 


10 




3 


9 


9.319 


1 


0.681 


10 




4 


10 


9.146 





0.854 


10 




5 


10 


8.852 





1.148 


10 




6 


7 


8.226 


3 


1.774 


10 




7 


7 


7.521 


3 


2.479 


10 




8 


5 


6.909 


5 


3.091 


10 




9 


6 


6.117 


4 


3.883 


10 




10 


6 


4.536 


4 


5.464 


10 



Tabla 4. Clasificación-predicción del modelo. 



Observado 


Pronosticado 


Porcentaje 




No moroso Moroso 


correcto 


No moroso 


51 29 


63.8 


moroso 


1 19 


95.0 


Porcentaje global 




70.0 



La tasa de clasificación de verdaderos negativos (no 
morosos pronosticados como no morosos) fue del 63%; 
mientras que la tasa más relevante en este caso, clasificar 
correctamente a los clientes morosos (verdaderos 
positivos) fue de 95%, demostrando gran capacidad 
predictiva del modelo. Cabe mencionar que el Pe se eligió 
con base en la experiencia de la simulación, sin embargo, 
formalmente se puede seleccionar usando la llamada curva 
de características de operación o curva ROC (Receiver 
Operating Characteristic), consultar (1) y (3). 
Tabla 5. Estadísticas para el modelo de regresión logística. 
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B 


E.T. 


Wald 


g" 


Sig. 


Exp(B) 


I.C. 95% para EXP(B) 


Inferior 


Superior 


Destino 






1.602 


2 


0.449 








Destino (1) 


0.897 


0.782 


1.316 


1 


0.251 


0.408 


0.088 


1.888 


Destino (2) 


0.715 


0.687 


1.086 


1 


0.297 


0.489 


0.127 


1.878 


Nuevo residente 


1.096 


0.687 


2.543 


1 


0.111 


0.334 


0.087 


1.285 


No. de impagos 


1.307 


0.568 


5.294 


1 


0.021 


3.696 


1.214 


11.255 


Salario 


0.217 


0.098 


4.872 


1 


0.027 


0.805 


0.664 


0.976 


Cargas 


0.207 


0.334 


0.385 


1 


0.535 


1.231 


0.639 


2.370 


Empleo 






1.882 


8 


0.984 








Empleo (1) 


0.483 


1.052 


0.211 


1 


0.646 


0.617 


0.078 


4.850 


Empleo (2) 


0.507 


1.084 


0.219 


1 


0.640 


0.602 


0.072 


5.044 


Empleo (3) 


0.540 


1.104 


0.240 


1 


0.625 


0.583 


0.067 


5.069 


Empleo (4) 


0.318 


1.080 


0.087 




0.769 


1.374 


0.166 


1 1 .400 


Empleo (5) 


19.286 


25577.4 


0.000 




0.999 


0.000 


0.000 




Empleo (6) 


0.939 


1.567 


0.359 




0.549 


2.557 


0.119 


55.122 


Empleo (7) 


18.062 


40192.9 


0.000 




1.000 


0.000 


0.000 




Empleo (8) 


0.527 


1.595 


0.109 




0.741 


0.590 


0.026 


13.452 


Constante 


2.402 


1.719 


1.952 




0.162 


1 1 .043 







Nota. El número entre paréntesis en el nombre de cada variable indica su categoría. 



La tabla 5 contiene en la columna B la estimación de 
en la columna Wald, el cual es un indicador para probar 
la significancia (Sig.) estadística de cada variable (1), 
en este caso se utilizó un nivel de significancia del 
5%. Se observa en la sexta columna que las variables 
significativas fueron número de impagos previos y el 
salario; mientras que la variable nuevo residente es 
significativa al 12%. La variable empleo no resulta 
significativa en ninguna de sus nueve categorías. Cabe 
mencionar que se hicieron pruebas eliminando estas 
variables del modelo, sin embargo, su desempeño 
fue más pobre, por lo que se mantuvieron. La tabla 
5 también incluye las estimaciones de interés de las 
ORs y cuya interpretación se da más abajo; además, 
se incluyen sus intervalos de confianza con propósitos 
ilustrativos; para una discusión más amplia ver (1 ), (3) y 
(4). Observe que los coeficientes ps negativos implican 



los coeficientes ps del modelo y el estadístico de prueba 
un OR menor a 1 como indicador de disminución del 
riesgo a medida que aumenta la variable predictora; un 
coeficiente p positivo mucho mayor que cero no influirá 
en el cálculo de la probabilidad, mientras que un valor 
negativo mucho menor que cero producirá por si sola 
una probabilidad estimada aproximadamente cero. Lo 
anterior debido al cambio de signo en la función logit. 

Si usamos los coeficientes estimados de la tabla 
5, columna B, y los valores de las variables para 
cada caso en la base de datos simulada, podemos 
estimar la probabilidad o riesgo de morosidad de cada 
cliente mediante el modelo de regresión logística. Por 
ejemplo, los cálculos para el caso 1 muestran el valor 
de predicción de riesgo de morosidad siguiente: 



Caso 


Morosidad 


Destino 


Nuevo 
residente 


Número de 
impagos 


Salario 


Empleo 


Cargas 


Predicción 
















18.047 


2 


1 


0.0623 


x 2 


1 





1 





1 2.040 


1 


1 


0.0773 



Predicción de riesgo de morosidad para el caso 1 : 



1 + exp [-(2.402 - 0.897 x Destino (1) - 0.217 x Salario - 0.507 x Empleo(Z) + 0.207 x Cargas)] 

1 



P(X t ) = 



1 + exp [-(2.402 - 0.897 X 1 - 0.217 X (18.047) - 0.507 x 1 + 0.207 X 1)] 



= 0.0623 
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Observe que el caso 1 representa datos de un cliente 
no moroso, cuyo destino del crédito fue para compra 
de automóvil, no es nuevo residente, no tiene impagos 
previos, tiene un salario elevado y con una carga familiar, 
el cual podríamos decir que es un perfil típico de un buen 
cliente. Observe que la probabilidad estimada por el 
modelo resultó menor que el Pe = 0.15, por lo que este 
primer caso es clasificado de manera correcta como cliente 
no moroso. Es necesario aclarar que en la variable destino 
la codificación de categoría de referencia es la última, 
en este caso la tercera categoría "otros"; así, la categoría 
representa a "traspaso de negocio" mientras que la 
categoría 1 representa "compra de automóvil", esto explica 
por qué se usó el se usó el valor uno para indicar la variable 
destino(l) en el cálculo anterior. Las probabilidades 
restantes se calculan de manera similar y se muestran en 
la columna 8 del apéndice. Debido a que la probabilidad 
estimada para los primeros dos casos es menor que el Pe = 
0.15, estos serán pronosticados no morosos, sin embargo, 
según los datos del apéndice, observe que el segundo caso 
corresponde en realidad a un cliente moroso, por lo que la 
clasificación del modelo para este caso es incorrecta, de 
hecho es el único "caso extremo-moroso" presente en los 
datos simulados. 

El valor 3.696 en la columna Exp (B) nos dice que por 
cada crédito impagado (incremento unitario) el riesgo de 
no devolver el crédito aumenta 3.696 veces. Asimismo, 
vemos que por cada $ 1,000 de incremento (unitario) 
en el salario, el riesgo (multiplicativo) de no pagar el 
crédito disminuye un 80.5%; es decir, si un cliente con 
un salario determinado presenta un riesgo del 50%, si 
mantenemos fijas todas las variables pero incrementamos 
el salario en $1,000, el riesgo (multiplicativo) se reduce 
aproximadamente al 40%; aún más, si el salario aumenta 
en incrementos de $ 5,000, el riesgo se reduce hasta un 
e 5x(-o.2i7) x 10 o « 34o/ . También vemos una correlación 
negativa en el renglón de nuevo residente, por lo tanto, el 
(Tk = 0.334 significa que los nuevos residentes (categoría 1) 
tienen una disminución del riesgo por un factor de 33.4% 
comparado con un residente antiguo. 



CONCLUSIONES 

Se mostró que un modelo de regresión logística para 
pronosticar la morosidad de un crédito bancario se 
desempeña adecuadamente como modelo clasificador 
sobre una base de datos simulada; en consecuencia, se 
asume factible que este tipo de modelos sea adecuado 
sobre casos reales o similares. El estadístico - 2Log de 
Verosimilitud mostró que en general el modelo es un 
buen predictor; además, se mostró que se ajusta bien a 
los datos, tanto en las frecuencias observadas contra las 
esperadas. Se obtuvo un excelente desempeño del 95% 
para predecir particularmente posibles clientes morosos. 
La tasa de predicción global del modelo fue buena de 
acuerdo a la literatura estándar, ya que se ajusta a un 70% 
de clasificaciones correctas. Las variables más significativas 
para pronosticar el riesgo de morosidad son número de 
impagos previos, salario y nuevo residente. Sorprende que 
la variable número de dependientes económicos o carga 
familiar resultó no significativa. Es importante recalcar 
que el modelo fue contrastado con los datos simulados, 
la metodología aquí presentada se puede ver como 
referencia para un trabajo futuro más realista, en tal caso, 
se debe tener cuidado en validar primero la base de datos 
mediante el estudio de las distribuciones adecuadas, así 
como validación final del modelo, o mejor aún, usar una 
base de datos reales. 
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APÉNDICE: Base de datos simulados para el modelo de 
riesgo para créditos bancarios 

La siguiente tabla contiene la base de datos para 100 
clientes obtenida mediante simulación para generación 
del modelo; dichos datos pueden ser recuperados del 
repositorio ftp.mat.uson.mx/~claudio o replicando el 



CASO 


D 


■1 


X 2 


X 3 


X 4 


X 5 


X 6 


p 


1 


n 

u 


u 








18.047 


2 


1 


0.062 


2 


1 

I 


n 

u 


1 





1 2.040 


1 


1 


0.077 


3 


1 

I 


9 
z 








16.237 


3 





0.159 


4 


1 

I 




1 


2 


1 1 .670 


8 





0.491 


5 


n 


9 
z 








1 1 .674 


8 


1 


0.389 


6 


n 

u 


n 

u 








17.786 


9 





0.087 


7 


o 


1 

1 


1 





12.591 


1 


1 


0.082 


8 





1 


1 





19.416 


3 


2 


0.023 


9 





1 


1 





18.149 


8 


1 


0.025 


10 





1 








18.223 


1 


1 


0.073 


11 





1 





1 


18.504 


2 


1 


0.210 


12 


1 











13.499 


1 


3 


0.217 


13 














1 7.870 


1 


2 


0.080 


14 





1 


1 





15.187 


1 


2 


0.059 


15 


1 


2 





1 


15.064 


2 


1 


0.535 


16 








1 


1 


16.303 


1 


2 


0.131 


17 





1 


1 





1 1 .882 


2 


2 


0.111 


18 





1 








14.164 


2 





0.131 


19 





1 








19.224 


3 


1 


0.056 


20 





2 


1 





14.898 


2 


1 


0.097 


21 





2 


1 


2 


16.009 


2 


3 


0.637 


22 





1 





1 


14.273 


2 


1 


0.401 


23 





2 


1 





17.359 


9 


1 


0.095 


24 





2 


1 





14.648 


3 


2 


0.120 



experimento descrito en la sección de simulación, por 
lo que pueden omitirse (se incluyen para propósitos 
de revisión). El significado de cada variable x 7 ,..., x 6 , en 
esta tabla es como se definió también en la sección de 
simulación. 
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INTRODUCCIÓN 

El intercambio de los bienes inmuebles se desarrolla 
principalmente en un mercado secundario notablemente 
imperfecto, donde la opacidad de los precios de operación 
son la constante y sin embargo sus participantes deben 
decidir los montos en que se intercambian los bienes 
¡nmuebles;asítambién,el mercado primarioes desarrollado 
por un reducido grupo de oferentes, dificultando al 
comprador la elección a partir de una limitada gama de 
características del inmueble. 

La ausencia de criterios objetivos que ayuden en la 
identificación de los elementos que definen estos valores 
es motivo para que los precios iniciales de oferta confundan 
a los participantes y generen ineficiencias en la realización 
de sus operaciones de compra-venta. En este sentido es 
importante estudiar el comportamiento del mercado en 
relación con los atributos principales de los inmuebles; 
esto visto desde la perspectiva del comprador que define 
su disposición a pagar por estas características del bien (1 ). 

La adquisición de inmuebles habitacionales puede 
deberse a motivaciones diferentes en el individuo: la 
principal, solucionar su necesidad de abrigo ante la 
naturaleza, pero también, de los inconvenientes que 
genera el convivir en concentraciones humanas, como es la 
necesidad creada de intimidad; así también, la adquisición 
de la segunda residencia como una oportunidad para 
cambiar su entorno por otro normalmente diferente 
que le ayude a disfrutar mejor sus momentos de ocio; y 
finalmente, pero no menos importante, la adquisición de 
inmuebles habitacionales con el propósito de obtener 
ingresos futuros que le permitan conservar el valor original 
del capital invertido y obtener un monto adicional (renta). 

La aplicación de la técnica de modelos de precios 
hedónicos presenta la característica de explicar el 
comportamiento del precio de los bienes inmuebles a 
partir de los precios de los diferentes atributos que lo 
conforman. Esto significa que las viviendas se constituyen 
de una serie de elementos que las hace atractivas para sus 
compradores y los precios de estos elementos pueden 
ser estimados a partir de información del mercado. Así 
pues, conociendo dichos elementos es posible explicar y 
construir el precio de cada una de las viviendas del lugar. 

Lo anterior permite mostrar cómo los precios de las 
viviendas en la ciudad son resultado de un conjunto de 
características intrínsecas y extrínsecas propias de cada 
inmueble, las cuales es posible establecer y modelar 
estadísticamente en el espacio urbano (2). Para esto se 
desarrollan los modelos que explican el comportamiento 
de los precios de la vivienda y se muestran los principales 
elementos que proporcionan valor a las viviendas en cada 
sector. 

CONCEPTOS FUNDAMENTALES 

La formación espacial de los valores urbanos pude definirse 
con base en la teoría económica espacial; esto, a partir 
de la formación espacial de los valores urbanos del 



mercado residencial, en donde es posible identificar 
tres componentes principales: la accesibilidad, las 
externalidades físicas y ambientales y los factores sociales. 

La accesibilidad se define como la característica que 
presenta el inmueble dentro del territorio urbano, la cual 
permite al individuo tener desplazamientos más cortos o 
de menor tiempo a los lugares de interés, normalmente 
al centro de trabajo, de educación o equipamientos 
generales de la ciudad (3). La accesibilidad puede medirse 
a partir de la distancia euclidiana, la distancia de recorrido 
o el tiempo de traslado entre dos puntos de la ciudad y 
es normalmente asociada como un factor importante de 
influencia en el valor del inmueble (4). 

El crecimiento de las ciudades genera externalidades 
físicas y ambientales, positivas como son las economías de 
aglomeración y negativas como la presencia de polución, 
congestión vehicular, disminución del espacio verde, entre 
otras; elementos que presentan efectos importantes en el 
comportamiento de los valores de los inmuebles urbanos 
(5). 




Por esto, los entornos urbanos bien planeados y 
ejecutados, que consideran aspectos tecnológicos, del 
medio ambiente y edificaciones de calidad, que armonizan 
las características de la zona, son productos inmobiliarios 
cada vez más preferidos por el consumidor. 

La ubicación de la vivienda en una zona ocupada 
por individuos con alto nivel de renta, niveles educativos 
elevados o excelente nivel cultural, favorece de forma 
importante el incremento del valor de las viviendas (6). 
Es decir, el atractivo de una zona suele estar ligado a 
diferentes fenómenos económicos y sociales que generan 
atracción y disposición del comprador a pagar más por la 
vivienda; lo anterior, es producto principal de la imagen de 
jerarquía social asociada al lugar. 
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METODOLOGÍA 

Objeto de estudio 

El área de estudio debe ser un referente importante que 
pueda ser ubicado fácilmente, que facilite el acceso a la 
información necesaria para desarrollar la metodología, 
que presente características similares a las supuestas 
en la hipótesis y, sobre todo, que sea de interés para el 
investigador (Figura 1). 




Figura 1 . Ciudad de Hermosillo, Sonora, México. 
Fuente: Google Earth. 



El sitio y objeto de estudio seleccionado es el mercado 
de venta de vivienda de la ciudad de Hermosillo, Sonora, 
México; se localiza 29°05'56" de latitud N y 1 10°57'15" de 
longitud O con una altitud de 216 msnm (Figura 2) en 
una mancha urbana sensiblemente plana de 160.93 km 2 , 
donde 125.68 km 2 corresponden a superficies de lotes y 
predios urbanos, de los cuales el 37.25% se encuentran 
baldíos (7). 




Figura 2. Distribución de lotes baldíos y edificados. 

Fuente: Programa de Desarrollo Urbano del Centro de Población 
de Hermosillo (2006). 



El municipio de Hermosillo registra una población de 
784,342 habitantes y un mercado de vivienda compuesto 
por 213,304 viviendas particulares habitadas, donde 
el 95.24% posee agua entubada, el 94.59% sistema de 
drenaje y el 97.49% servicio de energía eléctrica (8). 

Bases de datos 

Se recolecta información de operaciones de vivienda 
directamenteen campeen prensa,enportaleselectrónicos, 
con profesionales de la valuación y unidades de valuación, 
integrando una base de 450 inmuebles distribuidos en la 
ciudad. En éstos se describen características como: 

❖ Precio de oferta. 

❖ Superficie de construcción. 

❖ Superficie de terreno. 

❖ Número de recámaras. 

❖ Número de baños. 

❖ Número de plantas. 

❖ Calidad de los servicios e instalaciones. 

❖ Cercanía a vialidad principal. 

❖ Cercanía a centro comercial. 

❖ Cercanía al centro de población. 
♦:♦ Número de huellas de cochera. 

Además se agrega información de la Dirección de 
Catastro Municipal Hermosillo correspondiente a los 
valores catastrales de los terrenos urbanos de la ciudad, 
a partir de la cual se integran los atributos de valores 
de terrenos por zona homogénea y por banda de valor 
vigentes al 2012. 

Sectorización de la ciudad 

Parece prudente utilizar la subdivisión oficial de la ciudad 
de Hermosillo, correspondiente a 450 colonias, como las 
subdivisiones del mercado de la vivienda; sin embargo la 
pulverización del territorio con esta sectorización dificulta 
enormemente su utilización en el estudio. Así también, 
se explora la denominación oficial de las 97 regiones 
catastrales del municipio de Hermosillo; aun así, después 
de una simple inspección se observa que no representan 
zonas homogéneas en cuanto a características de calidad 
de las viviendas. 

Buscando una subdivisión del mercado de la vivienda 
que permita obtener modelos con grados de explicación 
aceptables, se procede a utilizar información oficial de la 
Dirección de Catastro Municipal de Hermosillo (9) relativa 
a los valores catastrales del suelo por zona homogénea, 
se procede a sectorizan en regiones de valores de suelo 
definidas a partir de quintiles, es decir, realizando una 
distribución equitativa de colonias por cinco rangos de 
valores catastrales del suelo. 

Con base en esta clasificación, se presenta la ciudad 
de Hermosillo dividida en sectores de rangos de valores 
catastrales de suelo, las cuales se muestran en la figura 3. 
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Figura 3. Sectorización de la ciudad de Hermosillo. 

En la figura 3 es posible observar comportamientos 
consistentes con la teoría de formación espacial de 
submercados de los valores inmobiliarios, es decir 
diferencias en valores para una misma característica 
de un sector a otro (10). Así también, las figuras 4, 5, 6, 
7 y 8 muestran correlaciones lineales aceptables entre 
la superficie construida y los precios de venta en los 
diferentes sectores. 
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Figura 4. Relación entre la superficie construida y el precio 
de venta. 
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Figura 5. Superficie construida contra precio de venta. 
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Figura 7. Superficie construida y precio de venta en rango 
de $691 a$1165/m 2 . 
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Figura 8. Relación entre la superficie construida y el precio 
de venta en rango de $ 1 1 66 a $3200/m 2 . 

Los índices de correlación mayores parecen asociarse 
a inmuebles ubicados en sectores de valores catastrales 
altos, lo que deja entrever que los valores catastrales 
bajos corresponden a edificaciones más heterogéneas 
propias de los procesos de autoconstrucción donde las 
expectativas de valores de venta suelen caracterizarse por 
una alta variabilidad. 
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MODELOS EXPLICATIVOS 
DEL VALOR DE LA VIVIENDA 

Los modelos de precios hedónicos inician su gestación con 
el trabajo de Hass en 1922, quien intenta obtener el valor 
de tierras dedicadas la actividad agrícola; posteriormente 
Court, en 1939, realiza una aplicación para obtener 
precios en el mercado automovilístico; Rosen en 1974 
realiza una aplicación en el mercado inmobiliario, donde 
conceptualiza el bien inmueble a partir de un conjunto 
de atributos que satisfacen las diferentes necesidades o 
gustos de la demanda, es decir, un bien inmueble urbano 
que está compuesto por una serie de características 
deseables por el comprador, las cuales pueden integrarse 
en la siguiente ecuación: 

Y = P Q + + p 2 X 2 + ¡3 3 X 3 . . . + p n X n + e 
Donde: 

Y variable dependiente (valor del inmueble). 

X variable independiente (característica del inmueble). 

P coeficiente de la variable producto del modelo 

de regresión. 
e error de estimación del modelo. 

La solución de la expresión del modelo propuesto se 
lleva a cabo con la aplicación de la técnica de regresión lineal 
múltiple aproximada por mínimos cuadrados ordinarios. 
Este proceso se lleva a cabo con el programa estadístico 
SPSS 18, donde se aplica la técnica de pasos sucesivos 
(stepwise). Además se revisan las condiciones que debe 
cumplir un buen modelo, como es el comportamiento 
normal de la variable dependiente, valores adecuados 
de multicolinealidad, relaciones lineales y coeficientes 
significativos, comportamiento homocedástico, entre 
otros. 

RESULTADOS 

MPH del Sector I de $1 1 4.00 a $ 435.00 

P.V .=-23,063.68+1 1 8,61 8.63 C. Vivienda+968.1 6 S.Construida 



] 



El modelo resultante con una explicación del 74% del 
precio de venta, se fundamenta en dos variables la Calidad 
de la vivienda y la superficie construida. La Calidad de la 
vivienda con un coeficiente (3 tipificado de 0.630, mientras 
que la Superficie construida con (3 tipificado de 0.394; esto 
es la Calidad de la vivienda es la principal característica 
que define el valor de la vivienda seguida de la Superficie 
construida. 

El modelo muestra que el incremento de 1 
nivel de Calidad de la vivienda, incrementa el precio 
aproximadamente en $1 18,000.00, mientras que cada m 2 
de Superficie construida tiene un impacto de $968.1 6 en el 
precio. Por otro lado, la pregunta obligada ¿Qué sucede 
con el resto de las variables independientes? son variables 
que no se incluyen en el modelo debido a que presentan 
una significancia pobre en el comportamiento del valor de 
la vivienda, o bien, presentan colinealidad con las variables 
independientes seleccionadas. 

MPH del Sector II de $436.00 a $540.00 

P.V .=-54,472.1 98+1 20,944.05C.V¡v¡enda+33,51 2.08A.Centro Comercial 

El modelo resultante con una explicación del 54% del 
valor de la vivienda se fundamenta en dos variables, la 
Calidad de la vivienda y la Accesibilidad a centro comercial. 
La Calidad de la vivienda con un coeficiente (3 tipificado 
de 0.722, mientras que la Accesibilidad a Centro comercial 
con (3 tipificado de 0.281; esto es la Calidad de la vivienda 
es la característica principal en la definición del valor de 
la vivienda y en menor grado la Accesibilidad a Centros 
comerciales. 

El modelo muestra que el incremento de 1 nivel 
de Calidad de la vivienda, incrementa el precio en 
aproximadamente $120,944.00, mientras que cada nivel 
de Accesibilidad a Centro comercial tiene un impacto de 
$33,512.06 en el precio. El resto de las variables tienen un 
comportamiento similar al del modelo anterior. 

MPH del Sector III de $541 .00 a $ 690.00 

P.V.=-51 ,581 .44+1 32,570.41 C. Vivienda+1 ,545.56 S.Construida 
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El modelo resultante con una explicación del 85% del 
precio de venta, se fundamenta en dos variables, la Calidad 
de ¡a Vivienda y la Superficie construida. La Calidad de la 
vivienda con un coeficiente (3 tipificado de 0.690, mientras 
que la Superficie construida con (3 tipificado de 0.332; esto 
es la Calidad de la vivienda es la principal característica 
que define el valor de la vivienda seguida de la Superficie 
construida. 

El modelo muestra que el incremento de 1 nivel 
de Calidad de la vivienda, incrementa el precio en 
aproximadamente $132,570.41, mientras que cada m 2 de 
Superficie construida tiene un impacto de $1,545.56 en el 
precio. El resto de las variables tienen un comportamiento 
similar al de los modelos anteriores. 

MPH del Sector IV de $ 691.00 a $ 1,165.00 

P.V .=-1 93,782.99+1 1 8,028.44 C.Vivienda+4,224.34 S.Construida 

El modelo resultante con una explicación del 73% 
del precio de venta, se fundamenta en dos variables, la 
Calidad de la vivienda y la Superficie construida. La Superficie 
construida con un coeficiente (3 tipificado de 0.654, 
mientras que la Calidad de la vivienda con (3 tipificado de 
0.297; esto es, en este sector la Superficie construida es la 
principal característica que define el valor de la vivienda 
seguida de la Calidad de la vivienda. 

El modelo muestra que el incremento de 1 nivel 
de Calidad de la vivienda, incrementa el precio en 
aproximadamente $118,000.00, mientras que cada m 2 de 
Superficie construida tiene un impacto de $4,224.00 en el 
precio. El resto de las variables tienen un comportamiento 
similar al de los modelos anteriores. 

MPH del Sector V de $ 1,1 66.00 a $ 3,200.00 

P.V =-236,1 75.94+1 51,666.05 C. Servicions e lnstaL+8,606.98 S. Construida 

El modelo resultante con una explicación del 77% 
del precio de venta, se fundamenta en dos variables, la 
Calidad de servicios e instalaciones y la Superficie construida. 
La Superficie construida con un coeficiente (3 tipificado de 
0.731, mientras que la Calidad de servicios e instalaciones 
con (3 tipificado de 0.249; esto es la Superficie construida es 
la principal característica que define el valor de la vivienda 
y en segundo lugar la Calidad de servicios e instalaciones. 

El modelo muestra que el incremento de 1 nivel de 
Calidad de servicios e instalaciones, incrementa el precio en 
aproximadamente $151,666.00, mientras que cada m 2 de 
Superficie construida tiene un impacto de $8,607.00 en el 
precio. El resto de las variables tienen un comportamiento 
similar al de los modelos anteriores. 

CONCLUSIONES 

Los resultados muestran que es posible explicar el precio 
de la vivienda en la ciudad de Hermosillo a partir de la 
aplicación de modelos de precios hedónicos, si bien el 
grado de explicación de los cinco modelos presentan 



diferentes aproximaciones, es posible mejorar su 
explicación y uniformidad utilizando una base informativa 
más amplia del mercado inmobiliario de la vivienda en la 
ciudad. 

El comportamiento de los valores de la vivienda en la 
ciudad puede ser explicado fundamentalmente con base 
en la Calidad de la vivienda, sobre todo en los sectores I, II 
y III, mientras que la Superficie construida es el elemento 
principal que caracteriza el precio en los sectores IV y 
V. Esto es, la calidad de los materiales utilizados en la 
edificación y el proyecto arquitectónico, son los elementos 
predominantes en los sectores I, II y II, mientras que el 
tamaño de las viviendas es fundamental en los sectores IV 
yV. 

Lo anterior parece indicar que el valor de los sectores 
bajos (I, II y III) queda definido en primera instancia por 
la calidad de los materiales y proyecto, mientras que los 
niveles altos, consideran esto invariante y las diferencias 
de valor están fundamentadas por la superficie construida. 

Por otro lado, es importante destacar que siendo 
la ciudad de Hermosillo una urbanización extendida, la 
accesibilidad al centro de la ciudad, vialidad principal 
y centro comercial, no muestran una importancia 
significativa. Lo anterior puede deberse a la existencia de 
diferentes puntos importantes en la ciudad, premisa bajo 
la cual parece desarrollarse y crecer, restando importancia 
al centro histórico y de oficinas públicas. Otra razón puede 
ser, la notable presencia de transporte privado en las 
familias de la ciudad, que prioriza el uso del automóvil 
sobre la accesibilidad de la vivienda. 
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INTRODUCCIÓN 

Este año se celebra un siglo de la teoría atómica de Bohr, la 
cual representa la transición entre la mecánica clásica y la 
mecánica ondulatoria o cuántica. Finalmente la mecánica 
cuántica ganó la batalla y en la actualidad es la que se usa 
para tratar fenómenos a nivel atómico. 

En 1913 el físico Niels Bohr (1) planteó los postulados 
que hoy llevan su nombre y que le permitieron proponer 
un esquema semiclásico, para explicar los espectros de 
emisión y absorción del átomo de Hidrógeno y sentar 
las bases para comprender estos fenómenos en átomos 
complejos. Esta nueva teoría se basa, por un lado, en 
el trabajo de Planck (2), en el cual se establecía que el 
intercambio de energía entre la radiación y la materia se 
realiza de forma discreta, fenómeno que se conoce como 
la cuantización de la energía, así como en el de Rutherford 
(3), que estableció la estructura moderna del átomo con 
el experimento de difracción de partículas alfa por átomos 
de oro. 

En este trabajo se aplica la teoría atomística de Bohr 
al problema de la interacción de dos electrones con un 
núcleo en la aproximación electrostática. Se resuelve 
numéricamente el problema de los tres cuerpos (4) usando 
la mecánica clásica y los postulados de Bohr, con el fin de 
buscar soluciones estables para garantizar la confiabilidad 
de los resultados encontrados. 

Primero se establece un algoritmo numérico, basado 
en el movimiento con aceleración uniforme para tiempos 
ultracortos, para encontrar la solución al problema 
planteado y que satisfagan los postulados de Bohr, 
posteriormente se aplica dicho formalismo al problema de 
dos electrones interaccionando con un núcleo. 

Enseguida se transcriben los postulados de Bohr (1 ) en 
su versión original en inglés: 

That the dynamical equilibrium of the systems ¡n the 
stationary states can be discussed by help of the ordinary 
mechanics, while the passing of the systems between 
different stationary states cannot be treated on that basis. 

Que el equilibrio dinámico de los sistemas en estados 
estacionarios puede ser discutido con la ayuda de la 
mecánica ordinaria, mientras que el paso de los sistemas 
entre diferentes estados estacionarios no puede ser 
tratado mediante esta base. 

That the latter process is followed by the emission of 
a homogeneous radiation, for which the relation between 
the frequency and the amount of energy emitted is the 
one given by Planck's theory. 

Que el anterior proceso es seguido por la emisión de 
una radiación homogénea, para la cual la relación entre la 
frecuencia y la cantidad de energía emitida es aquélla dada 
por la teoría de Planck. 

The different stationary states correspond to the 
emission of a different number of Planck's energy-quanta 
and, 

Los diferentes estados estacionarios corresponden a 
la emisión de diferentes números cuánticos de energía de 



Planck y, 

That the frequency of the radiation emitted during the 
passing of the system from a state in which no energy is 
yet radiated out to one of the stationary states, is equal to 
half the frequency of revolution of the electrón in the latter 
state. 

Que la frecuencia de la radiación emitida durante 
el paso del sistema de un estado en el cual no hay aún 
energía radiada a uno de los estados estacionarios, es igual 
a la mitad de la frecuencia de revolución del electrón en el 
último estado. 

EL ÁTOMO DE HIDRÓGENO Y 
EL ALGORITMO NUMÉRICO 

Un electrón orbitando alrededor de un protón constituye 
un átomo de hidrógeno, el elemento más ligero y 
abundante del universo. Estas partículas están sujeta a 
fuerzas de tipo electromagnético y gravitacional, esta 
última se desprecia ya que su magnitud es pequeña 
comparada con la magnitud de la fuerza electromagnética. 

Usando la segunda ley de Newton y considerando 
únicamente la fuerza de Lorentz, la ecuación del 
movimiento del protón es: 

mpOo(t) = +eÉ 1 [? (t),t] + ev (t) x B^Ct), t] (1) 

Donde a (t), v (t) , m p y +e son la aceleración, la 
velocidad, lajnasa y la carga del protón, respectivamente; 
Ei(r Ql t) y fíi(r ,t) representan el campo eléctrico y 
magnético producido por el electrón en la posición del 
protón r (t). Mientras que, ecuación del movimiento del 
electrón está dada por: 

m e a x {t) = -eÉofaífXt] - ev ± (t) x B [*i(O,t] O) 

Donde ¿^(t), v t (t) , m e y -e son la aceleración, la 
velocidad, lajnasa y la carga del electrón, respectivamente; 
É (r lt t) y B (f lf t) representan el campo eléctrico y 
magnético producido por el núcleo en la posición del 
electrón r t (t). 

Los campos electromagnéticos producidos por cargas 
en movimiento son matemáticamente complicados y más 
adelante se harán algunas simplificaciones para resolver 
(Dy(2). 

Para encontrar una solución numérica a este problema 
se debe establecer las condiciones iniciales. Las posiciones 
y velocidades iniciales del núcleo y del electrón son: 

? (0) = -fdj, 7^(0) = +d& v (0) = -fvj vi(0) = +v¡J (3) 

Donde / = m e /m p « 1/1836. Con esta selección de 
posiciones y velocidades, el centro de masa, que es el 
punto donde se supone concentrada toda la masa del 
sistema, y la velocidad del centro de masa son cero. 
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El momento angular total es L = r (0) x m p v (0) + 
r t (0) x mgViCO) , así el postulado de Bohr, se refiere a 
la versión moderna del postulado 4, la magnitud del 
momento angular total es un múltiplo entero de la 
constante de Planck dividida entre dos n, es: 



(1 + f)d 1 m e v 1 = nh 



(4) 



Donde n es un número entero positivo y h es la llamada 
constante de Planck dividida entre dos n. Esta relación 
establece la cuantización de la posición. Una vez que se 
conoce la velocidad, la posición solo puede tomar ciertos 
valores discretos de acuerdo a la ecuación 4. 

LA APROXIMACIÓN ELECTROSTÁTICA 

En esta sección asumimos que el campo eléctrico se puede 
calcular usando el campo electrostático o Coulombiano y 
despreciamos el campo magnético. Este aproximación se 
aplica en átomos ligeros, donde las velocidades esperadas 
de los electrones son mucho menores que la velocidad de 
la luz. De las ecuaciones 1 y 2 calculamos la aceleración 
inicial del electrón y del protón: 



%(o) = 



ri(0)-r (0) 



47re m e |f 1 (0)-r (0)| 3 



(5) 



3 



entonces decimos que la órbita es estable, y si se satisface 
la ecuación 4 entonces se clasifica como una de las órbitas 
de Bohr. 

Aí=1.5xl0" 19 s, ¿=1.5xl0" 15 s, v=2.2 Mm/s,T=1.32xlO" 17 s 
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Figura 1 . Trayectorias calculadas para el núcleo (color rojo), 
la cual ha sido aumentada 200 veces, y para el electrón 
(color azul). Las flechas indican la dirección del movimiento. 



a (0) = - 



?o(0)-?i(0) 



4nE m p \r (0)- 7^(0)1 3 



= +fa e t 



(6) 



Donde £ es la constante dieléctrica del vacío y 
a e = e 2 /4n£ m e dl(i + f) 2 . Hay que recordar que el campo 
eléctrico Coulombiano es proporcional al producto de 
las cargas e inversamente proporcional al cuadrado de la 
distancia que separa las partículas, su forma matemática 
es: 



E 1 {t) = - 



Ze 2 fiOQ-foCt) 
47T£ [r 1 (t) -r WI 3 



EL ALGORITMO NUMÉRICO 

Una vez que se tiene la aceleración al tiempo inicial t = o 
, se puede calcular las velocidades y las posiciones del 
electrón y del núcleo al tiempo t 1 = At, suponiendo un 
movimiento con aceleración uniforme y así sucesivamente. 
En un movimiento con aceleración uniforme la 
posición y la velocidad de la partícula están dadas por 

= x + v xo t + o.Sa x t 2 y v x (t) = v xo +a x t, y similarmente para 
la componente y. Si suponemos que esto se cumple para 
tiempos muy pequeño, entonces es posible usar estas 
ecuaciones para conocer las posiciones y las velocidades 
de las partículas involucradas a un tiempo posterior 
para así calcular las nuevas aceleraciones dadas por las 
ecuaciones 5 y 6. 

Si la órbita es estable, es decir, si ante una pequeña 
perturbación energética la órbita es casi la misma, 



De la figura 1 contando los puntos para completar 
una vuelta y teniendo el intervalo temporal para el cálculo 
de éstos, se obtuvo el periodo de oscilación t=1.32x10 17 
s y por tanto la frecuencia, el recíproco del periodo, es 
v=7.5x10 +16 Hz. 

Para el caso del átomo de Hidrógeno se puede 
demostrar que el cuasi-electrón, la pseudo-partícula 
que aparece al cambiar a las coordenadas del dentro de 
masa, satisface r 2 cp = cte, donde r y y son la magnitud 
de la posición y la velocidad angular del cuasi-electrón, 
respectivamente. Esto significa que el momento angular 
se conserva. Este resultado establece una ecuación 
para la magnitud r, la cual es resuelta en los libros de 
mecánica clásica (4). Este caso sirvió para probar el 
algoritmo numérico. En este problema se establece que las 
trayectorias pueden ser circulares, elípticas, etcétera. En el 
caso numérico que tratamos en esta sección, buscamos las 
condiciones numéricas para obtener trayectorias circulares 
y que sirviera como punto de partida en el siguiente caso: 
El átomo de Helio. En la siguiente sección analizaremos el 
átomo de Helio. 

EL ÁTOMO DE HELIO 

Dos electrones orbitando alrededor de un núcleo, que 
consta de dos protones y dos neutrones, forman un 
átomo de Helio. Este elemento es el segundo más ligero 
y el segundo más abundante del Universo. Al igual 
que en el caso anterior, se plantean las ecuaciones de 
movimiento para las tres partículas, enseguida se asume 



Raúl Gordo Hornos: El modelo otómico de Bohr: uno oplicoción . 



EPISTEMUS 




la aproximación electrostática y se desprecia la fuerza 
magnética. Posteriormente, se considera que a intervalos 
de tiempos pequeños el movimiento de cada una de las 
partículas es con aceleración uniforme. Estas ecuaciones 
permiten obtener la posición y la velocidad de cada una de 
ellas y de nuevo se regresan a las ecuaciones de movimiento 
para conocer la aceleración y así sucesivamente. Se buscan 
trayectorias estables que obedezcan los postulados Bohr. 



' s, t=5.12xl0" 17 s, v=2.50 Mm/s,T=3.91xlO +16 Hz 




Figura 2. Trayectorias del núcleo (color negro), con 
las condiciones iniciales seleccionadas el núcleo está 
estacionario, para los electrones ligados (color magenta) y 
(color naranja). Tanto la posición inicial de los electrones 
como sus velocidades se seleccionan del tal manera 
que satisfagan la condición de Bohr, así como la energía 
experimental del estado baso del átomo de Helio. 

Como se puede observar en la figura 2, las trayectorias 
de los electrones son elípticas y no circulares como en la 



figura 1. Esto puede ser así porque no se conoce la forma 
exacta teórica del estado base del átomo de Helio y esta 
se aproxima como la suma de dos átomos hidrogenoides. 

Cabe mencionar que con el paso del tiempo surgieron 
otros modelos, de tal manera que a finales de la década de 
1 920 se tenían cuatro teorías para explicar los espectros de 
emisión y absorción del átomo de Hidrógeno: La teórica 
semiclásica de Bohr, la teoría relativista de Sommerfeld, la 
ecuación de Schródinger y la ecuación relativista de Dirac. 

CONCLUSIONES 

El modelo de Bohr fue el primero que de una manera semi- 
clásica explicó los espectros de emisión y absorción del 
átomo de Hidrógeno. Se planteó un algoritmo numérico 
para estudiar la dinámica de átomos complejos, el cual 
está basado en la teoría de Bohr y en el movimiento 
uniformementeacelerado para tiempos infinitesimalmente 
pequeños. Se aplicó el formalismo al átomo de Hidrógeno 
con el fin de validar los resultados y luego se usó para 
estudiar el movimiento de los electrones y del núcleo en 
el átomo de Helio. Las trayectorias encontradas en este 
último caso son elípticas. 

Como un comentario final, es necesario enfatizar que 
las trayectorias de los electrones y del núcleo son sólo 
dentro del contexto del modelo semi-clásico de Bohr. 
La mecánica cuántica, en su forma actual, no habla de 
trayectorias, sino de regiones donde puede encontrarse a 
las partículas. 
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ARMONIZACIÓN DE MODELOS DE 
CALIDAD DEL SOFTWARE 



GABRIEL ALBERTO GARCÍA MIRELES*, 
FÉLIX ÓSCAR GARCÍA RUBIO 



La tendencia actual entre las organizaciones que desarrollan software es la adopción de más de un modelo 
de calidad para estar en condiciones de competir en diversos sectores del mercado. En esa situación, las 
organizaciones se enfrentan al problema de coordinar la implantación de múltiples modelos de tal manera 
que cada modelo logre los beneficios esperados, y al mismo tiempo, que los costos y esfuerzo de implantación 
sean mínimos. Como respuesta a este problema, surge la armonización de modelos de calidad del software. 
Actualmente los estudios se centran en la composición de modelos de calidad que abordan prácticas del 
proceso de desarrollo de software. Sin embargo, las organizaciones también deberían considerarla perspectiva 
de la calidad del producto. En este artículo, se presenta un avance de la investigación en la armonización de 
modelos de calidad que toma en cuenta ambas perspectivas, proceso y producto. Se describe la técnica que 
se propuso para realizar comparaciones de bajo nivel en este contexto. La técnica se aplica en los modelos 
CMMI-DEV, ISO 12207 e ISO 250 10 y, como resultado, se encuentra que los modelos que abordan la perspectiva 
de proceso contienen prácticas, que corresponden a las etapas de análisis y diseño de software, en donde se 
debería prestar más atención a la calidad del producto de software. 
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INTRODUCCION 

Las organizaciones que desarrollan software deberían 
consideraren susestrategiasdecrecimiento las necesidades 
del mercado para ser competitivas. Las necesidades de 
los clientes y de los usuarios pueden incidir en mejorar 
la calidad del software, liberar la versión del software en 
menor tiempo, o que el software se elabore dentro de los 
límites del presupuesto acordado. Una aproximación a la 
satisfacción de dichas necesidades se lograría a través de 
la evaluación del proceso empleado para desarrollar el 
software, y tras los resultados del diagnóstico, planear una 
iniciativa de mejora. 

Desde los años noventa del siglo XX, han estado 
surgiendo diversos modelos de calidad en los cuales 
se recopilan prácticas que los expertos consideran 
esenciales y que pueden contribuir a incrementar la 
calidad del software. Dichos modelos se fundamentan en 
los principios de la gestión de la calidad total (TQM, por 
sus siglas en inglés) en donde el enfoque hacia la calidad 
cubre todas las áreas de la organización, desde los niveles 
directivos hasta los ingenieros individuales y en donde 
todas las actividades que se realizan en el ciclo de vida 
del software tienen impacto en la calidad final obtenida. 
Estos modelos de calidad se basan en la convicción de que 
el proceso usado para desarrollar el producto de software 
influye en la calidad del producto y que un buen indicador 
de la calidad del producto se puede obtener si se conocen 
las prácticas o actividades que la organización lleva a cabo 
(1). 

Los modelos de calidad pueden presentarse en 
diferentes modalidades y con distintos requisitos de 
conformidad. Los modelos de calidad usados en procesos 
de evaluación o certificación son prescriptivos, es decir, 
establecen las prácticas o requisitos del proceso que 
tienen que ser adoptados en una organización o proyecto 



concreto. Algunos ejemplos de éstos son los estándares y 
los modelos de madurez. En el primer caso, los estándares 
contienen las prácticas recomendadas, enunciadas como 
requisitos, que deben ser implantadas en el proceso de 
desarrollo de software. Estos requisitos pudieran servir 
de base para establecer contratos entre un cliente y un 
proveedor. Modelos de calidad representativos de este 
grupo son el ISO 9000 o el ISO 1 2207. 

Los modelos de madurez, por su parte, presentan 
un mapa de ruta con el cual una organización puede 
determinar el grado en el cual su proceso de desarrollo 
es efectivo para lograr las metas del negocio, introducir 
innovaciones tecnológicas y usar los recursos limitados 
de la empresa eficientemente. Generalmente tienen de 
5 a 6 niveles de capacidad, y a partir de una evaluación, 
se determina el nivel que le corresponde a la empresa 
evaluada. Capacidad se refiere al grado en el cual el 
proceso de desarrollo de software se puede categorizar 
en distintos niveles, según se gestiona, define, administra 
cuantitativamente y optimiza. Modelos como CMMI e 
ISOI 5504 pertenecen a esta categoría. 

No obstante en los últimos años se asiste a una gran 
proliferación de modelos de calidad, lo que requiere que 
expertos analicen sus características para entender sus 
similitudesydiferenciasy tener la capacidad de recomendar 
el más apropiado para una empresa particular. Así, en 
la revisión realizada en (2), los autores han encontrado 
más de 315 estándares, guías y otros documentos 
prescriptivos que son mantenidos por 46 organizaciones 
diferentes. Otra experta analiza el tipo de dependencia 
que existe entre los modelos de calidad y para ello usa 
más de 35 modelos. Todos están enfocados en mejorar la 
calidad del software pero con distintos matices, según las 
organizaciones que los crearon y el dominio de aplicación 
destino. Sin embargo, tal cantidad de modelos complica la 
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toma de decisiones respecto del modelo de calidad más 
conveniente, y una mala decisión puede resultar en altos 
costos y esfuerzo. En el caso de los modelos de madurez, el 
esfuerzo de implantación puede llevar años para pasar de 
un nivel de capacidad a otro. Además cada modelo tiene 
un vocabulario, semántica, estructura y nivel de detalle 
distintos y con requisitos de evaluación diferentes. 

De hecho, la tendencia vigente es que las 
organizaciones, para acceder a distintos mercados, 
implanten simultáneamente distintos modelos de 
calidad (Figura 1). Esto hace que surja el problema de la 
implantación de múltiples modelos. En este estado, aunque 
se pretenden lograr los beneficios esperados al implantar 
cada modelo, la organización debería ser consciente de las 
limitaciones de recursos, tanto humanos como financieros, 
para emprender una iniciativa de esta naturaleza y de los 
conflictos en la coordinación entre múltiples modelos de 
calidad. Entonces, ¿Qué se puede hacer para lograr los 
beneficios que señalan los modelos y al mismo tiempo, 
reducir los costos de implantación y evaluación?, ¿Cómo 
afrontar los cambios que se dan en los estándares como 
resultado de su evolución inherente?, ¿Cómo mantener 
la alineación con los procesos implantados en la 
organización?, y ¿Cómo se pueden reutilizar los productos 
generados y la experiencia obtenida de la implantación de 
un modelo para reducir costos en la implantación de uno 
nuevo? 





Figura 1. Modelos de calidad más comunes que las 
organizaciones toman en cuenta para implantarlos. 

EL ENFOQUE DE LA ARMONIZACIÓN 

El marco de armonización surge como un enfoque para 
abordar la implantación simultánea de modelos de 
calidad. En este marco, se reconoce que las empresas 
abordan distintos estándares de calidad con la intención 
de satisfacer las metas del negocio. Por tanto, la empresa 
debe escogerlo de acuerdo a sus metas a largo plazo. Así, 



algunos expertos en el área (3) elaboraron un modelo 
de evaluación estratégica para alinear las metas de la 
compañía con los objetivos que persigue cada una de 
los modelos de calidad. Además, el marco considera que 
la introducción de nuevos modelos de calidad se realice 
coordinadamente con los modelos implantados con 
anterioridad. Otros aspectos fundamentales, considerados 
en el marco de la armonización, están relacionados con 
la selección de la tecnología y composición de la misma, 
así como el soporte a la evaluación y despliegue de las 
tecnologías. Una de las áreas que recientemente se está 
investigando es la composición de modelos, con especial 
atención a la comparación de modelos (4). 

La comparación de modelos se puede dar en un 
nivel alto de abstracción y en uno bajo. En el nivel alto se 
consideran las características de los modelos, su estructura, 
criterios de evaluación, objetivos y dominio. Con esos 
elementos se construyen taxonomías (o clasificaciones) 
de modelos de calidad con las cuales las organizaciones 
pueden identificar los modelos de calidad más apropiados 
a sus metas de negocio. 

Las comparaciones en un nivel de abstracción bajo o 
mopping, considera de cada modelo: metas, vocabulario, 
estructura, nivel de detalle de los elementos del modelo y 
requisitos de evaluación. Con base en esta información se 
determina el nivel de la estructura del modelo en el cual se 
debe hacer la comparación. Cuando una organización elige 
otro modelo de calidad para incorporarlo, es recomendable 
que conozca el efecto que tendrá sobre los modelos de 
calidad previamente implantados. Una manera de conocer 
el impacto es a través del mapping entre modelos al nivel 
que resulte adecuado para propósitos de certificación/ 
evaluación. Es en esta área en la que se concentra mucho 
del trabajo de la armonización. Por ejemplo, se tiene el 
marco de armonización, y diferentes comparaciones entre 
estándares. Los modelos más evaluados son el CMMI, ISO 
9001, ISOI 5504, ISOI 2207 (5). 

En el estado actual del avance de las investigaciones, 
la armonización se ha enfocado en la composición de 
modelos de calidad que abordan prácticas que deben 
implantarse. Es decir, la perspectiva de procesos en los 
trabajos de armonización es la dominante. Muchos de 
los estudios han considerado los modelos ISO9001, CMM, 
CMMI, ISOI 5504, SCAMPI y SPICE (5). 

Sin embargo, también es de gran importancia abordar 
la calidad de proceso teniendo en cuenta las características 
de calidad de producto y por tanto abordando de forma 
armonizada ambas perspectivas. 

PROPUESTA PARA LA ARMONIZACIÓN ENTRE 
ENFOQUES DE CALIDAD ORIENTADOS AL 
PRODUCTO Y AL PROCESO 

Como se mencionó en los primeros párrafos, una de los 
principios que subyacen en los modelos de calidad es la 
convicción de que la calidad del software depende del 
proceso empleado para desarrollarlo. Dentro del campo 
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de la ingeniería del software, se reconoce también la 
perspectiva de calidad orientada al producto (6). En esta 
se consideran las características de calidad que debe tener 
el software, tales como usabilidad, fiabilidad, rendimiento 
o seguridad. Esas características de calidad se han 
presentado en los modelos de McCall, ISO 9126, y el ISO 
25010, entre otros. 




La mayoría de los modelos de calidad orientados al 
producto de software tienen una estructura jerárquica. 
Las categorías de calidad contienen subcategorías y estás 
pueden estar asociadas tanto a propiedades identificables 
en los productos de software como a las opiniones y 
juicios de los usuarios. El modelo de calidad propuesto 
en ISO 25010 (7) identifica las siguientes características: 
funcionalidad, seguridad, usabilidad, mantenibilidad, 
portabilidad, rendimiento, compatibilidad y fiabilidad. 
En cada una de ellas se incluyen varias subcaracterísticas. 
Además, el modelo ISO 25010 considera que el software 
puede ser estudiado en dos etapas diferentes del ciclo 
de vida: desarrollo y operación. Durante el desarrollo se 
pueden especificar y evaluar las propiedades estáticas y 
dinámicas del software y, en la etapa de operación, seevalúa 
la calidad en uso del software cuando éste interacciona 
con los usuarios finales y otros sistemas del contexto de 
operación. Este modelo se puede usar para especificar y 
validar que los requisitos de calidad relevantes para los 
distintos usuarios hayan sido considerados en el proyecto. 

Aunque la relevancia de esta perspectiva de calidad 
del producto es notoria, poco trabajo se ha realizado desde 
la perspectiva de la armonización. De hecho, la mayoría de 
los trabajos se enfocan en proponer la combinación de 
modelos orientados al proceso (5). Se sabe que el objetivo 
principal de los modelos de calidad orientados al proceso 
incluye prácticas que la organización debe implementar. 
Por su parte, los modelos de calidad orientados al producto 
tienen una visión de las características del software 
relevantes para los usuarios de dicho producto. Además, 
algunos han considerado que proceso y producto son dos 
dimensiones distintas. 

Algunos investigadores han tratado ese problema 
empleando mopping para identificar los procesos que 



enuncian atributos de calidad y con ellos construyen 
árboles de decisión para recomendar modelos de calidad 
orientados al proceso. En nuestra investigación, estamos 
considerando la perspectiva de la armonización entre 
estándares de calidad de producto y de proceso para 
encontrar puntos de enlace (8). Dichos puntos de enlace 
son tareas (o actividades) que aparecen en los modelos 
de calidad enfocados en el proceso y que enuncian algún 
atributo o característica de calidad del producto o término 
relacionado (como requisito no funcional o requisito de 
calidad). Para realizar esta comparación, se toma como 
modelo base el modelo de proceso. Sobre éste se buscan 
los términos asociados al modelo de calidad de producto 
que nos interesa. La propuesta de mopping que se presenta 
puede identificar el nivel de la estructura en el cual se 
consideran las características de calidad. Esto último es 
relevante para establecer los requisitos de conformidad 
cuando se solicite la evaluación de los modelos de calidad 
implementados. 

El método que se propone incluye una actividad de 
análisis de los modelos en los cuales se identifican las 
metas, estructura y requisitos de los modelos de calidad 
que nos interesan. En la actividad de diseño del mopping 
se elabora el procedimiento de comparación de bajo 
nivel, según la estructura de los niveles y los requisitos 
que establece cada modelo. Tras definir el diseño, se 
efectúa el mopping en donde se extrae de cada modelo 
los enunciados relevantes que se han considerado como 
puntos de enlace. Finalmente, se evalúan los resultados 
para realizar el informe correspondiente. 

APLICACIÓN DEL MÉTODO DE 
ARMONIZACIÓN PROPUESTO 

Para mostrar cómo se pueden combinar modelos de 
calidad orientados al producto y proceso, realizamos dos 
ejercicios de mopping de bajo nivel entre los modelos 
de calidad. El primero de ellos fue entre el modelo de 
madurez de capacidades integrado (CMMI-DEV1.3) (9) 
y el ISO 25010 (7). El segundo incluye a los procesos de 
técnicos del ciclo de vida del software del ISOI 2207:2007 
(1 0) contra el ISO2501 (7). Tanto el CMMI-DEV e ISOI 2207 
corresponden a la categoría de modelos de calidad 
orientados al proceso y sólo consideramos los procesos 
principales para el desarrollo de software. De los modelos 
de calidad propuestos en ISO 25010, tomamos el modelo 
de calidad de producto de software y consideramos todas 
las características y subcaracterísticas del mismo. 

En este análisis buscamos elementos en los modelos 
que abordaran las características de calidad. Por ejemplo, 
la tabla 1 muestra un extracto de los enunciados que 
encontramos en los modelos. En la parte superior se señala 
un ejemplo del CMMI-DEV; mientras que en la parte inferior 
se muestra un ejemplo del ISO 12207. En ambos casos, 
se señala la localización según la estructura del modelo 
y cómo fue calificada su importancia respecto de lo que 
el modelo de procesos considera importante. El extracto 
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de CMMI-DEV se califica como relevancia baja porque el elemento que lo 
contiene es una nota informativa y estos elementos sólo describen cómo 
podría implementarse la práctica, mas no son obligatorios. En el caso del 
extracto del ISO 1 2207, al ser un requisito, es un elemento de alta relevancia 
que la organización que decida usar ISO 12207 deberá también incluir en 
su proceso de desarrollo prácticas para garantizar que las características de 
calidad sean abordadas adecuadamente en el proyecto. 

En la tabla 2 se muestran los diferentes enunciados que encontramos 
según la estructura del CMMI-DEV. Como se observa, la mayoría de 
las referencias a atributos de calidad tiene una relevancia baja, pues 
se encuentran enunciados en subprácticas, elementos informativos y 
ejemplos. Sólo encontramos una referencia a características de calidad 
en una práctica específica del área de proceso de desarrollo de requisitos. 
En el nivel de metas no hay mención explícita a características de calidad 
del producto. Además, se puede observar que la mayoría de los puntos de 
enlace identificados en las áreas de proceso denominadas desarrollo de 
requisitos y solución técnica corresponde a prácticas que deben realizarse 
durante la identificación de requisitos, análisis de requisitos y el diseño de 
la arquitectura del software. Pocos puntos de enlace se identifican en las 
áreas restantes (validación, verificación e integración del producto). En el 
caso del ejercicio con ISO 12207 también encontramos que los puntos de 
enlace se concentran en los procesos de definición de requisitos, análisis de 
requisitos y diseño arquitectónico del sistema. 



Tabla 1. Extracto de los puntos de enlace identificados en los modelos de 
calidad enfocados en los procesos. 



Modelo 


Elemento 


Descripción 


Relevancia 


CMMI-DEV 


Nota informativa 
de la práctica 
SP1.3 del área 

de proceso 
Integración del 
producto 


Los criterios pueden 
ser definidos por el 
comportamiento 
del rendimiento 
(funcionalidad y 
atributos de calidad) 


Baja 


ISO12207 


Requisito 
6.4.2.3.1.1 


La especificación 
de requisitos debe 
describir: ... requisitos 
de usuario, seguridad, 
ingeniería de factores 
humanos, requisitos de 
interfaz, operaciones, y 
mantenimiento 


Alta 



Tabla 2. Punto de enlace en CMMI-DEV. 



Área de proceso 


Relevancia 
alta 


Relevancia 
media 


Relevancia 
baja 


Validación (VAL) 








1 


Verificación (VER) 








1 


Integración del 
producto (Pl) 








4 


Solución técnica (TS) 








18 


Desarrollo de 
requisitos (RD) 





1 


22 
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Es importante recordar que los modelos de calidad 
orientados al proceso incluyen prácticas para que la 
organización elija los métodos adecuados para el tipo 
de software que desarrolla. Por otra parte, los estándares 
generalmente contienen los requisitos mínimos que 
se esperarían en el desarrollo de software. Bajo este 
contexto, la identificación de puntos de enlace es una 
primera aproximación al problema de mejorar la calidad 
del software combinando ambas perspectivas. Por una 
parte, se requiere investigar con mayor profundidad 
qué prácticas son relevantes para incrementar las 
características de calidad que los usuarios buscan en un 
producto de software particular. Además, se debe estudiar 
cómo se pueden alinear estos modelos de calidad con la 
práctica real del desarrollo de software. 

CONCLUSIONES 

En este artículo se ha presentado brevemente el área 
de armonización de modelos de calidad. Dado que 
en la armonización se pretende introducir de manera 
coordinada nuevos modelos de calidad, las organizaciones 
que deben satisfacer los requisitos de conformidad de 
distintos modelos podrían beneficiarse de los métodos 
y técnicas que existen, y que se están desarrollando. Con 
base en esta sugerencia, las empresas pueden optimizar 
el uso de recursos y lograr los beneficios esperados que 
cada modelo de calidad señala. Además, el proyecto de 
investigación que se está realizando busca combinar 
las dos perspectivas de calidad, producto y proceso, a 
través de enfoques de armonización. Estos resultados 
iniciales permiten identificar áreas que necesitan soporte 
metodológico para emprender una iniciativa de mejora 
integral de la calidad del software. 
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DIAGNÓSTICO DE LOS INSTRUMENTOS UTILIZADOS 
PARA LA MEDICIÓN DE LA PRESIÓN ARTERIAL EN LOS 
PACIENTES ATENDIDOS EN UN CENTRO MÉDICO 
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En los sistemas de salud se dispone de una diversidad de instrumentos, procedimientos, pacientes y 
competencias técnicas del personal que realiza mediciones. Los elementos pueden tener una influencia 
significativa en los datos obtenidos pudiendo presentar variabilidad con respecto a las condiciones reales de 
salud de los pacientes, lo que se puede reflejar en información incorrecta, influyendo en la toma de decisiones 
equivocadas al momento de diagnosticar y proporcionar terapias a los pacientes. Lo anterior puede deteriorar 
la salud de los pacientes atendidos y, por consecuencia, su calidad de vida. Este artículo muestra los resultados 
de una auditoría administrativa de las condiciones metrológicas que deben tener todos los instrumentos de 
medición de acuerdo a la LFMN, así como de sus características físicas de acuerdo a la normatividad mexicana. 
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INTRODUCCIÓN 

Los sistemas de medición son el conjunto de 
procedimientos, instrumentos, personal y condiciones 
ambientales necesarias para realizar una medición. 
Por lo tanto, el control de estos factores son de gran 
importancia, porque permite tener confiabilidad de los 
datos obtenidos; se debe garantizar que estos datos de 
medición representan las condiciones reales de salud de 
los pacientes. Así pues los profesionales de salud deben 
utilizar un sistema de medición confiable que les permita en 
la medida de su competencia tomar decisiones correctas, 
logrando asegurar que los diagnósticos y tratamientos 
proporcionados se encuentren en concordancia con la 
condición de salud de los pacientes atendidos. 

En el Plan Nacional de Salud de México (1 ) se establece 
que la protección de la salud de los mexicanos requiere de 
estrategias integrales que fortalezcan y amplíen la lucha 
contra los riesgos sanitarios y favorezcan la cultura de la 
salud y el desarrollo de oportunidades para elegir estilos 
de vida saludables, mediante la prestación de servicios de 
salud con calidad y seguridad. 

Las mediciones de la presión arterial, temperatura 
corporal y la frecuencia de los latidos del corazón en los 
pacientes son mediciones realizadas por los profesionales 
de salud de forma rutinaria y son consideradas como 
parámetros que determinan la condición de salud de los 
pacientes. Desde la antigüedad, la presencia del pulso 
arterial se ha entendido como un signo fundamental de la 
vida (2). 

Adicionalmente, la caracterización de parámetros 
de la composición de la sangre y orina, al igual que las 
mediciones eléctricas que sirven para determinar el estado 
del corazón, del cerebro y de los músculos, son prácticas 
comúnmente empleadas por los profesionales de salud 
en el diagnóstico y posterior tratamiento de los pacientes, 
donde el diseño e ¡mplementación de los sistemas de 
aseguramiento de las mediciones generan confiabilidad 
de los datos obtenidos. También reflejan la información de 
salud real de los pacientes, permitiendo a los profesionales 
de la salud la aplicación de sus conocimientos para la 
generación de una calidad de vida saludable a los pacientes 
atendidos en los diferentes sistemas de salud. 

En el Foro de Metrología para el Sector Salud (3) se 
identificó la necesidad de establecer laboratorios de 
calibración específicos para el área médica a nivel nacional. 
Esta necesidad se hace más evidente ya que, de acuerdo al 
Departamento de Comercio de los EUA, la comercialización 
y uso de instrumentos de medición rehabilitados o de 
segunda mano en clínicas y hospitales en México se está 
incrementando. Sin embargo, no se tiene conocimiento 
de que estos equipos sean evaluados metrológicamente 
antes de ser utilizados en servicios médicos. Además 
señalan la necesidad de establecer mecanismos que 
permitan asegurar la calidad y evaluar los diferentes 
instrumentos de medición utilizados para brindar servicios 
de salud, ya que las mediciones realizadas con estos 



instrumentos tiene una gran influencia directamente en 
la calidad del diagnóstico y/o efectividad del tratamiento 
proporcionados. 

PLANTEAMIENTO DEL PROBLEMA 

Para efecto de esta investigación se verificaron las 
característica metrológicas de los instrumento utilizados 
para la medición de la presión arterial, esfigmomanómetros 
(Figura 1). La información proporcionada por las 
mediciones realizadas con estos instrumentos debe 
reflejar la condición de los parámetros de salud de los 
pacientes y debe permitir a los profesionales de salud 
tomar decisiones confiables en los diagnósticos y 
tratamientos necesarios para cubrir las necesidades de 
salud de los pacientes atendidos. Los servicios de salud 
tienen una posición única entre otros servicios debido a 
la propia naturaleza de los riesgos altamente involucrados 
(4), los cuales generalmente puede ocasionar la pérdida 
de la vida o el deterioro de la salud de los pacientes, lo que 
genera una disminución de su calidad de vida. 




Figura 1. Esfigmomanómetro. 



Los sistemas de salud deben contar con instrumentos, 
procedimientos de medición estandarizados y programas 
de aseguramiento de la calidad de las mediciones de 
acuerdo con las características de los pacientes. Las 
fuentes de variabilidad en una medición podrían deberse al 
mismo procedimiento (instrumentos de medida, factores 
ambientales), al observador (metodología de medición, 
déficit sensitivos) y al sujeto observado (5). El presente 
artículo tiene como objetivo general identificar las 
condiciones metrológicas de los instrumentos utilizados 
para la medición de la presión arterial, lo cual permita 
definir las estrategias que nos permitan incrementar la 
confiabilidad de las mediciones realizadas. 

Es importante que todos los instrumentos utilizados en 
cualquier medición cuenten con la información visible y las 
características físicas que demuestren que sus resultados 
de medición son trazables y confiables, como se establece 
en el artículo 1 1 de la LFMN (6). 
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JUSTIFICACIÓN 

Cuando se toman en cuenta los datos obtenidos a partir 
de una medición no confiable, se pueden presentar dos 
tipos de errores en el análisis de la información generada, 
los cuales son: 

DATOS FALSOS POSITIVOS: Estos datos generan 
diagnósticos erróneos cuando se determina que el paciente 
presenta hipertensión y en realidad no la tiene. Esto 
implicaría proporcionar al paciente medicamentos que 
no requiere, pudiendo ocasionar otros padecimientos de 
salud, así como un alto costo al suministrar medicamentos 
no necesarios. 

DATOS FALSOS NEGATIVOS: Estos datos generan una 
información ala especificar al paciente con un diagnóstico 
de sano cuando en realidad es hipertenso, lo cual implicaría 
que, al no proporcionarle los tratamientos adecuados, los 
pacientes pueden presentar una enfermedad degenerativa 
que afecte su calidad de vida. 

El Centro Nacional de Excelencia Tecnológica en Salud, 
CENETEC (7), estableció que la importancia de la correcta 
medición de la presión arterial reside en el hecho de que 
con un error sistemático de medición, como podría ser 
la subestimación de 5 mm/Hg de la presión arterial real, 
podría traducirse en que 21 millones de personas no 
están recibiendo tratamiento antihipertensivo; lo anterior 
significa que serán candidatos a la presentación de 
complicaciones por hipertensión arterial. Por el contrario, 
un error sistemático de medición de 5 mm/Hg por arriba 
de la presión arterial real, podría clasificar erróneamente 
a 27 millones de personas como hipertensas cuando no lo 
son y exponerlos a los riesgos que implica la prescripción 
de medicamentos antihipertensivos a personas que no lo 
requieren. 

La hipertensión es una de las enfermedades más 
comunes en todo el mundo que afligen a los seres humanos 
debido a la morbilidad y mortalidad. La medición de la 



presión arterial es el método principal para la detección 
de la hipertensión en los pacientes y de los padecimientos 
correlacionados. 

En México, de acuerdo a la estadística de mortalidad 
de 2008 elaborada por el Sistema Nacional de Información 
de Salud (8) las causa de mortalidad por enfermedades 
hipertensivas se encuentran en el séptimo lugar de las 
causales con un 2.9 % (15,698 defunciones). La primera 
causa de mortalidad es la diabetes mellitus con un 14,0 
% (75,572 defunciones), caracterizada por el incremento 
de padecimientos del corazón, cerebro vasculares y 
la insuficiencia renal. En el año 2000, la prevalencia de 
diabetes fue de 7,5% y la de hipertensión fue de 30,7%; 
así mismo casi la mitad de los adultos con diabetes 
tenían hipertensión (9). Actualmente el 8 % de los adultos 
mayores a veinte años la padecen (1 ). La segunda causa de 
mortalidad son las enfermedades isquémicas del corazón, 
con un 11.1 % (79,579 defunciones) que se caracteriza 
por una disminución de aporte de oxígeno al corazón por 
consecuencia de obstrucción y/o estrechamiento de las 
arterias coronarias. La tercera causa es por enfermedades 
cerebro vasculares con un 5.6 % (30,212 defunciones), de 
un total de 538,288 defunciones registradas en México 
durante el 2008. 

La diabetes y las enfermedades isquémicas del corazón 
y las cerebro vasculares son los principales factores de 
riesgo de muerte en los pacientes. Por lo tanto la medición 
de la presión arterial confiable nos proporciona una 
información que permite poder realizar diagnósticos 
adecuados y a tiempo a los pacientes y evitar los riesgos 
a la salud provocados por estos padecimientos. La 
hipertensión arterial también influye en el desarrollo de las 
enfermedades isquémicas del corazón y las enfermedades 
cerebro-vasculares (1). 

El funcionamiento adecuado de los establecimientos 
de atención médica está en relación directa con los 
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conocimientos, habilidades, destrezas, motivación y ética 
del personal de salud, para que utilice adecuadamente la 
tecnología que se pone a su disposición como parte de la 
infraestructura y equipamiento (1 0). 

Dentro de los servicios de salud, la capacidad de los 
prestadores de servicios de hacer juicios confiables es, 
de hecho, una habilidad esencial y crítica (11), por lo cual 
las mediciones de la presión arterial deben identificar 
las condiciones reales de los pacientes, para que los 
diagnósticos generados estén en concordancia con el 
estado de salud de los pacientes. 

METODOLOGÍA 

Para lograr los objetivos de la investigación desarrollada 
se determinaron las condiciones metrológicas de los 
esfigmomanómetros utilizados para la medición de la 
presión arterial en las diferentes áreas de un centro médico 
mediante la siguiente herramienta de auditoría. 

Auditoría de instrumentos de medición 

A fin de adquirir la información de interés se realizó una 
auditoria de los instrumentos de medición de la presión 
arterial existentes dentro de las instalaciones del centro 
médico, verificando el cumplimiento con las características 
metrológicas establecidas en la guía de verificación (Figura 
2). 



TIPO 

INSTRUMENTO: 



MARCA: 



SERIE: 



MODELO: 



CALIBRACIÓN: SI 



NO 



FECHA DE PUESTO 
EN SERVICIO: 



CONDICION: 



FECHA 

CALIBRACIÓN: 



EMPRESA: 



UBICACIÓN: 

DIMENSIÓN DE 
MANGUITO (cm): 

LOCALIZACIÓN: 



# MEDICIONES (Día): 

LARGO 



ANCHO 



PROGRAMA DE 
CALIBRACIÓN: 



SI 



NO 



# USUARIOS: 

# INSTRUMENTO: 



Figura 2. Guía de verificación de auditoría. 



RESULTADOS 

A continuación se muestran los resultados obtenidos 
mediante la aplicación de la auditoría de los instrumentos en 
el centro médico. Durante el desarrollo de la investigación 
se utilizó la Estadística Descriptiva como herramienta 
para el análisis de los datos por medio de la generación de 
gráficas, por ser una de las formas más rápidas y eficiente 
para la realización de los análisis correspondientes. 

El centro médico nos proporcionó el inventario de 
instrumentos de medición con los que cuentan en cada 
área especificada. De acuerdo a este inventario, que da un 
total de cincuenta instrumentos de medición, se realizó 
una auditoría de veinte instrumentos (40 %) de tipo 
monitor (Figura 3). 



ÁREAS 


CANTIDAD DE 
INSTRUMENTOS 


URGENCIAS 


6 


UCI ADULTOS 


2 


UTI ADULTOS 


4 


GINECOLOGÍA 


2 


PEDIATRÍA 


13 


MEDICINA 


6 


CIRUGÍA 


8 


UCIP 


3 


UCIN 


6 



Figura 3. Inventario de instrumentos de medición de la 
presión arterial por ubicación. 



Además de los instrumentos inventariados, se 
localizaron otros en las áreas mencionadas de tipo 
aneroides de pedestal y empotrados en la pared 
no inventariados, dando un cantidad de diecisiete 
instrumentos de medición; de igual manera se les aplicó la 
guía de verificación correspondiente. 

No se auditaron todos los instrumentos por 
encontrarse localizados en áreas de acceso controlado o 
por las condiciones de salud de los pacientes a los cuales 
se les realizaban las mediciones. 

Se presentan los resultados obtenidos a partir de 
la aplicación de la auditoría en donde se agrupó la 
información específicamente en dos características 
metrológicas normalizadas que deben cumplir los 
instrumentos de medición correspondientes a la 
trazabilidad de la mediciones (sistema de calibración de los 
instrumentos) y dimensiones de los manguitos (brazalete). 
Por consiguiente, la obtención de la información que 
nos permitiera un análisis de los resultados más precisa 
y una mejor determinación de las conclusiones y 
recomendaciones, correspondiente a los datos requeridos 
por la herramienta, se pueden observar a continuación. 
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Es importante remarcar que los resultados, de 
acuerdo a la estadística descriptiva utilizada, nos muestra 
la tendencia generalizada de que los instrumentos no 
presentan ningún control administrativo y necesario para 
demostrar la trazabilidad de las mediciones realizadas y 
que se puede propagar esta ¡neficiencia a los otros tipos 
de instrumentos de medición requerido en los sistemas de 
salud. 

Resultado 1 . Sistema de calibración 

Se encontró que los veinte instrumentos tipo monitor 
y los diecisiete tipo aneroide no contaban con ninguna 
referencia física (etiqueta de calibración) que demostraran 
que los instrumentos habían sido calibrados interna o 
externamente. Por consiguiente, no contaban con la fecha 
de calibración, nombre del laboratorio de calibración 
utilizado, así como tampoco un programa de calibración 
autorizado e implementado. 

La totalidad de los instrumentos auditados no 
contaban con ninguna información visible que demostrara 
la trazabilidad de las mediciones realizadas, lo cual 
demuestra que los datos obtenidos no son confiables, 
desde un punto de vista metrológico. 

Resultado 2. Dimensión de manguito 

El manguito es la cámara o bolsa de insuflación por 
compresión (parte que comprime la arteria externamente 
lo cual interrumpe el paso de la sangre formada por una 
cubierta inelástica pero flexible). Por las características 
de su función es importante que sus dimensiones y las 
características de los pacientes sean las adecuadas. 

Se encontró una gran diversidad de dimensiones 
de los manguitos en las diferentes áreas auditadas. Las 
dimensiones van desde 14.5 x 27.0 cm, que son los de 
dimensiones más grandes, hasta de 2.5 x 7.0 cm, los cuales 
son los más pequeños, (Figura 4). 

Comparando las dimensiones de los manguitos 
encontrados con las especificaciones establecidas con 
respecto al tipo de paciente o dimensiones del brazo 
establecido en la figura 5 para que los resultados obtenidos 
por las mediciones sean confiables. 



14.5 X 27 




Figura 4. Dimensiones de manguito. 

Las diferencias entre el tamaño del manguito y la 
dimensión del brazo pueden generar errores falso positivos 
o falso negativos según sea el caso, como se indica en el 
punto 3. 





Circunferencia 


Manguito 


del brazo 
(cm) 


Ancho 
(cm) 


Largo 
(cm) 


Neonato 


5.0-7.5 


3 


5 


Preescolar 


7.5-13.0 


5 


8 


Escolar 


13.0-20.0 


8 


12 


Adulto 


24.0 - 32.0 


13 


24 


Adulto grande 


32.0-42.0 


17 


32 


Muslo 


42.0 -50.0 


20 


43 



Figura 5. Tamaños de manguito. CENETEC. 



No se observaron manguitos con las dimensiones 
exactas establecidas. Las que más se acercan son las de 
los manguitos de 13 X 23 cm y 14 x 25 para ser utilizados 
por pacientes adultos con dimensiones del brazo entre 24 
- 32 cm que fueron un total de doce instrumentos (32%) 
y quince instrumentos (40%) respectivamente, utilizados 
en las áreas de urgencias, medicina, cirugía y pediatría. 
Los demás manguitos (28%) se encuentran fuera de las 
especificaciones establecidas por el CENETEC, por lo cual 
pueden presentar datos no confiables. 

CONCLUSIONES 

La aplicación de la metodología tiene como objetivo 
demostrar si las decisiones tomadas por el personal 
de salud en el diagnóstico y en los tratamientos 
proporcionados a los pacientes son realizadas mediante 
el uso de instrumentos con características metrológicas 
confiables. Es importante aclarar que estas conclusiones se 
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establecen desde un punto de vista metrológico, es decir 
apegados a las características metrológicas y normalizadas 
que deben tener todos los instrumentos de acuerdo a la 
LFMN (6) y la normatividad mexicana. 

Los resultados muestran, a partir de la aplicación de 
la herramienta, que existe suficiente evidencia obtenida 
a partir de la estadística descriptiva utilizada para el 
análisis de los datos, que las mediciones no son confiables. 
Los instrumentos auditados no proporcionan datos 
metrológicamente válidos, no se encontró evidencia de 
que los instrumentos se encuentren calibrados, lo cual es 
una condición para demostrar que los instrumentos no 
pueden proporcionar datos reales de este parámetro de 
salud. Es importante considerar que las conclusiones son 
establecidas desde un punto de vista metrológico y no 
médico. 

No se tiene un programa de calibración, se observaron 
instrumentos que no tenían número de serie, el cual es 
proporcionado por el fabricante, ni identificación única 
lo que dificulta conocer la rastreabilidad del instrumento 
cuando se considera importante. No se encontró 
información que indicara cuándo fueron puestos en 
servicio ni el número de mediciones realizadas por cada 
instrumento, por lo cual no se podría determinar la deriva 
del mismo. En general, todos los instrumentos tienen una 
condición buena en sus características físicas. Asimismo, se 
encontraron instrumentos fijos en la pared entre dos camas, 
los cuales no se utilizan por no tener accesorios apropiados 
para la distancia entre el paciente y el instrumento, 
además colocados en pedestales y no fijos lo cual pueden 
sufrir daños por el manejo inapropiado. Se encontraron 
dimensiones de manguitos de los instrumentos diferentes 
a las especificaciones establecidas por el CENETEC, que es 
el organismo responsable, dependiente de la Secretaría 
de Salud en México, del control de uso de este tipo de 
instrumentos de medición. 

RECOMENDACIONES 

De acuerdo a los resultados y a las conclusiones 
obtenidos, se recomienda diseñar un sistema de 
calibración con respecto a los instrumentos de medición, 
donde se elabore: 

Un procedimiento de calibración aplicable para todos 
los tipos de instrumentos utilizados donde se asegure 
la trazabilidad de los patrones utilizados. Se tiene que 
considerar que se cuenta en el centro médico con un 
instrumento patrón que se puede utilizar para la calibración 
de los instrumentos, previa calibración correspondiente, 
validación del procedimiento de calibración interno, 
acondicionamiento del espacio físico y capacitación de 
personal específico para el aseguramiento de calibraciones 
por desarrollar, así como el diseño de un programa de 
calibración interna donde se encuentren incluidos todos 
los instrumentos utilizados en las áreas correspondientes. 

Un procedimiento de identificación única de cada 
instrumento, a fin de que permita rastrear su uso y poder 



determinar su estatus de confiabilidad en cuanto a los 
datos generados por las mediciones, así como el número 
de mediciones realizadas por cada instrumento y que 
éstos se puedan utilizar como indicativo para determinar 
la fecha de próxima calibración. 

Un procedimiento técnico mediante la comparación 
de los resultados obtenidos por las mediciones de los 
instrumentos de medición de diferentes marcas y modelos, 
lo cual permita identificar los instrumentos más robustos y 
confiables, es decir, que no presenten errores de medición 
significativos, y que esto influya en la medida de lo posible 
en la adquisición de nuevos instrumentos requeridos 
por las diferentes áreas mediante el establecimiento de 
criterios técnicos de aceptación y no precisamente por 
criterios administrativos. 
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DISEÑO DE ALTERNATIVAS DE CONSERVACIÓN DE 
CARRETERAS EN EL SUR DEL ESTADO DE SONORA 



ELSA ELIZABETH MORALES MORALES*, JOSÉ FRANCISCO GONZÁLEZ 
OCHOA, ISRAEL MIRANDA PASOS, ANA LILIA LEAL CRUZ 



En este trabajo se analiza el comportamiento de la estructura de pavimento existente, identificando las 
variables que intervienen en la realización de los estudios de evaluación, con la finalidad de proponer una 
alternativa de rehabilitación de dicha estructura. La infraestructura carretera existente se construyó a un costo 
elevado en décadas pasadas; el incremento del flujo vehicular y el no tener una adecuada programación de 
los trabajos de mantenimiento han dado como resultado que actualmente se cuente con una red de caminos 
deteriorada, un mayor costo de transporte y un impacto en general negativo para la economía. 

La identificación de las variables que intervienen en el proceso de evaluación de proyectos, así como los criterios 
que deben ser aplicados a cada una de ellas, conforman la base conceptual de toda metodología que tiene por 
objeto determinar una serie de indicadores de tipo económico. Estos indicadores constituyen una herramienta 
importante para que los responsables de la programación de inversiones de infraestructura carretera tomen 
decisiones para la optimización de los recursos disponibles. Una conservación adecuada de la carretera permitirá 
reducir los costos de operación de los vehículos que participan en el traslado de bienes y personas, además de 
incrementar la velocidad de operación que en términos económicos se traduce en una mayor rotación de capital 
en cualquiera de sus manifestaciones (dinero, fuerza de trabajo y mercancía). 
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INTRODUCCIÓN 

En el tramo de Guaymas - Estación Don, el cuerpo B se estima que se 
construyó en 1950, su estructura de pavimento consistía en una capa 
de terracería con una superficie de rodamiento de revestimiento. En 
1960 se modificó la estructura de pavimento, aprovechando la capa de 
revestimiento existente para que funcionara como una capa de base 
hidráulica y como superficie de rodamiento se tendió una capa de carpeta 
de mezcla elaborada en el lugar (mezcla en frío). Esta capa se deterioró por 
su uso y en 1 984 se construyó una capa de carpeta de concreto asfáltico 
elaborada en caliente, que es la existente en la actualidad. El cuerpo A 
se construyó en el período que va de 1985 a 1988, el cual consistía en 
una estructura de pavimento formada por una capa subrasante, base 
hidráulica y carpeta de concreto asfáltico. Esta información se obtuvo 
con base en la experiencia del personal que ha trabajado en el sector 
carretero. 

En ambos cuerpos se han realizado trabajos de conservación, 
como son: bacheos, riegos de sello, riegos de protección, carpetas de 
renivelación, estabilización de base con asfalto, estabilización de base 
con cemento portland, la cual se lleva a cabo para mitigar el deterioro de 
la superficie de rodamiento. 

Las características geométricas de la carretera en el subtramo del 
km 80+000 al km 106+000 señalan que es un camino tipo A4S (de dos 
cuerpos, con dos carriles de circulación por sentido y con faja separadora), 
con ancho de carril de 3.50 m, acotamiento a ambos lados de 1.50 m, 
con faja separadora con ancho variable de 3.00 m a 10.00 m, presenta 
también una pendiente máxima de 2.5%, bombeo del 2%, grado de 
curvatura máximo de 3.5° y se encuentra ubicado en terreno lomerío. La 
velocidad de proyecto es de 1 00 km/h. 

Con el análisis del subtramo en estudio se propondrán las alternativas 
de los trabajos de rehabilitación de la estructura de pavimento a fin de 
disminuir los problemas ocasionados por las fallas que se presentan 
en la superficie de rodamiento. De esta manera se tendrán beneficios 
en ahorros en tiempos de recorrido y se proporcionará seguridad y 
comodidad a los conductores. 

DESCRIPCIÓN DE LAS CARRETERAS EN MÉXICO 

Algunas carreteras que se construyeron hace más de sesenta años, para 
condiciones de tránsito radicalmente diferentes a las actuales, muestran 
debilidad estructural y problemas para su correcta conservación. En 
efecto, los vehículos de antaño transmitían esfuerzos relativamente 
pequeños, cuyo alcance vertical era también escaso, quizá no superaba 
los 30 ó 40 cm. En comparación, los vehículos de carga actuales producen 
esfuerzos mucho mayores, con valores significativos y profundidades 
más grandes, en el orden de 1 m e incluso más. 

A continuación se presentan dos casos. El primero es cuando 
los materiales empleados en aquellos años para la construcción, 
especialmente en terracerías, eran de una calidad que hoy debe 
considerarse como inaceptable para capas que quedan bajo la 
influencia de las nuevas cargas. En las carreteras del país, las capas de 
terracerías son arcillosas, de baja resistencia y muy sensibles a cambios 
volumétricos debido a la variación en sus contenidos de agua, lo que 
conduce, obviamente, a carreteras de superficie muy deformable. El 
segundo consiste en que el bajo volumen de los vehículos que entonces 
se consideraban pesados, producía efectos de fatiga relativamente poco 
notables. En la actualidad, esos materiales débiles están al alcance del 
efecto de penetración de los modernos arreglos vehiculares y dejan 
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ver dramáticamente su baja resistencia; pero además, la 
repetición de cargas mucho más pesadas induce efectos de 
fatiga devastadores y causan deformaciones permanentes 
intolerables. 




Estas condiciones imponen a la red básica mexicana 
condicionantes de conservación muy propia y, por 
supuesto, diferentes a las prevalecientes en otras redes 
carreteras en las que ya se ha realizado un esfuerzo de 
modernización que México aún no ha completado. Así 
pues, es necesaria una nueva estrategia de construcción 
de las carreteras para que se incorporen en el futuro a 
la red mexicana. Anteriormente, la filosofía de diseño 
de la sección estructural fue lograr una zona superior 
relativamente resistente, aceptando en forma progresiva 
materiales francamente débiles en la zona inferior, que 
se consideraban a salvo de la influencia de las cargas. 
Actualmente, para la conservación de dichas rúas, se 
requieren verdaderas acciones de reconstrucción en lo 
profundo, pues aquellas zonas débiles quedaron dentro 
de la zona crítica de esfuerzos. El cambio de diseño que 
se pondera para la época actual tiende a lo contrario: 
secciones convenientemente robustas en capas profundas. 
El pavimento no es una estructura que falle de un minuto 
a otro, capas superficiales débiles significan duraciones 
cortas, de manera que el criterio expuesto puede 
manejarse dentro de otro de inversiones diferidas. Fallas 
en lo profundo no se resuelven más que con costosísimas 
operaciones de reconstrucción. 

La red básica actualmente tiene una longitud de 30,000 
km que cubren un área de la geografía nacional en la que 
seguirán manteniéndose fuertes demandas de transporte. 
Esta red deberá de ser reforzada estructuralmente, 
ampliada y mejorada para afrontar demandas más 
exigentes. En muchos casos las vías actuales deben de ser 
sustituidas por otras nuevas más acordes con las crecientes 
demandas. 

Las grandes inversiones que exigen las carreteras 
modernas requieren una creciente participación del capital 



privado. El aspecto del financiamiento y de los mecanismos 
de recuperación correspondientes será cada vez más 
relevante. El punto álgido a buscar en estos mecanismos 
de financiamiento es el de lograr paquetes financieros que 
permitan disponer de recursos para proyectos en los que 
la recuperación económica será escasa en los primeros 
años de funcionamiento. Si la planeación nacional ha 
de adelantarse a la demanda o, por lo menos, ha de 
conducir a obras de aparición muy cercana a la necesidad, 
seguramente seguirá ocurriendo que carreteras muy 
importantes para el país comiencen con niveles de tránsito 
bajo pero susceptibles a desarrollos adecuados; de esta 
manera podrán afrontarse las primeras etapas de vida de 
las obras con cuotas o costos de mantenimiento razonables 
que no se transformen en un factor decisivo de su puesta 
en operación. Independientemente del desarrollo de 
una red productora de riqueza con capacidad estructural 
adecuada, buena cobertura y buena conservación, habrá 
de seguirse desarrollando una red alimentadora municipal 
y rural que, con sus propias funciones, redondee el 
panorama nacional del transporte carretero (1). 

UBICACIÓN DE LA ZONA EN ESTUDIO 

Se localiza en las cartas topográficas Obregón, Guaymas 
y Huatabampo, escala 1:250,000, limitada por las 
coordenadas del tramo: Estación Don - Cd. Obregón, en el 
km 68+900 en la latitud Norte 26°1 4' y en la longitud Oeste 
109°01', ASNM53 (Figura 1)(2). 




■I'JMpO B 



Figura 1 . Localización del tramo en estudio (SCT, 201 0). 

ANÁLISIS DE LA SUPERFICIE DE RODAMIENTO 
DEL TRAMO EN ESTUDIO 

Los estudios básicos para la evaluación de la superficie de 
rodamiento son: calificación o índice de servicio actual, 
levantamiento de deterioros, estudio geotécnico y pruebas 
de laboratorio de suelos, análisis del tránsito y diseño de 
los espesores de la estructura de pavimento. 
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Calificación o índice de Servicio Actual (ISA) 

La calificación que se otorga a un tramo de la vialidad 
se referirá siempre a toda la sección considerada, por lo 
que ésta debe reflejar el estado promedio que guarda 
el tramo dentro de la propia sección. Las condiciones 
de transitabilidad de la superficie de rodamiento se 
calificaron en función de la comodidad y/o seguridad que 
experimenta el usuario al transitar por el camino. Se utilizó 
una escala del al 5, en el cual el se considera como 
camino intransitable; del al 1 camino muy malo; del 1 al 
2 camino malo; del 2 al 3 camino regular; del 3 al 4 camino 
bueno; del 4 al 5 camino muy bueno y 5 camino excelente 
(Tabla 1) (3). 



Tabla 1 . índice de Servicio Actual (ISA). 



DE 

KILÓMETRO 


A 

KILÓMETRO 


CALIFICACIÓN 
PROMEDIO 


80+000 


85+000 


3,1 


85+000 


90+000 


2,8 


90+000 


95+000 


2,9 


95+000 


100+000 


2,8 


100+000 


106+000 


2,9 



Se realizaron recorridos a lo largo del tramo en estudio 
para observar la superficie de rodamiento con la finalidad 
de apreciar y evaluar la extensión y severidad de los daños 
que presenta la superficie de rodamiento y la estructura de 
pavimento. 



Levantamiento de deterioros 

Los deterioros se clasificaron, para los fines de este 
trabajo, dentro de los siguientes cuatro grupos: textura, 
deformación, agrietamiento y mantenimiento realizado. 
El grado de severidad se midió de acuerdo a la siguiente 
escala: Muy ligera (con calificación de 1); Ligera (con 
calificación de 2); Moderada (con calificación de 3); Severa 
(con calificación de 4) y Muy severa (con calificación de 
5). Igualmente se consideró la intensidad del deterioro 
(porciento de área afectada), la cual se califica como: 
Menor de 10% del área (Esporádica); 10% al 20% del área 
(Aislada); 20% al 50% del área (Frecuente); 50% al 80% del 
área (Muy frecuente) y del 80% al 1 00% del área (Total) (4). 

En la superficie de rodamiento en el tramo estudiado 
se observaron los siguientes deterioros: la textura presenta 
una severidad de 3, las deformaciones se presentaron 
esporádicamentecon severidad entre 1 y 2,el agrietamiento 
con severidad de 3 y mantenimiento (bacheo superficial y 
bacheo profundo) con severidad entre 2 y 3. 

Los últimos trabajos de conservación en los tramos 
en estudio a nivel carpeta asfáltica se realizaron hace 
dos años; a la fecha no se ha realizado ningún trabajo de 
reconstrucción o colocado un refuerzo en la estructura de 
pavimento; se han realizado trabajos de riego de sello; el 
flujo vehicular se ha ¡do incrementando, presentándose un 



alto porcentaje de vehículos pesados; el índice de servicio 
se encuentra en la etapa media al nivel de rechazo y esto 
se ha logrado mantener por arriba de este nivel por los 
trabajos de conservación rutinaria que se han realizado. 
La carpeta asfáltica presenta agrietamientos moderados 
por la oxidación y/o endurecimiento del cemento asfáltico 
lo cual es causado por su vida en servicio. Esto ocasiona 
que se vuelva permeable, originando deformaciones, 
desprendimiento y agrietamientos piel de cocodrilo. 

Medición de la capacidad 
estructural del pavimento 

Para conocer el comportamiento elástico de la estructura 
del pavimento bajo la acción de las cargas se efectuaron 
mediciones de rebote, utilizando una viga Benkelman y un 
camión de volteo con peso en el eje trasero de 8.2 Ton y 
presión de inflado de 5.8 kg/cm 2 . 

Con los datos obtenidos se determinó la deflexión 
percentil 80 por el método de California y la deflexión 
característica por el método del Instituto Norteamericano 
del Asfalto, para medir el espesor de refuerzo que necesita 
la estructura de pavimento actual (Tabla 2) (5). 



Tabla 2. Levantamiento de deflexiones. 



SUBTRAMO 
KM A KM 


PERCENTIL 80 
MÉTODO DE 
CALIFORNIA 
plg x 10-3 


MÉTODO DEL 
INSTITUTO 

DEL ASFALTO 
plg x 10-3 


80+340 


80+820 


30,6 


34,8 


83+940 


84+420 


44,5 


53,3 


84+280 


84+760 


39,3 


44,9 


87+940 


88+420 


27,5 


31,6 


89+760 


90+240 


38,0 


41,7 


90+080 


90+560 


31,6 


34,8 


93+780 


94+260 


32,9 


36,2 


95+080 


95+560 


42,8 


49,8 


97+480 


97+960 


32,5 


37,1 


99+040 


99+520 


39,9 


46,1 


100+000 


100+480 


43,1 


50,2 


102+140 


102+620 


43,3 


52,0 


105+360 


105+840 


25,5 


32,4 



Estudio geotécnico y pruebas de laboratorio 

Con el objeto de conocer el tipoy la calidad de los materiales 
que actualmente forman la estructura de pavimento, así 
como los espesores de las diferentes capas que lo forman, 
se realizaron sondeos del tipo pozo a cielo abierto, se 
determinó el peso volumétrico del lugar, el contenido de 
humedad natural; se tomaron muestras representativas 
para enviarlas al laboratorio para realizar las pruebas de 
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Tabla 3. Espesores existentes de las capas de pavimento. 



UBICACIÓN 


ESPESOR 


DEL SONDEO 


CARPETA 


BASE 


SUBRASANTE 


Km. 


(cm) 


(cm) 


(cm) 


80+520 


16,50 


14,00 


30,00 


81+040 


19,40 


20,30 


30,00 


82+100 


13,80 


14,30 


30,00 


82+640 


17,40 


11,50 


30,00 


84+200 


13,20 


16,50 


30,00 


85+180 


16,10 


28,30 


30,00 


86+540 


13,00 


11,00 


30,00 


87+520 


16,10 


17,40 


30,00 


88+040 


15,80 


12,40 


30,00 


89+500 


11,20 


11,90 


30,00 


90+580 


11,50 


15,30 


30,00 


91+120 


14,30 


9,40 


30,00 


92+300 


13,80 


13,80 


30,00 


93+260 


13,90 


13,20 


30,00 


94+120 


13,80 


14,80 


30,00 



calidad correspondientes (Tabla 3). 

Para localizar las zonas homogéneas, se realizaron las 
siguientes actividades: obtención del ISA, levantamiento 
de daños, medidas de deflexión y ejecución de los sondeos 
de exploración en los que se obtuvo información de la 
estructura de pavimento. 

Con el objeto de conocer el tipo y la calidad de los 
materiales que actualmente forman la estructura de 
pavimento, así como los espesores de las diferentes capas 
que lo forman, se realizaron sondeos del tipo pozo a cielo 
abierto, se determinó el peso volumétrico del lugar, el 
contenido de humedad natural; se tomaron muestras 
representativas para enviarlas al laboratorio para realizar 
las pruebas de calidad correspondientes. 

De los análisis de laboratorio se hacen los siguientes 
comentarios: 

Capa de terracería: el material corresponde a una arena 
arcillosa (SC), con buenas características mecánicas para su 
uso, y el grado de compactación varía de 77.2% a 95.4%. 

Cuerpo de terraplén: está formado por los materiales 
de préstamos laterales; el material corresponde a una 
arena arcillosa (SC), con buenas características mecánicas 
para su uso, y el grado de compactación varía de 79.9% a 
99.0%. 

Capa subrasante: corresponde a una capa de 30 cm; se 
construyó de manera similar a la de terracerías, ya que el 
material es adecuado también para utilizarse en esta capa, 
y el grado de compactación varía de 84.3% a 1 04.9%. 

Capa de base hidráulica: corresponde a un material 
con características adecuadas para formar esta capa, es 
material procedente de banco de grava arena de arroyo; 



UBICACIÓN 


ESPESOR 


DEL SONDEO 


CARPETA 


BASE 


SUBRASANTE 


Km. 


(cm) 


(cm) 


(cm) 


95+100 


13,80 


13,80 


30,00 


96+400 


13,20 


13,20 


30,00 


97+860 


16,50 


16,50 


30,00 


98+520 


13,20 


13,20 


30,00 


99+380 


10,40 


10,40 


30,00 


100+300 


14,50 


14,50 


30,00 


101+240 


16,50 


16,50 


30,00 


102+620 


13,40 


13,40 


30,00 


103+600 


11,50 


11,50 


30,00 


104+700 


9,90 


9,90 


30,00 


105+260 


10,60 


10,60 


30,00 


106+980 


16,70 


16,70 


30,00 


107+420 


12,90 


12,90 


30,00 


108+960 


16,20 


16,20 


30,00 











tiene una clasificación de una grava arena limosa (GM- 
SM), los espesores varían de 8.0 cm a 35.0 cm y el grado de 
compactación varía de 86.3% a 107.2%. 

Capa de base asfáltica: el material corresponde a la 
recuperación de la carpeta asfáltica y la capa de base 
hidráulica mezclada con emulsión asfáltica; los espesores 
varían de 9.0 cm a 25.0 cm y el grado de compactación 
varía de 93.8% a 103.2%. 

Carpeta asfáltica: presenta espesores que varían de 
3.4 cm a 33.3 cm, se encuentra oxidada, agrietada y está 
llegando al final de su vida útil. 

Estudios del tránsito 

Para determinar el análisis del tránsito se utilizaron los 
datos viales que edita la Secretaría de Comunicaciones y 
Transportes, del año 2001 a 2012, y el número de clave 
de la autopista Estación Don - Nogales, que corresponde 
al 26445. La estación está ubicada en el km 128+000. El 
TPDA es de 6,652 vehículos/día, la composición vehicular 
es de vehículos tipo A = 73.6%, tipo B = 5.1%; tipo C = 
21 .3%; la tasa de crecimiento es de 0.43% y el coeficiente 
de distribución del tránsito es de 40% (6). 

Para el diseño de las alternativas se consideraron los 
datos viales, los parámetros de resistencia obtenidos en el 
laboratorio y se aplicaron los métodos de diseño como el 
Método del Instituto de Ingeniería de la UNAM, el Método 
de California, el Método AASHTO, el Método del Instituto 
Norteamericano del Asfalto y el Método de Fatiga de la 
UNAM; lo anterior para determinar los espesores de las 
capas que conformarán la estructura de pavimento. El 
resumen de los resultados se muestra en la tabla 4. 
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Tabla 4. Resumen de espesores de carpeta asfáltica. 





PERÍODO 
DE DISEÑO 
(AÑOS) 


ESPESOR (cm) 


ESPESOR 
RECOMENDABLE (cm) 


TRAMO DE 
ANÁLISIS 


MÉTODO 
DEL INSTO. 
ING.DELA 
UNAM 


MÉTODO 
AASHTO 


MÉTODO 
DE 

CALIFORNIA 


MÉTODO DEL 
INSTO. NORT. 
DEL ASFALTO 


MÉTODO 
DE FATIGA 
DELA 
UNAM 


CARPETA 


BASE DEL 
ASFALTO 


km 80+000-106+000 


15 


3,7 


5,6 


23,4 


10,4 


5,0 


9,6 


23,4 



DISEÑO DE ALTERNATIVAS DE SOLUCIÓN 
PARA LA REHABILITACIÓN 

Las tres alternativas que se consideran viables para la 
solución de las estructuras de pavimento analizadas 
apuntan hacia un mejoramiento de la base hidráulica que 
puede ser con cemento portland o emulsión asfáltica; la 
otra alternativa es mejorar la superficie de rodamiento 
utilizando carpeta de renivelación para el sellado de las 
grietas y posteriormente colocar una carpeta asfáltica. 

Primer alternativa de rehabilitación 

Cortar la carpeta asfáltica y la base hidráulica existente 
mediante disgregado con la máquina recuperadora 
RR-350 a la profundidad de 25.0 cm para formar una 
base estabilizada con cemento portland al 9.0% PVSM, 
compactada al 100%, y construir una carpeta de concreto 
asfáltico con asfalto modificado de 5 cm de espesor, 
compactada al 95%. 

Segunda alternativa de rehabilitación 

Cortar la carpeta asfáltica y la base hidráulica existente 
mediante disgregado con la máquina recuperadora RR-350 
a la profundidad de 25.0 cm, adicionarle emulsión asfáltica 
superestable en proporción de 80 L/m 3 para formar una 
base asfáltica de 25.0 cm de espesor, compactada al 1 00%, 
y construir una carpeta con concreto asfáltico con asfalto 
modificado de 6.0 cm de espesor, compactada al 95%. 

Tercer alternativa de rehabilitación 

Se realizará un bacheo de caja y posteriormente se tenderá 
una carpeta de renivelación de 8.0 cm con concreto 
asfáltico compactada al 95%, tender una malla geotextil 
y construir una carpeta de concreto asfáltico con asfalto 
modificado de 1 0.0 cm de espesor, compactada al 95%. 

CONCLUSIONES Y RECOMENDACIONES 

El problema principal de la red carretera nacional 
se encuentra en su conservación, rehabilitación y 



modernización, ya que gran parte de ella se está operando 
con índices de servicio actual bajos, cercanos al nivel de 
rechazo, produciéndose en consecuencia un importante 
incremento en el costo global del transporte. Esta 
situación, se origina por falta de recursos suficientes y 
oportunos, con asignaciones del orden de una cuarta parte 
de las requeridas. Por lo anterior es grave en la red troncal 
básica. 

Es importante que la conservación, rehabilitación 
o modernización de la red nacional, se refiera al manejo 
de los pavimentos, donde se pretenden agrupar los 
métodos de auscultación, los criterios para la selección y 
el tratamiento de los materiales, las normas que se utilizan 
para su caracterización, las soluciones de estructuración, 
los métodos de diseño de espesores, los procedimientos 
constructivos y los criterios o posibilidades de inversión. A 
este respecto ya se tienen avances considerables para la 
creación de un Sistema Mexicano de Administración de 
Pavimentos. 
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MÉTODO ALTERNO PARA DETERMINAR EL COSTO 
DE PRODUCCIÓN EN LÍNEAS DE ENSAMBLE 



JAIME DÍAZ SANTANA 



En este artículo se presento una metodología para determinar el costo de la 
producción en líneas de ensamble mediante el sistema de costeo por procesos. Este 
método es aplicable en empresas que trabajan bajo pedidos específicos de productos 
cuyas corridas de producción, por pedido, abarcan varios meses. Lo anterior es posible 
mediante la transformación virtual del sistema de producción real (por pedidos u 
órdenes de trabajo) a un sistema de producción por procesos. 
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INTRODUCCION 

El costeo de la producción en líneas de ensamble se apega 
más al método de costeo por órdenes de trabajo; sin 
embargo, también puede efectuarse como una variación 
del costeo por procesos. Las industrias más representativas 
para aplicar esta variación de costeo utilizan mano de 
obra directa en forma intensiva en la fabricación de 
sus productos, como es el caso de las maquiladoras de 
manufactura. En estas empresas la producción se realiza en 
líneas de producción llamadas líneas de ensamble (1). Las 
líneas de ensamble se dividen en estaciones de trabajo; en 
cada estación labora un obrero; cada obrero efectúa una o 
varias operaciones y, de acuerdo a la demanda a satisfacer, 
se balancea la línea (labor que corresponde a la ingeniería 
industrial), es decir, se determinan cuántas estaciones 
de trabajo ¡guales (en cuanto al contenido de trabajo) se 
requieren y cuáles operaciones se deben efectuar en cada 
una de ellas (2). Para el caso que nos ocupa, cada conjunto 
de estaciones con igual contenido de trabajo formarían el 
equivalente de un proceso. 

Mediante el método de costeo de la producción por 
procesos, los costos correspondientes se calculan cada 
mes en cada proceso (3). La producción terminada de cada 
uno de ellos se transfiere al siguiente proceso, el cual la 
recibe con su correspondiente costo y le agrega más valor. 
La producción terminada en el último proceso refleja el 
costo total de los productos que estarán disponibles para 
su venta o se envían directamente a los clientes de acuerdo 
a sus pedidos. 

Cada mes, en cada proceso, es necesario calcular los 
grados de avance de los diferentes elementos del costo de 
producción (materia prima directa, mano de obra directa 
y gastos indirectos de fabricación) de las unidades que no 
se terminaron de fabricar completamente, es decir, del 
inventario final de producción en proceso. Lo anterior es 
para determinar el número de unidades equivalentes y los 
costos de producción. 

En cada proceso, el costo unitario de producción se 
calcula como un promedio del costo total de producción 
del proceso entre las unidades de producto totalmente 
terminadas y éstas se calculan mediante la suma de las 
unidades que ya tienen integrada la materia prima, la 
mano de obra y los gastos indirectos de fabricación que 
le corresponden ser agregadas en un proceso, todo 
al 100%, y la producción que quedó sin terminar, la 
cual debe convertirse a su equivalente de producción 
terminada, es decir, debe determinarse a cuántas unidades 
100% terminadas equivale esta producción que quedó 
inconclusa. 

En cuanto a la producción inconclusa, consideremos 
el siguiente ejemplo: si al final de un mes quedaron 4 
unidades de producto sin terminar a las cuales sólo se les 
agregó, en promedio, el 50% de la materia prima del total 
requerido, éstas equivalen a 2 unidades 100% terminadas, 
es decir, a dos unidades equivalentes de producto 
terminado en cuanto a materia prima se refiere. 



CASO ESPECIFICO 

Se utilizará un ejemplo muy simplificado puesto que 
el objetivo de este artículo es el de ¡lustrar cómo puede 
determinarse el costo de producción en líneas de ensamble 
mediante el sistema de costeo por procesos con base en 
costos promedio. 

Supondremos que para fabricar una camisa 
primeramente se efectúa el corte de la tela y luego la 
formación de las mangas. El formado de cada manga 
requiere de tres operaciones: formado del puño (OP1), 
cosido de la parte superior de la manga (OP2), y cosido 
de la parte superior de la manga con el puño (OP3); no 
se incluirán las materias primas hilo y botones. Es lógico 
que siguen otras operaciones para completar el producto: 
la camisa; sin embargo, nos referiremos solamente a la 
formación de la manga, sea la izquierda o derecha, para 
¡lustrar el caso que nos ocupa. Finalmente diremos que 
para satisfacer la demanda de camisas requerida por los 
clientes, el balanceo de la línea de producción requiere 
de dos estaciones de trabajo para corte de las piezas de la 
camisa y de cinco estaciones de trabajo para la formación 
de mangas. Las dos estaciones de corte formarían el 
proceso de corte y las cinco estaciones de formado de 
mangas formarían el proceso de formado de mangas; y en 
éste, en cada estación de trabajo hay un obrero y cada uno 
de ellos realiza las mismas tres operaciones de formado de 
la manga (Figura 1). 
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Figura 1. Distribución esquemática de las estaciones 
de trabajo para formar la manga de la camisa. 

Es lógico suponer que del proceso de corte se envían 
las otras piezas de las camisas a otros procesos para su 
ensamble, de acuerdo a como se haya determinado en el 
balanceo de la línea de producción. Sin embargo, se reitera 
que nos referiremos solamente al proceso de formado de 
mangas, donde sólo se agrega costo por concepto de la 
mano de obra directa empleada y por los gastos indirectos 
de fabricación incurridos. 



EPISTEMUS 



Jaime Díaz Santana: UNISON /EPISTEMUS 14/ Año 7/2013/pág.: 59-63 




En cada estación de trabajo del proceso del formado 
de mangas se efectúan las tres operaciones, anteriormente 
explicadas, sobre las materias primas MP1 Y MP2 (telas para 
formar el puño y la manga de la camisa, respectivamente), 
las cuales provienen del proceso de corte, y las realiza un 
solo obrero. En la tabla 1 se presentan los tiempos de 
operación y los porcentajes de cada tiempo de operación 
con relación al tiempo total de mano de obra para formar 
la manga. 

Tabla 1 . Tiempo de las operaciones para formar la manga. 



Operación 


Tiempo (min) 


Tiempo (%) 


OP1 


0.3 


25 


OP2 


0.6 


50 


OP3 


0.3 


25 


Total 


1.2 


100 



Supondremos que en un mes determinado se ha 
generado la siguiente información en el proceso de 
formado de mangas. Una unidad de producto en este 
proceso corresponde a una manga formada. 

El Inventario Inicial de Producción en Proceso (IIPP) 
al inicio del mes fue de 55 unidades inconclusas con los 
siguientes costos de las materias primas, que incluyen 
los costos de mano de obre directa y gastos indirectos 
de fabricación; y que fueron transferidas por el proceso 
de corte al proceso de formado de mangas: MP1 = $440, 
MP2 = $785. El valor agregado por este último proceso 
a las anteriores materias primas del inventario inicial fue 
de $680 por concepto de mano de obra directa y $1 30 de 
gastos indirectos de fabricación. 

Durante el mes se agregaron, procedentes del 
departamento de corte, materiales para producir 13,000 
unidades de producto (mangas), con un costo de $1 04,000 
para MP1 y $195,000 para MP2. La mano de obra utilizada 
en el mes y los gastos indirectos de fabricación fueron 
$160,000 y $30,000, respectivamente. 

El inventario de producción en proceso al final del mes 
(IFPP) fue de una unidad incompleta en cada estación de 
trabajo, con los materiales y mano de obra utilizados, y 
grados de avance de esta última para formar una manga, 
tal y como se muestra en las tablas 2 y 3. 



Tabla 2. Materiales utilizados en el inventario final de 
producción en proceso. 



Estación de trabajo 


Materiales utilizados 


del proceso 


MP1 MP2 


1 Si Si 


2 


Si No 


3 


Si No 


4 Si Si 


5 Si Si 



Total 



Tabla 3. Operaciones realizadas sobre el inventario final de 
producción en proceso y grados de avance de la mano de 
obra en cada estación de trabajo. 



Estación de 

trabajo 
del proceso 



Operaciones 
realizadas 



Tiempo de Avance de 
operaciones mano de obra 
(min) por unidad 

















Min 


% 


1 


Si 


Si 


No 


03 


06 


00 


0.9 


75 


2 


Si 


No 


No 


03 


00 


00 


0.3 


25 


3 


Si 


No 


No 


03 


00 


00 


0.3 


25 


4 


Si 


Si 


No 


03 


06 


00 


0.9 


75 


5 


Si 


Si 


No 


03 


06 


00 


0.9 


75 


















275% 



En el inventario final también quedaron sin utilizar, es 
decir, sin trabajo alguno, 63 piezas de MP1 y 65 de MP2 
(Tabla 4), las cuales, sumadas a los materiales que tienen 
algún grado de proceso: cinco de MP1 y tres de MP2 (Tabla 
2), dan un total de 68 piezas de cada uno. Dicho en otras 
palabras, quedaron en el inventario final, con o sin trabajo 
agregado, materiales para el formado de 68 mangas. 

Con la información proporcionada se procede a aplicar 
el procedimiento para costear la producción del proceso 
de formado de mangas, mediante el procedimiento de 
costo promedio (4). 



Tabla 4. Número de piezas de MP1 y MP2 que quedaron sin 
trabajo alguno. 







Estación de trabajo 




Material 






Total 


MP1 


10 


15 12 18 8 


63 


MP2 


10 


16 13 18 8 


65 
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PASO 1 : Flujo de unidades físicas 



En este paso se determina el número de unidades que 
en el período se transfieren al siguiente proceso o que se 
terminan y estarán disponibles para su venta. Para el flujo 
de las unidades físicas que intervinieron en el proceso en 
el mes de cálculo no se consideran los grados de avance de 
las unidades. El procedimiento se presenta en la tabla 5, en 
unidades de producto en todos los renglones. 



Tabla 5. Determinación de la producción transferida al 
siguiente proceso en unidades de producto. 



Inventario inicial de producción en proceso 


55 


+ Agregadas en el período 


1 3,000 


Disponibles en producción 


13,055 


- Inventario final de producción en proceso 


68 


Producción transferida al siguiente proceso 


1 2,987 



PASO 2: Unidades equivalentes 



Todos los materiales que quedaron, tanto en el IIPP 
como en el IFPP, se consideran con el 100% de avance en 
cuanto a sus costos de MPD, MOD y GIF que se agregaron 
en el proceso de corte. 

En el proceso de formado de mangas no se agregan 
nuevas materias primas. A las materias MP1 y MP2 recibidas 
del proceso de corte se les incrementa el valor agregado 
con el costo de la mano de obra directa y los gastos 
indirectos de fabricación que se generaron en el mes. 

Quedaron en proceso, con o sin trabajo agregado, 68 
unidades de producto. En cuanto a la mano de obra directa, 
la suma de grados de avance que se muestra en la tabla 3, 
de las unidades que tuvieron algún trabajo es del 275%. 
Las 68 unidades inconclusas, en lo relativo a mano de obra 
directa (MOD), tienen un avance promedio de 275%/ 68 
ud = 4.0441 18% y las unidades equivalentes en MOD son 
igual a 68 ud (0.040441 1 8) = 2.75 unidades. 

Para este ejemplo diremos que la aplicación de los 
gastos indirectos de fabricación al costo de los productos 
se realiza con base en la mano de obra empleada en 
su fabricación. Cuando lo anterior ocurre, el por ciento 
de avance de la producción inconclusa y las unidades 
equivalentes, en cuanto a la mano de obra se refiere, son 
idénticos en lo relativo a los GIF. 

La producción que este proceso transfiere al siguiente 
se encuentra al 100% de avance en cuanto a la materia 



prima y trabajo que tenía que agregarse en el proceso de 
corte y en el de formado de mangas. Los resultados de este 
paso se presentan en la tabla 6. 



Tabla 6. Unidades equivalentes de la producción transferida 
y del inventario final de producción en proceso. 



MATERIA PRIMA 
MP1 1 MP2 


MANO DE 

OBRA 
DIRECTA 


GASTOS 
INDIRECTOS 
DE 

FABRICACIÓN 


Producción transferida 
100% de avance 


12,987 


12,987 


12,987 


12,987 


Más IFPP (68 unidades) 










MP1 

68 unidades con 
100% de avance 


68 








MP 

68 unidades con 
100% de avance 




68 






Mano de obra directa 

68 unidades con 
4.04411 11 8% de avance 






2.75 




Gastos indirectos de 
fabricación 

68 unidades con 
4.04411 11 8% de avance 








2.75 


Total unidades equivalente 


s 13,055 


13,055 


12,989.75 


12,989.75 



En el costeo de la producción con base en los costos 
promedio se considera que en cada proceso el costo del 
inventario inicial de producción en proceso, más el costo 
que se agrega en un período, se convierte en el costo del 
inventario final de producción en proceso más el costo 
de la producción que se termina o que se transfiere al 
siguiente proceso. Estos dos últimos elementos se calculan 
con base en los costos promedio de cada uno de los 
elementos del costo de producción: materia prima directa, 
mano de obra directa y gastos indirectos de fabricación. 
Los costos promedio se obtienen mediante una división de 
los dos primeros costos (del inventario final + el agregado) 
entre las unidades equivalentes correspondientes que se 
trabajaron en el período. Los cálculos correspondientes se 
presentan en las tablas 7, 8 y 9. 

PASO 3. Costos unitarios (Tabla 7) 



Tabla 7. Determinación de los costos unitarios. 





IIPP 

($) 


Agregado 

($) 


Total 

($) 


Unidades 
equivalentes 


Costo 
unitario 

($) 


MP1 


440.00 


104,000.00 


104,440.00 


13,055.00 


8.00 


MP2 


785.00 


195,000.00 


195,785.00 


13,055.00 


15.00 


MOD 


680.00 


160,000.00 


160,680.00 


12,989.75 


12.37 


Total 


2,035.00 


489,000.00 


30,130.00 


12,989.75 


2.32 
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PASO 4. Costo del inventario final de producción en 
proceso (Tabla 8) 



Tabla 8. Determinación del costo del inventario final de 
producción en proceso. 





Unidades 
equivalentes 


Costo 
unitario ($) 


Total ($) 


MP1 


68.00 


8.00 


544.00 


MP2 


68.00 


15.00 


1,019.79 


MOD 


2.75 


12.37 


34.02 


GIF 


2.75 


2.32 


6.38 


Total 






1,604.19 



PASO 5. Costo de la producción transferida al 
siguiente proceso (Tabla 9) 



El anterior resultado puede obtenerse también como 
se muestra en la tabla 1 



Tabla 10. Cálculo alterno para determinar el costo de la 
producción transferida. 



Costo delllPP 


$ 2,035.00 


+ Costo agregado 


$ 489,000.00 


Costo disponible 


$491,035.00 


-Costo del IFPP 


$ 1,604.19 


Costo producción transferida 


$489,430.81 



Si bien el ejemplo resuelto se refiere a un solo proceso, 
la metodología es la misma para los otros procesos que 
resultarán y, una vez resuelto el último, se tendría el costo 
unitario promedio de cada camisa producida en el período. 



Tabla 9. Determinación del costo de la producción 
transferida. 





Unidades 
equivalentes 


Costo 
unitario ($) 


Total ($) 


MP1 


12,987.00 


8.00 


103,896.00 


MP2 


12,987.00 


15.00 


194,765.21 


MOD 


12,987.00 


12.37 


160,645.98 


GIF 


12,987.00 


2.32 


30,123.62 
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INTRODUCCION Y MOTIVACION 

Por ser de naturaleza local, todas las definiciones que 
daremos en este artículo son válidas en cualquier 
subconjunto abierto A de R n ; sin embargo, generalmente 
escribiremos IR n en vez de A por simplicidad de notación. 

Empezaremos por motivar el concepto de diferencial 
de una función. 

Recordemos que por definición una función /: R— ► IR 
es diferenciable en a e R si existe un número f(a) tal que 

f(a+ h) - f(a) 



lim 



h 



f'(a). 



Si definimos la función lineal A:R— ^R 
A (h) = f(aj • h entonces la igualdad (1 ) equivale a 

f(a+h)-f(a)-X(h) 



lim 



h 



0. 



(1) 



por 



(2) 



La interpretación de la igualdad (2) es que la función 
A+/(a)es la mejor aproximación lineal de /en el punto a. 
La expresión (2) nos da la ¡dea para generalizar la definición 
de diferenciabilidad para una función / : R n — ^ R m : 

Una función / : R n ^R m es diferenciable en p e R n si 
existe una función lineal A : R n — ^ R™ tal que 



lim 



l/(p+ h)-f(p)-X(h)\ 



0. 



Notemos que como h es un punto en R n y 
f{p+h)-f{p)-X(h) es un punto en R™ , es indispensable 
considerar las normas. Se demuestra en [7] (Teorema 
2-1, pág. 14) que de existir la función lineal A, es única. 
Denotaremos a A por Df(p) y llamaremos a la función 
lineal Df(p) :R n ^R m la diferencial de la función / en 
el punto p. 

También se demuestra en [4] (Teorema 1 6, pag. 1 25) que 



(Dfíip) 



í & (P) 

OX\ 

\ t i(p) 



f ' (p) \ 

OXn 



dx r 



/ 



La matriz anterior se llama matriz Jacobiana. Ver [7], 
pag. 15. 

Para obtener el valor de (D f) (p) en un vector v en R n 
usando la matriz Jacobiana, simplemente consideramos el 
vector v como vector columna y lo multiplicamos por la 
izquierda por dicha matriz. 



Diremos que una función / : R n — ^ R m es diferenciable 
en R n si / es diferenciable en cada punto p en R n . 

El caso que nos interesa en este artículo es cuando m=1 . 
En este caso tenemos que la diferencial Df(p) : R n — » R de 
la función /: R n — R está dada por 



(3) 



Esto es, Df(p) = V p f donde por definición 
v p = (-^-(p),..., ^- (p)). Para cada v en R n tenemos que 
Df(p) (v) es un escalar que representa la rapidez de cambio 
de la función / en el punto p en la dirección del vector v. 

Ahora procederemos a motivar los conceptos de 
1 -forma y producto cuña entre 1 -formas. 

Sea/: R n — Runa función diferenciable. Si {e l5 . . . e n } 
es la base canónica de R n tenemos que cualquier vector 

v en R n se expresa como v= v 1 e 1 -\ h v n e n con v J 

escalares en R. Entonces 

n 

Df(p)(v)=Df(p)(Y]v>e J ) 

3=1 

n 



3=1 



= E 



donde la segunda igualdad es por la linealidad de 
Df(p)y\a última igualdad es por (3). 



Obtenemos así la igualdad 



Df(p)(v) = J] 



3=1 



3 9 f 

r : 

dx 3 



(P). 



(4) 



Dado que también las proyecciones canónicas 

x^:R n ^R 

definidas por x % (p 1 , . . .p n )= p l son funciones 
diferenciables tenemos que 



dx 1 



i=i 



i=i 
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Jules Henri Poincaré 

donde la primera igualdad es por (4) y donde á¿¿ es la 
función delta de kronecker definida por 



1 si i = j 
si i ^ j 

Obtenemos de esta manera la igualdad 

Dx\p)(v) = v\ (5) 
Por las igualdades (4) y (5) tenemos que 

df 



Df(p)(v)=Y]Dxi(p)( 



3=1 



Reordenando los términos en la igualdad anterior 
y denotando a Dx j por dx j obtenemos la conocida 
expresión de cálculo 



3=1 



ÉL 



(6) 



Observemos que el conjunto ordenado 
{dx 1 (p), . . . dx n (p)} es la base dual para el espacio 
vectorial dual (R n )* asociada a la base ordenada canónica 
{e l5 . . . , e n } de (R n pues por construcción se cumple que 
dx l (p){e 3 ) = Síj para cada p en R n . (Consultar [2] pág. 98, 
para la definición de base dual). 



Motivados por la expresión (6) asociada a la función 
diferenciable f, una pregunta natural que surge es la 
siguiente: 

Dada cualquier "1 -forma" w, es decir una expresión de 
la forma 

w = gidx 1 + • • • + g n dx n 

donde las 9i son funciones de R n en IR diferenciabas, 
¿existe una función /: R n R con segundas derivadas 
parciales continuas tal que Df = g 1 dx 1 +---+g n dx n 7 
Esta es precisamente la pregunta inicial a la que se hace 
referencia en el resumen de este artículo. 

De existir dicha / debe cumplir que 



dx j 



X! 9jdx j 

dado que los dx J son base la igualdad anterior se da 
si y sólo si 

df 



dxi 



■ Qj para cada j = 1, . . . ,n 



(7) 



Por otra parte, la condición de que / tiene segundas 
derivadas parciales continuas implica que las segundas 
derivadas parciales mixtas de / deben coincidir (Ver [4] 
Teorema 1 5, pág. 1 1 9 o [7] Teorema 2-5, pág.24), esto es, 
para cadal < i,j < n se cumple 

d _(®L\ 9 ( d f\ 
üi\dx l J dx l \dxi/' 



dx? 



Entonces por (7) y (8) tenemos que 



dgi _ dgj 



es necesario que se cumpla la condición 

% dgj 
dx j dx l 
para que exista la / requerida. 



, es decir, 



paral < h3 <™ (9) 



Reescribiremos la condición (9) introduciendo nueva 
simbología. Si definimos el "producto" A en el conjunto 

{dx 1 , . . . dx n } mediante la regla dx l Adxi = — dx J Adx 1 
podemos expresar la condición (9) como la única ecuación 



n n 
3=1 i=l 



O- 



(10) 



Así si w = Yjí=i gjdx j es una 1 -forma, definiendo el 
operador d por 



dw ^^Dgj Adx j 



(11) 



obtenemos que la condición (10) se traduce en lo 
siguiente 
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n n n r>. 

dw = Dg¡ A dx j = ( ¿ dx¿ ) A = 

¿=1 ¿=1 i = l 

donde la última igualdad de la expresión anterior es 
por (10). 

Así, con este nuevo lenguaje, obtenemos que una 
condición necesaria para que la 1-forma w = Yj ? i=i 9jdx j 
tenga una "primitiva" / ( esto es, Df = w ) es que dw = 0. 

Veremos mas adelante que la condición dw = no es 
una condición suficiente y refinaremos los conceptos de 
producto cuña A y operador d, pero antes hagamos un 
pequeño resumen de lo que se hizo en esta sección: 

Vimos que el concepto de operador deferencial D 
para funciones / de R n en IR nos conduce de manera 
natural a definir el concepto de 1-forma. Una vez 
definido este concepto surgió la pregunta natural de 
que si damos cualquier 1-forma w, ¿la 1-forma tiene 
siempre una primitiva /? Para tratar de contestar esta 
pregunta fue conveniente introducir un "producto cuña 
A" en el conjunto de 1-formas {dx 1 , . . . dx n } tal que 
dx % Adxj = —dxi Adxi. Lo que a su vez nos condujo 
a introducir el concepto de "derivada de 1-formas"o 
operador diferencial de 1-formas. Así si w = ^gidx 1 
es una 1-forma, definimos dw = y] Dqí Adx 1 , es decir, 

dw = SEÍ^;^) Adx* y resultó, con este nuevo lenguaje, 
que una condición necesaria para que una 1-forma dada 
tenga una primitiva es que dw = 0. 

CAMPOS VECTORIALES 

Antes de definir formalmente el concepto de 1-forma 
presentaremos en esta sección el concepto de campo 
vectorial y algunas de sus propiedades. Las definiciones de 
campo vectorial y de 1-forma son totalmente análogas y 
comparten algunas propiedades en común desde el punto 
de vista algebraico. La razón por la que presentamos 
primero el concepto de campo vectorial es porque su 
definición es un poco mas simple que la definición de 
1 -forma y su interpretación geométrica es mas evidente. 

Definición 1 . Sea p G R n . Definimos el conjunto 

R n p = {( P ,v)\ v e R n } 

Llomoremos ol conjunto R™ espacio tangente o R n en p. 

Para cada p G R n , R™ es un espacio vectorial sobre R 
con operaciones definidas por 

(p,v 1 )+ (p,v 2 ) = (p,v 1 +v 2 ) 

a(p,v) = (p,av) 

Vi : v 2 G R n , a G R. Denotaremos a (p,v) por v p de 
aquí en adelante. 



Definición 2. Un campo vectorial es una función 

talqueF(p) G para cada p G R n . 

La condición F(p) G Rp para cada p G R n de la 
defnición 2 es equivalente a decir que se cumple que 
7ro F{p) = p para cada p G R n donde 

71 : U K ~> Rn 

peR n 

es la proyección canónica definida por yr (y J = p para 
cada p G R n . 

Dado que por definición un campo vectorial 
F es tal que F(p) G ' entonces tenemos que 
F (p) = F 1 (p) ( ei ), + ••• + F n {p)(e n ) p donde F 1 (p),... F n (p) 
son escalares en R que dependen del punto p. Es decir, 
tenemos que las F l son funciones de R n en R. Este hecho 
nos permite formular la siguiente definición: 

Definición 3. Decimos que un campo vectorial F 
es continuo o diferenciable si cada una de las funciones 
F % : R n —> R lo son en el sentido usual. 

La definición 2 y el hecho de que para cada p G R n 
el conjunto R^ tiene estructura de espacio vectorial nos 
permite inducir (punto por punto) operaciones entre 
campos vectoriales: 

Definición 4. Sean Fy G campos vectoriales. Definimos las 
operaciones 

(F + G){p) = F{p)+G{p) 

<F,G>(p)=<F(p) J G(p)> 
(fF)(p) = f(p)F(p) 

donde p está en R n ? f : R n —> R es diferenciable y < , > 
denota el producto punto en R n . 

Observemos que F + G y fF son de nuevo campos 
vectoriales y que < F f G > es una función de R n en R . 
Mas aún, si F y G son campos vectoriales diferenciabas 
(continuos) y f es una función diferenciable (continua) 
entonces no es difícil demostrar que los campos vectoriales 
F+G y fG son diferenciables (continuos) y que la función 
< F f G > es diferenciable (continua). 

ÁLGEBRA EXTERIOR 

En esta sección escencialmente lo que haremos es 
definir formalmente los conceptos de k-forma y operador 
cuña A entre k-formas, pero en un sentido estrictamente 
algebraico. Dejaremos la definición geométrica formal de 
k-forma para la siguiente sección. 
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Definición 5. Sea V un espacio vectorial de dimensión n 
sobre R. Seak > 1 un entero. Definimos A k (V*) como el 
conjunto 

A k (V*) = {a :Vx • • -x V— »IR \a es multilineal y alternante} 

k - veces 

donde alternante significa que a(v a ^ ,...v a (k)) = 
Signo (o)ol(v\ ,...,v k ) para toda permutación a en el grupo 
simétrico Sk, Vi G V .( Ver [1] para el concepto de grupo de 
permutaciones). 

Para k=0 definimos A k (V*) simplemente como el 
conjunto de números reales R. 

Esimportanteseñalarque para cualquier k > se tiene 
que A k (V*) es un espacio vectorial sobre R.( Consultar [7], 
pag 72). 

Definición 6. Sean (fi,...,(fk elementos en V* donde V* 
es el espacio dual de V. Definimos el producto exterior 

(Pi A" - A(p k : V x- • x V — ► IR por 

k - veces 

¥>i A- • • Acp k [v ly m m v k ) = detiiflvj))^^ 

donde Vi£ V,((pi(vj))l<ij<k es la matriz cuadrada 
de tamañno kxk con entradas (fi(vj) y det significa 
determinante. 

El hecho de que el determinante es una función 
multilineal y alternante implica que ip 1 A' — A(fk es 
un elemento de A k (V*). ([7], pag.75). Esto nos permite 
"generalizar" y al mismo tiempo formalizar el producto A 
que discutimos en la primera sección: tomando V = 
en la definición anterior tenemos que A 1 (V*) = ( R£ )* 
y así los elementos dx 1 (p)^ . . . ,dx n (p) pertenecen a 
A 1 ^*). Aplicando la definición 6 y haciendo un simple 
cálculo obtenemos que 

dx\e r )(p)/\dx\p)(e s ) = - dx\p)(e r ) A dx\p)(e s ) 

donde e r y e s son elementos cualquiera de la base 
canónica de R n , lo cual significa que 

dx\p)Adx j (p)= -dx j (p)Adx\p) 

Justo como necesitábamos. Esto es, el elemento 
dx\p)Adxi(p) ya no es sólo un "símbolo", sino que 
es un elemento bien defiido del espacio vectorial A 2 (V*) 
con V = R™ . Retomaremos estos hechos en la siguiente 
sección. 

A continuación analizaremos a A k (V*). Dado que A^ 
(y*) es un espacio vectorial sobre IR es natural preguntarse 
cual es su dimensión y exhibir una base. Este espacio es 



de hecho de dimensión finita y una base se construye a 
partir de los productos cuña de los elementos de la base 
dual asociada a cualquier base de V, mas precisamente 
(Consultar [7], Teorema 4-5, pág. 75 para ver la 
demostración de la siguiente construcción), 

Si {v l5 . . . ,^ n }es una base de Vy {vi , . . . , }es 
la base dual asociada a la base {v±, . . . ,v n } para V*. 
Entonces 

{<A---A?4 |l<¿i<---<¿ife<n} 

es una base de A k (V*). Así dim A h (V*)=( n k ) y 
además 

A k (V*)={ Y, a H...i k v*ii ^■■■Av* k \a il ,_ Ak ER] 

ii<---<i k 

donde la expresión ii< - - - <i k que aparece debajo 
de la sumatoria anterior significa que los sumandos de 
dicha sumatoria están indexados por los elementos de 
la base {v*i A- - - Av* k \ 1 <i\< ■ ■ • <¿fe<n} de 

a\v*). 

La manera natural de extender el producto exterior A 
es la siguiente 

Definición 7. Se define A : A k (V*) X A Z <y*) A h + l 
(V*)por 

aA(3(v 1 ,...,v k ,v k + 1 ,...,v k +i) 

= (¡k+í)T ^ ae S k+l sgna a{v a{1): . . . ,v a[k) )f3(v a{k+1) , . . . ,v a{k+l) ) 

para a£A k (V*), (3 G A l (V*) , v u ...,v k+1 en V. 



FORMAS 

En la sección 1 vimos que la condición dw=0 es 
necesaria para que una 1 -forma W tenga una "primitiva" 
/ y mencionamos que dicha condición no es sufiiente. 
En realidad si lo es cuando el dominio de definición de la 
1 -forma es todo R n . Más precisamente, veremos al final de 
esta sección que el hecho de que una A;-forma w tal que 
dw=0 tenga una primitiva depende de una condición 
topológica del dominio abierto de definición de la A:-forma 
W. Dicha condición es que el dominio abierto en cuestión 
sea contraible (ver [6], pág. 374, para la definición de 
espacio contraible). Este resultado se conoce como Lema 
de Poincaré. 
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En analogía con la definición de campo vectorial dada 
en la definición 2 tenemos la siguiente: 

Definición 8. Seo k>0 un entero. Uno k-formo es uno 
función 

w:R n^ y a*((r;)*) 

peR n 

tolque w(p)eA k {{Rp)*) poro codo peR n . 

En el caso de 1-formas tenemos que A 1 ((R^)*) = (Rp* . 
Nótese la clara analogía de la definición 2 de campo 
vectorial con la definición 8 para el caso de 1-formas, esto 
es, tenemos que la definición 8 es la misma que la definición 
2 de campo vectorial salvo que hemos cambiado (RIJ) * por 

Nuevamente como en el caso de campos vectoriales 
tenemos que la condición dada en la definicióon 8 de que 

w(p) E A k {{Rp) *) para cada pGR n equivale a decir que se 
cumple que 7r o w(p) = p para cada p £ R n donde 

es la proyección canónica definida por tt (a) = p si 

a€A*((RJ)*> 

Otra observación es que si w es una k-forma con 
dominio R n entonces por la definición 8 tenemos que 
w(p) es un elemento en A fc (([Rp*). Por otra parte, si 
• • • (e n ) p } es una base de R£ y WÁp), ■ ■ ■ ,<P n (p)} es 
la base dual de {(ej^, . . . (e n ) p } para (R™)* entonces de la 
construcción de una base para A fc ((R™)*) que dimos en la 
sección 2 tenemos que w(p) se expresa como 

»!<• "<ik 



donde los Wi 1 ...i k (p) son escalares en R que dependen 
del punto P. Esdecir, tenemosquelasU^ 1 ..i¿ c sonfunc¡ones 
de R n en R. Al igual que en el caso de campos vectoriales 
este hecho nos permite dar la siguiente definición 



Definición 9. Uno k-formo se Homo continuo o diferencioble 
si codo uno de los funciones Wí^.ai, lo son en el sentido 
usuol. 



Nuevamente como en el caso de campos vectoriales 
las formas se suman y se multiplican punto a punto, mas 
precisamente tenemos 

Definición 10. Seonwyr/dosk-formosyseaeuno l-formo. 
Definimos 

(w+r))(p) = w(p)+r)(p) 

(fw)(p)=f(p)w(p) 
(wAc)(p) = w(p) Ae(p) 

dondep<ER n yf : R n —> Res uno función diferencioble. 

En analogía con el caso de campos vectoriales 
tenemos que w+rj y fw son k-formas y que wAe es una 
(/c+l)-forma. Más aún, si w,rj y e son formas diferenciabas 
(continuas) y / es diferenciable (continua) entonces no 
es difícil demostrar que las formas w+r], fw y wAe son 
diferenciabas (continuas). Notemos también que una 
0-forma es simplemente una función f de R n en R. Así, 
tenemos que fAw = fw. 

Haciendo una ligera modificación de la expresión (3) 
de la sección 1 para Df(p) tenemos la siguiente: 
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Definición 11. Sea f : R n — )► R uno función diferencioble. 
Definírnoslo 1 -formo 

peR n 

de lo siguiente monero, o codo p en R n le osociomos 
lo función lineol df(p) : R n —> R definido por 
df(p)(v p ) = Df{p){v)poroco P dov p enR n r 

Notemos que la definición 11 es en realidad la 
definición formal del operador d para O-formas discutido 
en la sección 1 que manda a cada 0-forma / en la 1 -forma 

df. 

El siguiente lema nos permitirá dar una expresión 
explícita para la 1- forma df de la definición 1 1 y también 
dar la expresión de cualquier k-forma w en términos de los 

dxi 

Lema 1 2. Seonx 1 : R n — >> Ríos proyecciones canónicos poro 
codo i y consideremos los 1 -tormos asociadas dx 1 : ...^dx n . 
Entonces el conjunto 

{dx 1 (p),...,dx n (p)} 

es lo bose duol poro A 1 ^) = (Rp)* asociada a la base 
canónico {(e^p, . . . ,(e n ) p }de[R^. 

Demostración: 

(dxi(p))((e 3 ) p ) =Dx%p)(e j ) 

t-»fe> 

donde la segunda igualdad es por (3) de la seccón 1 . 



Corolario 1 3. por el lema 12 coda k-forma w con dominio en 
R n se puede escribir como 

w = 5j Wi V ..i k dx ll /\' • • Adx ik 

dondew ll ... lk : ^ n -> R } 1 <¿ 1; . . . i k <n. 

En particular si w es una 1 -forma con dominio R n 
entonces w es de la Forma 

w = g t dx x + g n dx n 

donde las gi son funciones de IR n en R, lo cual coincide 
con la definición informal de 1 -forma que dimos en la 
sección 1 . 

Procedamos ahora a obtener la expresión de la 1 -forma 
<i/asoc¡ada a una función diferenciable f ; |R n — >> IR de la 
definición 11. Para esto necesitamos primero el siguiente 
lema: 

Lema 14. 

(dx\p))(v p ) = v % donde v p =v 1 (e 1 ) p H v n (e n ) p 

Demostración: 

dxi(p)(v P ) = dx*(p){j: n . =1 v i (ej) P ) 
= Z P j=í v i dx*(p)((e J ) P ) 

= v l 
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donde la tercera igualdad es por el lema 1 2. 

Teorema 15. Sea j : ^ n _^ [R una función diferenciable, 
entonces 

df= ^-dxi+-.- + ^-dx n 
Demostración: 

(df(p)(v P ) = (Df)(p)(v) 

= Z n i=1 v l Df(p)(ei) 

donde la última igualdad es por el lema 1 4. 

Así tenemos que la expresión para df dada por el 
teorema 15 coincide con la expresión "informal" que 
obtuvimos para <if en la sección 1. 

Procederemos ahora a definir el pull-back asociado a 
una función / : R n R m diferenciable y a una /c-forma 
w con dominio R m . Pero antes necesitamos introducir dos 
definiciones. 

Haciendo una ligera modificación de la definición de 
la diferencial Df(p) dada en la página 2 en la sección 1 
tenemos la siguiente: 

Definición 16. Seo f : R n R m una función diferenciable. 
Para cada p fijo en R n f induce la función lineal 
/♦ : R n P ~+ % definida por /* (v p ) = (Df(p)(v)) m para 
cada v p G Rp : 

Definición 17. Para cada p fijo en R n la función 
/ + : Rp Rfl P ) de la definición 16 induce una función lineal 

/*: A*((R^)*)^A*((r p )*) 

definida de la siguiente manera, si cp£A k ((R^)*), 
esto es 



k - veces 

es una función multilineal y alternante, entonces por 
definición 



Rp ^R 



es el elemento en A k ((r™)*) tal que 

r Win w P ) = v(/*(K)p),...,/*((vfc) P )) 

para cada en 

Usando las dos definiciones anteriores estamos listos 
para presentar la definición de pull-back para una /c-forma: 

Definición 1 8. Dada una k-forma 

w : R m [J A fc ((R-)*) 

pGlR m 

y una función diferenciable f : R n — >• R m definimos la 
k-forma 

/»:R"^U A fc ((R"p)*) 

P eR n 

por (/*H)(p) = />(/(?)) ) £ A fc ((Rp)*) para 
cada p en R n . A la k-forma f* w se le llama el "pull-back" de 
w bajo la función f. 

Es importante notar lo siguiente: si n<m entonces 
R n C R m de manera natural. Si tomamos /como la función 
inclusión i : R n — )► R m en la definición 1 8 y consideramos 
una k- forma diferenciable w en R m entonces i*(w) es 
una /c-forma diferenciable en R n que puede pensarse 
como la restricción de w al dominio R n , aun cuando R n no 
es un abierto en R m . Esta es una propiedad fundamental 
de las formas que en general no comparten los campos 
vectoriales: cualquier forma en R m se puede restringir a 
un subespacio S de R m siempre y cuando tenga sentido 
decir que la función inclusión de S en R m es una función 
diferenciable. En general esto ocurre cuando S es una 
subvariedad diferenciable de R m (consultar [5], pag. 3 
y pag. 9 para la definición de variedad diferenciable y 
subvariedad diferenciable). 

En el siguiente teorema enunciamos las propiedades 
algebraicas del pulí back, para su demostración consultar 
[7] (Teorema 4-8, pag. 83). 

Teorema 1 9. Sea f : R n R m una función diferenciable. 
Así f = {f l ,...J m ) con fi:R n ^ R . Entonces para 
cualesquiera w,w 1 ,w 2 k-formas con dominio R m y 
cualquier función g : R m — >> R diferenciable se cumplen las 
siguientes propiedades: 

a) /*(^)=s; =1 £^- 

(2) /*K+w 2 )=/*K)+/*K). 

(3) r(gw)=(gof)f*(w). 

(4) /*ka™ 2 )=/*K)a/*K). 



k - veces 
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Daremos a continuación la definición formal de 
operador diferencial d para cualquier /c-forma. 

Definición 20. Sea W una k-forma (diferenciable) con 
dominio R n . Entonces por el corolario 13 W se puede escribir 
como 



ik 



w = ^ Wi 1 ...i k dx ll A- • • Adx 1 



donde Wí 1 ...íj £ ' R n Rson funciones diferenciabas. 
Definimos una (k+1)- forma dw con dominio en R n por 

dw = dwi 1 ...i k Adx il A"-Adx ik 



donde 



dx j 



por el teorema 75. 



En el siguiente teorema enunciamos las propiedades 
principales del operador diferencial d, su demostración se 
encuentra en [7] (Teorema 4-10, pág.84). 

Teorema 21. Para cualesquier par de k-formas W^W^ , 
cualquier escalar c en IR, cualquier función diferenciable fde 
R m en R n y cualquier l-formar] se cumple: 

(1) d(w + cWi) =dw + cdw 1 (Linealidad). 

(2) d(wAr¡) = dw A7¡-\-(-l) kl wAdr¡. 

(3) d(dw) = 0(es decir, d 2 (w) = 0). 

(4) f*(dw)=d(f*(w))(d conmuta con el pull-back). 

Definición 22. Una k-forma W es cerrada sidw = 0. 

Una k-forma w es exacta si existe una (k- 1)-forma r¡ tal 
que w = dr\. 

Observación 23. el inciso (3) del teorema 21 dice que toda 
forma exacta es cerrada. 

Observación 24. No toda forma cerrada es exacta, por 
ejemplo en R 2 - {0,0} podemos definir la 1 -forma 



w 



-y—d 



X + 



X 



x 2 +y 2 



-dy 



x z +y ¿ 



entonces se puede demostrar que no existe una 
0-forma f defiida en R 2 - {0,0} tal que df* = w. Para la 
demostración de este hecho consultar [7], pág. 85. 



Sin embargo el Lema de Poincaré nos da condiciones 
suficientes para que una fe-forma cerrada definida en un 
dominio abierto de R n tenga una primitiva: 

Teorema 25. (Lema de Poincaré) Si A C R n es un conjunto 
abierto contraíble entonces toda k-forma cerrada es exacta. 

Demostración: Ver cita bibliográfica (7), pág. 86, o bien 
(3), pág. 23. Para una prueba másgeneral consultar (8), pág. 
155 o también (5), pág. 109. 



Notemos que no es necesario pedir que el dominio de 
definición de una fe-forma cerrada w sea contraíble para 
que U> tenga una primitiva: dado que R n es contraíble 
tenemos, por el lema de Poincaré, que cualquier fe-forma 
cerrada wcon dominio R n tiene una primitiva, es decir, 
existe una (/c -1)- forma r¡ tal que drj = w.S\p es cualquier 
punto en R n , podemos considerar la restricción de w al 
dominio abierto U definido por JJ = [R n -{p} el cual no 
es contraíble. La restricción de w, w\u es nuevamente 
una fe-forma cerrada con dominio U y es tal que tiene una 
primitiva dada por r¡ \ u . 

Dicho de otro modo, si i : JJ R n denota la inclusión 
de U en R n entonces aplicando el inciso (4) del teorema 
21 tenemos que la fe-forma restricción i*{w) con dominio 
U es tal que 

d(i*w) =i*(dw) =¿*(0) = 
es decir, la fe-forma i*(w) es cerrada. Además 

d(i*rj) i* (drj) i*(w) 
es decir la primitiva dei*(w) es la (fe -1 Horma i*(r¡). 

Con este análisis concluimos los objetivos de este artículo. 
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procesos tratan con minerales, son las características fisicoquímicas de éstos las que determinan las 
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jugado un importante papel en el desarrollo nacional, que no existen en nuestro país laboratorios de 
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INTRODUCCIÓN 

A pesar de que desde hace mucho tiempo se conoce sobre 
la aplicación de la mineralogía a los procesos industriales 
(1), éste y otros aspectos de las rocas mineralizadas han 
sido resumidos en un ingenioso término acuñado en los 
últimos años: la geometalurgia (2, 3, 4). 

La geometalurgia es una herramienta que integra 
información geológica, minera, metalúrgica y económica, 
para maximizar el valor actual neto de una mina, al 
mismo tiempo que se minimiza el riesgo operativo 
y técnico (5). Puesto que parte de la geometalurgia 
consiste en la caracterización de los minerales de mena 
y ganga, esto enlaza la variabilidad mineralógica y las 
texturas del depósito con un modelo de bloques en tres 
dimensiones (3D) y la relaciona con los costos. Lo anterior 
permite tener en cuenta la variabilidad de la mena en 
las previsiones de calidad y producción, en el diseño 
de la infraestructura y en el diagrama de flujo a lo largo 
de la vida productiva de una mina. La geometalurgia ha 
tomado gran auge en los últimos años y ha reivindicado a 
la geología económica en el sentido de que revaloriza su 
influencia en la ingeniería de procesos. El mayor trabajo 
de investigación en geometalurgia se hace actualmente 
en algunas universidades australianas, donde se cuentan 
numerosas tesis de doctorado sobre este tema, incluyendo 
el desarrollo de software especializado del tipo de Mineral 
Liberation Analizer (MLA) adaptado al microscopio 
electrónico de barrido (MEB-EDS). 

BASES DE LA GEOMETALURGIA 

La importancia de la geometalurgia basa su principio en 
la comprensión de que las menas son combinaciones de 
minerales, no de elementos químicos y, puesto que todos 
los procesos tratan con minerales, son las características 
fisicoquímicas de éstos las que determinan las condiciones 
de un proceso de recuperación industrial. De una manera 



más específica, a continuación se mencionan algunas de 
las bases que sustentan la geometalurgia: 

Las menas son combinaciones de minerales. Son 
entonces las características de éstos las que determinan las 
condiciones de los procesos industriales. La determinación 
y control de procesos basada sólo en una química general 
puede ser engañosa. 

• Los depósitos minerales están formados por más de 
un tipo de mena. 

• Los minerales de ganga son mucho más abundantes 
que los de mena. Estos minerales pueden determinar, 
por ellos mismos, condiciones de procesos de minado, 
quebrado, molienda, flotación, etcétera. 

En algunas menas el metal importante ocurre en más 
de un mineral. Los minerales tienden a presentar diferente 
comportamiento durante los procesos; las recuperaciones 
de un elemento pueden ser diferentes para diferentes 
minerales. Por otra parte, es el tipo de mineral de mena 
recuperado el que determina por ejemplo, la ley de un 
concentrado que alimenta a una fundición. 

No reconocer las aseveraciones previas, implica 
problemas en los procesos {Le. baja recuperación). De 
donde las mezclas de menas de diferentes unidades 
geometalúrgicas, pueden ser la solución adecuada durante 
la vida de una mina. 

El modelo geometalúrgico es el resultado del 
conocimiento global del depósito mineral en lo referente 
a su mineralogía, química, comportamiento a los procesos 
y recuperación. Esto significa que no deberían ocurrir 
sorpresas significativas durante la vida de una mina. 

INFORMACIÓN QUE PUEDE APORTAR 
LA GEOMETALURGIA 

Toda la información que se deriva de un estudio 
geometalúrgico se puede utilizar en procesos relacionados 
con la evaluación de recursos, los métodos de explotación y 
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los métodos de beneficio de un mineral. A continuación se 
mencionan algunas de estas posibilidades de información: 

Relación de especies minerales, características 

texturales entre ellas, y leyes (Figuras 1 y 2). 

Relaciones entre minerales de mena y ganga. 
Concentración de elementos deseables o 

indeseables durante un proceso. 

Dureza del material (facilidad de molienda). 

Conocimiento de la mejor malla de liberación 

de las partículas minerales (Figura 3). 

Recuperación metalúrgica. 
Características de productos de desecho 

(terreros y jales de plantas de beneficio, escorias 

y polvos de fundiciones). 
Relaciones de la mena con el consumo de 

reactivos. 

Facilidades para la perforación durante la 
explotación o exploración. 
• Facilidades para la fragmentación durante la 
explotación. 

Características de productos de fundición. 




Figura 1. Concentrado de pirita. Fotografía al microscopio 
mineragráfico mostrando la inclusión de un telururo de 
oro (flecha) en un grano de pirita. El tipo de mineral y su 
textura, plantea problemas en la recuperación del oro. El 
grano mide aproximadamente 0.15 mm. Foto: Efrén Pérez 
Segura. 




Figura 2. Imagen al microscopio electrónico de barrido 
(MEB) del mismo grano de pirita de la foto de la figura 1. 
Se observa la misma inclusión de telururo de oro y otras 
de la misma composición no observadas al microscopio 
mineragráfico, algunas siguiendo el plano de crucero de la 
pirita. La composición del mineral es de Te (66 %) y Au (34 
%). Foto y análisis: Francisco Brown. 
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Figura 3. Diferentes grados de liberación de una partícula 
(negro). Los números representan el por ciento de mineral 
de interés en la partícula. Figura tomada de Petruk (6), 
modificada. 

TÉCNICAS DE TRABAJO 

Se considera que para la explotación y beneficio de 
minerales deben definirse unidades metalúrgicas que 
no forzosamente corresponden a los tipos litológicos o 
de alteración que se conocen en un yacimiento dado. La 
definición de una unidad geometalúrgica se debe hacer 
de común acuerdo entre el geólogo, el ingeniero de minas 
y el ingeniero metalurgista y depende del tipo de roca, 
de alteración, de la mineralogía de mena y ganga y de la 
morfología y texturas minerales. Por su importancia, una 
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unidad geometalúrgica debe representar por lo menos un 
5 o 1 0% del volumen de recursos económicos explotables. 
Todos los datos de un estudio geometalúrgico se pueden 
introducir como variables geoestadísticas en el software 
especializado de las minas conocido genéricamente como 
"el modelo". En "el modelo" se fundamenta el plan de 
minado y las variables geometalúrgicas se pueden traducir 
en cálculos económicos por bloques, desde antes de que 
el mineral sea enviado a una planta de beneficio. 

Como toda actividad de investigación que trabaja 
con muestras, lo más importante es la representatividad 
de la muestra que, en este caso, debe ser extrapolable a 
decenas o centenas de miles de toneladas de mineral. 
Por lo tanto, quienes se involucran en un estudio 
geometalúrgico participan activamente desde la selección 
y toma de muestras, hasta un reporte final. Las técnicas de 
trabajo van desde la escala hectométrica hasta la escala 
microscópica. Una metodología general involucra los 
siguientes aspectos: 

Cartografía geológica de detalle en el yacimiento 
(Reconocimiento de litología, tipos y grados de alteración, 
tipos y grados de mineralización). 

Descripción detallada de núcleos de barrenos de 
diamante (Reconocimiento de los mismos parámetros 
anteriores y selección de tramos para muestreo). 

Muestreo de rocas mineralizadas representativas y de 
tipos de mineral para compósitos (Con el objeto de definir 
unidades geometalúrgicas y realizar trabajo de detalle al 
microscopio). 

Estudio de láminas delgadas y secciones pulidas al 
microscopio polarizante (Definición de las características 
mineralógicas, litológicas y texturales al microscopio) 

Estudios de difracción de rayos X, principalmente 
método Rietveld (Composición de minerales específicos, 
principalmente arcillas y proporciones generales no 
identificadas al microscopio petrográfico). 

Estudios al microscopio electrónico de barrido (MEB) 
con EDS (Energy Dispersive System) (Composición específica 
de los minerales y texturas a esta escala). 

Utilización de un programa informático como MLA 
(Mineral Liberation Analizer) o equivalente, adaptado al 
MEB-EDS (Programa que permite un conteo automatizado 
de partículas por composición y cálculo del grado de 
liberación de las partículas). 

Reconciliación entre composición química 
global de muestras y fracciones calculadas con base 
en las composiciones de los minerales específicos y la 
composición modal (mineralógica) de la muestra. 



POSIBILIDADES DE VINCULACION 

La industria minera en México ha jugado un importante 
papel en el desarrollo nacional y el estado de Sonora se ha 
consolidadocomounodelosestadoslíderesen producción 
minera. Por otra parte, la Universidad de Sonora cuenta 
con personal académico altamente calificado en las áreas 
de geología económica, ingeniería de minas y metalurgia 
extractiva. Cuenta también con carreras de licenciatura y de 
maestría en las mismas áreas y con una infraestructura en 
laboratorios que incluye microscopía óptica y microscopía 
electrónica de barrido (MEB-EDS), laboratorios de análisis 
de minerales y software para el modelado minero como 
el que se usa en las minas en explotación. Considerando 
que no existen en el país laboratorios de investigación o 
de servicios que realicen trabajos de geometalurgia, se 
considera que esta disciplina constituye una excelente 
oportunidad para realizar trabajo multidisciplinario de 
vinculación con la industria minera. Para esto es necesaria 
la conformación de un equipo humano que pueda 
ofertar a la industria este tipo de servicios altamente 
especializados. La ventaja de este esquema es que ya 
se cuenta con personal académico especializado, con 
estudiantes a nivel licenciatura y maestría para realización 
de tesis y, lo más importante: no es necesaria la adquisición 
de infraestructura física, ni mucho menos la conformación 
de un aparato burocrático nuevo. 
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DE METEOROS Y METEORITOS 



RICARDO VEGA GRANILLO 



Un meteorito es un cuerpo extraterrestre que sobrevive a su impacto con la superficie de la Tierra sin 
ser destruido. Los meteoritos tradicionalmente se han dividido en dos grandes categorías: 1) Meteoritos 
no-diferenciados o condritas, y 2) Meteoritos diferenciados que incluyen acondritas, meteoritos de hierro y 
meteoritos de hierro-roca. La gran mayoría de los meteoros se agotan o desintegran por el calor antes de tocar 
la superficie de la Tierra, pero algunos alcanzan a llegar hasta ella produciendo un meteorito. Los grandes 
meteoritos pueden golpear a la Tierra con una fuerza considerable, produciendo un cráter de impacto. Los 
cráteres de impacto conocidos presentan una depresión o hundimiento de forma circular u ovalada. En la 
Tierra, la mayoría de los cráteres de impacto más antiguos han desaparecido casi en su totalidad debido 
a la erosión, aunque algunos han sobrevivido. Por los años sesentas, numerosos investigadores, hicieron 
estudios detallados de los cráteres y aportaron una evidencia clara de que habían sido originados por 
impactos meteóricos. Se considera que hace unos 65 millones de años, el impacto de un gran meteorito 
provocó la extinción de más del 75% de plantas y animales terrestres. Se cree que este meteorito cayó en 
lo que actualmente es la Península de Yucatán formando un cráter de unos 180 km de diámetro. En ese 
contexto, el estudio y conocimiento de los meteoritos puede ser fundamental para la supervivencia de la 
especie humana. 
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DEFINICIÓN Y CLASIFICACIÓN 
DE LOS METEORITOS 

Un meteorito es un cuerpo extraterrestre que sobrevive a 
su impacto con la superficie de la Tierra sin ser destruido. 
Mientras está en el espacio es llamado meteoroide. Cuando 
entra en la atmósfera, la resistencia del aire causa que la 
materia del meteorito se caliente y emita luz formando una 
bola de fuego, también conocida como meteoro o estrella 
fugaz. 

Los meteoritos tradicionalmente se han dividido en 
dos grandes categorías: 1) Meteoritos no-diferenciados 
o condritas, y 2) Meteoritos diferenciados que incluyen 
acondritas, meteoritos de hierro y meteoritos de hierro- 
roca (1) (Figura 1). Las condritas están formadas por 
cóndrulos o pequeñas esferas integradas principalmente 
por silicatos, son muy antiguas y se formaron antes de 
que existiera el sistema solar, al juntarse polvo y granos 
de una nebulosa primitiva que originó a los asteroides. Se 
estima que estas rocas forman un 86% de los meteoritos 
que han caído en la Tierra. Las acondritas son materiales 
provenientes de asteroides o planetas en los que ya existía 
una separación o diferenciación de los elementos químicos 
y que, en consecuencia, estaban formados por corteza y 
núcleo. Estos meteoritos más evolucionados forman un 8% 
de los que han impactado en la Tierra. Los meteoritos de 
hierro están compuestos en su mayoría por una aleación 
de hierro-níquel y constituyen el 5% de los meteoritos; y 
los meteoritos de hierro-roca que contienen cantidades 
grandes tanto de material rocoso como metálico, son los 
más raros y sólo constituyen el 1 % de los detectados en la 
Tierra (2). La clasificación más moderna de los meteoritos 
los divide en grupos de acuerdo a su composición química 
y mineralógica, así como por su estructura. 



Meteoritos no identificados 



Condritas 



Condritas Condritas Condritas Condritas 
Carbonáceas Ordinarias Rumuruti deenstatita 



Meteoritos diferenciados 

H 



Acondritas 



Meteoritos de Meteoritos 
hierro-roca de hierro 

, l , 



Acondritas Palasitas Mesosideritas 
Primitivas 



Figura 1 . Clasificación primaria de los meteoritos. 



EFECTOS DE LA CAÍDA DE METEORITOS 

La gran mayoría de los meteoros son cuerpos del tamaño 
de un grano de arena y cuando entran en la atmósfera 
sus átomos se vaporizan, chocando entre sí y contra los 
átomos de la atmósfera, produciendo la luz que vemos en 
las llamadas estrellas fugaces. En cuerpos mayores, como 
de 1 centímetros a algunos metros, la visibilidad se debe 
al calor producido por la presión del objeto en los átomos 
del aire. Este aire comprimido se calienta y entonces 
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calienta los átomos del meteoro. La gran mayoría de los 
meteoros se agotan o desintegran por el calor antes de 
tocar la superficie de la Tierra, pero algunos alcanzan a 
llegar hasta ella derivando en un meteorito. Esto depende 
de su tamaño, su estructura y su velocidad. Entre mayor 
es el tamaño, más probabilidades tendrán de sobrevivir, 
si tienen una estructura débil se fragmentarán y cada 
fragmento se desprenderá más fácilmente, pero si son 
muy veloces, la intensidad del desgaste en la atmósfera 
será mayor y tenderán a consumirse. 

Aproximadamente 500 meteoritos cuyo tamaño varía 
desde el de un frijol hasta pelotas de basquetbol alcanzan 
la superficie de la Tierra cada año, y algunos pueden ser 
mucho mayores. Los grandes meteoritos pueden golpear 
a la Tierra con una fuerza considerable, produciendo 
un cráter de impacto. La clase de cráter dependerá del 
tamaño del meteorito, de su composición y del ángulo de 
entrada. La fuerza de la colisión de un meteorito puede 
producir una destrucción amplia. A lo largo de la historia se 
ha registrado un daño ocasional a viviendas, ganado y aun 
a personas, como es el caso del meteorito de Cheliabinsk 
que recientemente explotó y se fragmentó en los cielos 
de Rusia. En el caso de fragmentos de cometas, los cuales 
están compuestos principalmente de hielo, puede ocurrir 
un daño considerable, aunque en ese caso, no sobreviva 
ningún fragmento del meteorito original. Por ejemplo, en 
Rusia se produjo un impacto de este tipo en la región de 
Tunguska. 



El evento en Tunguska fue una explosión natural que 
ocurrió en una región desolada de Siberia en 1908. Se 
considera que la explosión fue probablemente producida 
por la onda expansiva de un cometa que explotó entre 6 
a 10 km sobre la superficie de la Tierra. La energía de la 
explosión fue posteriormente estimada como la de una 
bomba de TNT de 10 a 15 megatones (3). Se estima que 
la explosión causó la caída de 80 millones de árboles en 
una superficie de 2,1 50 km 2 . El evento de Tunguska es una 
demostración histórica de que un evento catastrófico para 
la humanidad es una posibilidad real. 

Los cráteres de impacto conocidos presentan una 
depresión o hundimiento de forma circular u ovalada. 
En la Tierra, la mayoría de los cráteres de impacto más 
antiguos han desaparecido casi en su totalidad, aunque 
todavía persisten huellas de algunos de ellos. La erosión 
por el agua y el viento, los depósitos de arena, el suelo 
transportado por el viento o el agua y los flujos de lava, 
tienden a oscurecer y ocultar estos cráteres. Por otra parte, 
en la Luna podemos observar numerosos impactos de 
meteoros, algunos tan grandes que son visibles a simple 
vista. Sin embargo, en nuestro planeta aún queda alguna 
evidencia de los grandes impactos de los cuales se han 
reconocido alrededor de 183 (Eorth Impoct Datábase, 
University ofNew Brunswick). 

Un objeto que cae a la Tierra desde el espacio avanza a 
una velocidad mínima de unos 1 2 km cada segundo. Debido 
a que la energía de un objeto depende de su velocidad, un 
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objeto muy grande viajando a tal velocidad puede producir 
explosiones que se parecen a las explosiones nucleares. 
Los aparatos que detectan las ondas sísmicas registran casi 
cada año un impacto multikilotón en la Tierra, usualmente 
dicho impacto ocurre en algún océano. 

Cuando el objeto golpea, puede desintegrarse 
totalmente y los fragmentos, mezclados con los de la roca 
donde impacto, son arrojados hacia fuera, produciendo 
el cráter. Fragmentos del impacto y de la roca impactada 
llegan a derretirse y luego, ya enfriados, pueden 
encontrarse alrededor de algunos cráteres de impacto en 
forma de pequeñas bolas de vidrio a las que se les llama 
tectitas. Partículas del meteoro del tamaño de polvo 
pueden persistir en la atmósfera por varios meses. Esas 
partículas pueden afectar el clima, mediante la dispersión 
de la radiación solar y por reacciones químicas en la 
atmósfera superior. 

METEORITOS Y SU REGISTRO GEOLÓGICO 

Daniel Barringer, que vivió entre 1860y 1929, fue uno de 
los primeros en identificar una estructura geológica como 
un cráter de impacto (4). El cráter estudiado por 
este investigador se encuentra en Arizona, 
al sureste de la ciudad de Flagstaff y es 
llamado "Cráter Meteórico Barringer" 
o 'Cráter Meteoro 1 . En su tiempo, sus 
¡deas no fueron muy aceptadas, y 
cuando lo fueron, no se reconoció 
el hecho de que los impactos 
en la Tierra eran comunes 



considerando la edad del planeta y su registro geológico. 

En 1920 el geólogo americano Walter Bucher estudió 
numerosos cráteres en los Estados Unidos y concluyó que 
habían sido creados por grandes eventos explosivos, pero 
creía que habían resultado de explosiones volcánicas. 
Sin embargo, los geólogos John Boon y Claude Albritton, 
que revisaron los estudios de Bucher concluyeron que 
los cráteres eran probablemente resultado de impactos 
meteóricos. El asunto permaneció más o menos como 
una especulación hasta los sesentas, cuando numerosos 
investigadores, entre los que destaca Eugenio Shoemaker, 
hicieron estudios detallados de los cráteres y aportaron 
una evidencia clara de que habían sido originados por 
impactos, identificando los productos que solamente se 
asocian con impactos, entre los que destaca la presencia 
en la superficie de una variedad de sílice llamada coesita, 
que sólo se produce bajo muy altas presiones (a más de 
80 km de profundidad). Aun así, su trabajo fue discutido y 
atacado, sin embargo, los aterrizajes del Apolo en la Luna, 
aportaron evidencia de que los numerosos cráteres visibles 
en ella eran también producto de impactos meteóricos. 
Debido a que en la Luna los procesos de erosión son 
mínimos, los cráteres han persistido casi de manera 
indefinida. Por lo tanto, era de esperarse que la 
Tierra tuviera al menos la misma proporción 
de impactos que la Luna, volviéndose 
claro que nuestro planeta había sufrido 
muchos más impactos de los que se 
pueden observar actualmente en su 
superficie. 
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METEORITOS Y EXTINCIONES MASIVAS 

En 1980, un equipo de investigadores dirigidos por 
el físico Luis Álvarez y su hijo el geólogo Walter Álvarez 
(5), estudiaban rocas sedimentarias en el límite entre los 
periodos Cretácico y Terciario (límite KT), el cual se ubica 
hace unos 65 millones de años. Estos dos periodos dividen 
dos eras principales de la Tierra: la Mesozoica y la Cenozoica, 
en la que nos hallamos. La separación de estas eras se basa 
en un cambio global de los organismos de la Tierra, tanto 
plantas como animales, incluyendo a los microorganismos 
marinos. Este cambio está marcado por una extinción 
masiva de especies. Más del 75% de las plantas y animales 
perecieron, incluyendo a los dinosaurios. 

Álvarez y sus colaboradores encontraron en el límite 
entre el Cretácico y el Terciario (límite KT) una pequeña 
capa, como de un centímetro de espesor, con una 
concentración del elemento Iridio cientos de veces más 
grande que lo normal (6). Debido a que este elemento 
es común en muchos asteroides, los investigadores 
propusieron que esta capa estaba relacionada con el 
impacto de un asteroide (5). Supusieron que el impacto 
del asteroide provocó su fragmentación y conversión 
en polvo, dicho polvo ascendería a la atmósfera y luego 
caería como una lluvia fina, depositando la capa. A partir 
de la cantidad de Iridio presente en el límite KT, los 
investigadores estimaron que un asteroide de 10 a 14 km 
de diámetro había colisionado con la Tierra. Esas primeras 
publicaciones impulsaron la investigación en otras partes 
del mundo. Así, la capa rica en Iridio se ha encontrado en 
más de cien sitios distintos por todo el mundo. Una prueba 
más de que esta capa provenía de un impacto es que 
contiene coesita, forma de silicio creada sólo a enormes 
presiones. 

La explosión que resultó del impacto de tan gran 
meteoro debió ser cientos de millones de veces más 
devastadora que la de las armas nucleares más poderosas 



que se han explotado. Pudo haber producido un 
huracán de furia inimaginable y haber lanzado 
cantidades masivas de polvo y vapor en la 
atmósfera superior y aún al espacio, formando una 
tormenta de fuego global, cuando los fragmentos 
incendiados de la explosión cayeron de regreso 
a la Tierra. También pudo producir una nube de 
polvo que envolvió al mundo tapando la luz solar 
por meses, lo que debió traer como consecuencia 
una disminución de la fotosíntesis afectando así las 
cadenas alimenticias, incluso las del fitoplancton 
marino, de las cuales los dinosaurios eran el último 
eslabón. 

Un problema de la hipótesis de Álvarez y otros 
(5, 6) fue que no se había encontrado un cráter que 
coincidiera con ese evento de fines del Mesozoico. 
En 1978, investigadores que trabajaban para 
PEMEX, descubrieron mediante estudios geofísicos, 
una enorme estructura circular de casi 180 km 
de diámetro (7). La mitad de dicha estructura se 
encuentra en la parte norte de la Península de 
Yucatán (Figura 2), la otra mitad debajo del mar y su centro 
es cercano al poblado de Chicxulub. Los investigadores 
presentaron un reporte en una conferencia, pero fue poco 
atendido debido a que Yucatán está formado por rocas 
relativamente recientes las cuales cubren de forma más 
o menos horizontal a las rocas más antiguas, por lo tanto, 
no era posible tener una evidencia directa del impacto 
meteorice Por esta razón, los investigadores buscaron 
muestras de un antiguo barreno de PEMEX y hallaron rocas 
que exhiben claras evidencias del impacto (8). En 1995 
otros profesionistas que estudiaron imágenes de satélite 
pudieron descubrir una serie de cenotes o hundimientos 
de las rocas superficiales, que en la superficie seguían el 
contorno de la estructura a profundidad y consideraron 
que ocurrieron por el hundimiento de las paredes del 
cráter (9). Pudo entonces establecerse que el cráter tenía 
un diámetro total de cerca de 1 80 km. 
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Figura 2. Ubicación del impacto meteórico en Chicxulub. 
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Muchos estudiosos de la vida antigua de la Tierra 
o paleontólogos están de acuerdo en que un asteroide 
golpeó al planeta hace unos 65 millones de años, pero 
muchos no están de acuerdo en que dicho impacto haya 
sido la única causa de las extinciones. La edad del cráter de 
Chicxulub se ha establecido en unos 300,000 años antes 
del límite entre las eras Mesozoica y Cenozoica (10). Este 
hallazgo sugiere que uno o varios impactos pudieron 
contribuir, aunque no causar directamente la extinción 
masiva ocurrida a fines de la Era Mesozoica. 

REFLEXIÓN FINAL 

Al principio del nacimiento de la Geología como una 
ciencia, a mediados del siglo XVII, había una discusión 
filosófica acerca de la formación de la Tierra. Una hipótesis 
proponía que se había formado por una serie de grandes 
catástrofes, mientras que otra proponía que se debía a 
procesos normales, que actúan hoy y han actuado de 
igual modo, a través de miles de millones de años, a esta 
última teoría se le llama Uniformitarismo. Considerando 
el impacto de grandes meteoritos y el registro geológico, 
se puede pensar que ambas teorías actúan de manera 
conjunta. La verdad es que los impactos meteóricos 
pueden visualizarse desde el punto de vista uniformitarista 
como algo que ha ocurrido a través de la historia planetaria 
y seguirá ocurriendo en el futuro, lo que nos obliga a estar 
prevenidos y ampliar la investigación para preveer posibles 
impactos. 
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USOS Y OCURRENCIA DE LOS PRINCIPALES 
METALES QUE SE PRODUCEN EN SONORA 



EVA LOURDES VEGA GRANILLO*, ALMA PATRICIA SÁMANO TIRADO 



México posee gran cantidad de recursos minerales producto de la naturaleza, por lo que su explotación a 
través de la minería y el cuidado del medio ambiente deben de coexistir de una manera armoniosa. 

A través del tiempo la minería ha desempeñado un papel fundamental en la vida económica del país, ya 
que su desarrollo ha permitido que en la actualidad sea la tercera fuente generadora de ingresos. Asimismo, 
la minería permite que se desarrollen actividades, como la agrícola y la industrial (automotriz, acerera, 
cinematográfica, etcétera), entre otras. A nivel mundial, México ocupa el primer lugar en producción de 
plata, existiendo otros minerales metálicos importantes, como el oro, plomo, hierro, zinc, cobre, uranio y 
torio. 

Se considera a la minería como una actividad económica primaria debido a que los minerales se 
toman directamente de la naturaleza, encontrándose sobre la superficie o bien en el subsuelo a diferentes 
profundidades. Algunos de los principales centros mineros se localizan al norte del país, pero en los últimos 
años, Sonora se ha revelado como el estado productor de minerales más importante de México. Los 
principales minerales son: cobre y metales preciosos, como el oro y la plata. Por ello en el presente artículo, 
se mencionan los usos más comunes de estos metales, así como su ocurrencia en el estado, y su importancia 
económica a escala regional. 
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USOS DEL ORO Y LA PLATA 



El oro (Au) es el metal precioso por excelencia, es de color 
amarillo, brillante, pesado, maleable y dúctil. No reacciona con 
la mayoría de los productos químicos, pero sí es sensible al 
cloro y al agua regia (ácidos nítrico y clorhídrico concentrados). 
Es considerado por algunas personas como el elemento más 
bello, es tan dúctil y maleable, de tal forma que una onza (28.35 
g) de oro podría formar una sábana de 28 m 2 . También es un 
excelente conductor del calor y de la electricidad, no le afecta 
el aire, y se usa principalmente en la acuñación de monedas y 
en la joyería (1). 

La plata (Ag) es un metal de transición de color blanco, es 
un poco más duro que el oro y, al igual que éste, es brillante, 
blando, dúctil y maleable. Tiene la más alta conductividad 
eléctrica y conductividad térmica de todos los metales, pero 
por su alto precio no se usa en aplicaciones eléctricas. La plata 
es un mineral escaso en la naturaleza, se le puede hallar como 
plata pura o formando minerales, principalmente sulfuros (2). 

Desde el año IV a.C. el oro fue utilizado en artefactos 
mortuorios y en adornos como sombreros. En esos tiempos se 
creía que ingerir alimentos en platos de oro podría prolongar 
el tiempo de vida y retardar el envejecimiento de las personas. 
Durante la gran peste negra en Europa algunos alquimistas 
pensaron en curar enfermos dándoles a tomar oro finamente 
pulverizado. 

El uso de estos metales como adornos es una práctica 
milenaria, donde el orfebre se dedica a elaborar o crear figuras 
u objetos decorativos (Figura 1). Al elaborar joyería fina se 
reconocen tres tipos de oro: el oro alto o de 1 8 kilates (1 8 partes 
de oro por cada 6 de otro metal), oro medio o de 14 kilates (14 
partes de oro por 1 de otros metales), y oro bajo o de 1 kilates 
(1 partes de oro por cada 1 4 de otros metales). Dependiendo 
de con qué está aleado, el oro puede tener diferentes colores, 
pudiendo ser oro amarillo, oro rojo, oro rosa, oro azul, etcétera 
(3). 



Figura 1 . Ornamentos aztecas de oro. 




Pectoral azteca Aretes de oro 

http://www.raulybarra.com/notijoya/archivosnotijoya7/7joyeria_oro. 
azteca.htm 
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La plata también se usa para acuñar monedas ya que 
es muy dúctil y maleable, pero un poco más dura que el 
oro; tiene la más alta conductividad eléctrica de todos los 
metales, incluso mayor que la del cobre, pero no se usa de 
forma masiva por su alto precio. Se calcula que alrededor 
del 70% de la producción de plata se utiliza en la industria 
(fotográfica, química, médica, y electrónica), quedando 
el resto para fines monetarios y orfebres. Con la plata se 
elaboran armas blancas como espadas, lanzas, puntas de 
flecha, con el yoduro de plata se crea lluvia artificial (Figura 
2), se utiliza en aleaciones para piezas dentales, en baterías 
eléctricas de plata-zinc, etcétera (2). 

La plata puede encontrarse en forma nativa o 
combinada con azufre, arsénico, antimonio o cloro. Se 
obtiene de minas de cobre, cobre-níquel, oro, plomo y 
plomo-zinc de México, Canadá, Estados Unidos y Perú, 
siendo México el mayor productor mundial, habiendo 
extraído 2,747 toneladas en el año 2000 equivalentes al 
1 5% de la producción mundial de plata. 



Figura 2. Usos de la plata. 




Espadas de plata-bronce 




Lluvia artificial con yoduro de plata 

http://www.darinveracruzano.com/wp-content/up- 
loads/201 1/04/lluviaartificial.jpg 



En la antigüedad se creía que el oro y la plata eran 
regalos de la naturaleza formados por la influencia de la 
luna y el sol, con propiedades que no tenían otros metales; 
de aquí se derivan los primeros intentos de los alquimistas 
de transformar otros metales en oro y plata. 

El uso más generalizado del oro es para acuñar 
monedas (aleado con plata y/o cobre), por su resistencia a 
la corrosión se usa en la industria, la joyería y la electrónica 
(Figura 3). Este metal se ha usado desde el inicio de la 
electrónica en estado sólido, ya que con él se hizo el primer 
transistor. Casi todos los aparatos electrónicos actuales 
tienen oro, el que se utiliza en forma masiva principalmente 
en contactos con chapa de oro y alambres conductores. 
El uso del oro se prefiere sobre otros materiales por su 
gran resistencia a la corrosión, porque es muy dúctil, 
tiene la capacidad de deformarse sin romperse, tiene baja 
resistencia eléctrica y una elevada conductividad térmica. 
Aunque su costo es mayor, su uso en pequeñas cantidades 
es justificable y necesario (4). 



Figura 3. Uso del oro en electrónica. 




A pesar de que el oro y la plata no son nutrientes para 
nuestro cuerpo, su uso podría ayudar eliminando toxinas 
del cuerpo y retardando el envejecimiento de la piel, por 
ejemplo. Se puede afirmar que el oro no es tóxico aunque 
puede producir alergias en la piel, y permanece inalterable 
en los líquidos biológicos como la sangre. 

Un uso exótico que se le está dando en nuestros días 
a estos metales es como ingredientes de lujo en platillos, 
tanto en spray como en láminas o polvos, logrando 
presentar la comida de una manera sorpresiva y sofisticada 
(Figura 4). La Unión Europea considera como aditivos 
alimentarios autorizados y permite el uso en alimentos 
de algunos metales, como el oro, la plata, el titanio y el 
aluminio (5). 
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Figura 4. El oro y la plata en alimentos. 






Pan de dulce con plata Tomates con oro Oro en vino 

http://rvfconsultores.blogspot.mx/201 1/07/tendencias-alimentos-banados-en-oro-y.html 



Aunque la plata no es tóxica, la mayoría de sus sales 
son venenosas y pueden ser carcinógenas, sus compuestos 
pueden ser absorbidos por el sistema circulatorio, 
depositarse en diversos tejidos y provocar argirismo 
(coloración grisácea de la piel y mucosas que no es dañina 
pero si antiestética). Los efectos de la plata en la salud no 
se conocen; desde tiempos de Hipócrates se sabe de su 
efecto germicida y se comercializan hoy en día remedios 
para una gran variedad de dolencias pero no hay estudios 
clínicos que demuestren su utilidad como antibiótico (2). 

USOS DEL COBRE 

El cobre (Cu) es un metal de transición de color rojizo y 
brillo metálico, que junto con el oro y la plata forman la 
llamada "familia del cobre". Después de la plata, es uno de 
los mejores conductores de electricidad. Debido a su alta 
conductividad eléctrica, y ya que es muy dúctil y maleable, 
es el material más usado para fabricar cables eléctricos, así 
como componentes eléctricos y electrónicos. 

Se halla en una gran cantidad de alimentos como 
ostras, mariscos, carne, nueces y en el agua potable, por 
ello, es casi imposible que carezcamos de ese metal (Figura 
5). El Cu actúa, junto a la vitamina C, en la formación de 
colágeno de tejidos conectivos y de hueso. 

López y otros (6) analizaron la concentración de zinc, 
fierro y cobre en alimentos de consumo frecuente en 
México (104 de origen vegetal y 32 de origen animal); el 
cobre no fue hallado en todos los alimentos, siendo la 
máxima concentración reportada de 3.37 mg/100 g en 
muestra de hígado de res. 

Figura 5. El cobre está contenido en diversos alimentos. 




El cobre se usa en aleaciones que en general tienen 
menor conductividad eléctrica, como bronce (Cu-Sn), 
alpaca (Cu-Ni-Zn) y latón (Cu-Zn), entre otras. Fue uno 
de los primeros metales, junto con el bronce usados por 
el hombre, tanto que en la historia a esas épocas se les 
conoce como Edad del Cobre y Edad del Bronce. Ha servido 
para hacer utensilios, motivos ornamentales, puntas de 
flecha, herramientas, armas, en la construcción, objetos 
decorativos, etcétera (Figura 6). 

Figura 6. Usos del cobre. 




Utensilios 




Herramientas 
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El desarrollo de numerosos equipos eléctricos usando 
Cu, tales como generadores, motores y transformadores, 
le ha dado una importancia destacada en la humanidad y 
ha aumentado su demanda grandemente. La aleación con 
otros metales posterior a la revolución industrial también 
aumentó su uso. 

El Cu es el material más usado para fabricar cables 
eléctricos, tanto para la industria como en los hogares. 
También se utiliza en los cables telefónicos, que son la 
base para la comunicación por Internet. Su uso incide 
asimismo, en los medios de transporte (radiador, frenos, 
cojines, cables, motores eléctricos, etcétera), elaboración 
de monedas, maquinaria, e instrumentos musicales, entre 
otros. 

OCURRENCIA EN SONORA 

La geología y tectónica de Sonora permitieron el 
emplazamiento de ricos cuerpos minerales, cuya 
explotación lo ha colocado como uno de los estados 
mineros más importantes de México. Al ser la minería una 
actividad primaria, permitió el asentamiento de ciudades 
en la antigüedad y sentó las bases del Sonora moderno. 
En este contexto, la producción de minerales metálicos 
y no metálicos pone al estado como líder nacional, 
constituyendo importantes fuentes de empleo y otros 
ingresos económicos. Por ejemplo, Sonora es líder nacional 
en la producción de cobre de las minas de Cananea y 
Nacozari, así como de oro de dos minas de clase mundial 
ubicadas en los municipios de Caborca y Sahuaripa (7) 
(Figura 7). 

Figura 7. Minas de cobre y oro en Sonora. 




Espectacular tajo de Cananea, 
principal mina de cobre de México 
http://www.guiarte.com/fotos-usuarios/20011.html 





Tajo de la mina de cobre de Nacozari, la 
segunda en importancia del país 
http://www.grupo-vysisa.mx 




Panorámica de la mina de oro 
La Herradura en el Mpio. de Caborca 

http://farm4.staticflickr.com/3207/3087833491_d3f3b77eca.jpg 




Mina de oro Mulatos, Mpio. de Sahuaripa 

http://www.insetecmexico.com/portal2/notacompleta. 
php?id=361 
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Lo que extraen las minas de Sonora al año equivale a 
alrededor del 50% de la producción nacional y se calcula 
que es más de 31,500 millones de pesos. Esta producción 
incluye el 84.0% del total de cobre y el 28.6% del total 
de oro producidos en México, por citar algunos de los 
minerales más importantes. 

La economía de algunos municipios como Nacozari, 
Cananea, La Colorada, Cumpas, Sahuaripa, Álamos, 
Caborca, Altar, Trincheras, Santa Ana, Magdalena, entre 
otros, está basada en gran parte en la minería. Esta rama 
establecida da empleos directos a más de 1 3,000 personas, 
o sea, más del 2% de los empleos disponibles, cuyos 
sueldos suman alrededor de 1,500 millones de pesos al 
año, agregándole que los trabajadores mineros son de los 
mejor remunerados (7). 

La riqueza minera, el dinamismo del precio de los 
metales, una política laboral estable y un gobierno que 
apoya la llegada de nuevas compañías, hacen que la 
minería se desarrolle de manera notable en Sonora. 
Actualmente, el estado produce el 32.3% de la producción 
minera del país; es líder en la extracción de oro (31.1% 
del total), de cobre (79%), y de plata junto con Zacatecas, 
Chihuahua y Durango (79% del total). En cuanto a minerales 
no metálicos, Sonora aporta el 100% de la producción 
de molibdeno, grafito, carbón antracítico, wollastonita y 
molibdeno (8). 

AÑOS RECIENTES DE LA MINERÍA 

Por todos es conocido que Sonora es un estado 



tradicionalmente minero, pero a principios de la década 
de los noventa, hubo un auge por buscar yacimientos 
de oro, ya no las pequeñas vetas de alta ley, sino todo lo 
contrario: depósitos de gran volumen aunque de baja 
ley. Este auge incentivó la llegada de grandes empresas 
mineras, sobretodo canadienses y norteamericanas, que 
se enfocaron en desarrollar las minas de La Colorada, Santa 
Gertrudis, Amelia, La Choya, Lluvia de Oro, El Boludo, San 
Francisco y La Herradura, algunas de ellas agotadas en la 
actualidad (Figura 8). 

En años recientes se han evaluado distintos proyectos 
mineros, entre los que se destaca la búsqueda del oro de 
Mulatos y El Cháñate, del cobre de Mariquita, Milpillas, Luz 
del Cobre y Piedras Verdes, de la plata de Álamo Dorado, y 
del molibdeno del Crestón (9). 

A manera de conclusión se puede decir que, a la 
destacada producción de cobre en Sonora, se le ha sumado 
en los últimos diez años, la extracción de minerales 
metálicos como oro, plata, molibdeno, y la de no metálicos 
como wollastonita y carbón antracítico. Por ejemplo, 
respecto al oro, el territorio estatal es en este momento 
el principal productor nacional, y en cuanto a plata, 
mantiene volúmenes de producción importantes gracias 
a la fundición y refinería del complejo minero La Caridad. 
Las exploraciones geológicas reportan la existencia de 
otros yacimientos poco o nada explotados de wollastonita, 
zeolita, perlita, talco, arcilla, mica, carbonato de sodio, 
borato, mármol, travertino y feldespato, así como de rocas 
utilizadas en la industria de la construcción, como canteras, 
calizas, puzolana y morusa (7). 
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Figura 8. Principales minas de Sonora. 



1 La Herradura (oro) 

2 El Cháñate (oro) 

3 Mulatos (oro) 

4 Cerro Colorado (oro) 

5 El Mudo (oro) 

6 El Ombligo (oro) 

7 Mexicana de Cananea (cobre) 

8 Mexicana de Cobre (cobre) 

9 María (cobre) 

10 Milpillas (cobre) 

1 1 Piedras Verdes (cobre) 

12 La Escuadra (plata) 

13 Alamo Dorado (plata) 

14 Yavaros(sal) 

15 Salina de Lobos (sal) 

16 Salinera Kino (sal) 

1 7 Molimex (molibdeno) 

1 8 Ca/era Nacozari (cal) 

1 9 Calhidra de Sonora (cal) 

20 Productos Calcáreos (cal) 

21 Cementos Portland 

22 Cementos del Yaqui 

23 Fierro Carbó (fierro) 

24 Fíoty? Bacanora (fierro) 

25 El Volcán (fierro) 




El Yeso (yeso) 
El Taymunco (yeso) 
El Gallardofyeso) 

Bacadehuachi (Sulfato de sodio) 1 5 fc 
La Guásimas (arcilla) 20 11,, 

Porcelanite (feldespato) ¡-luyi-i,;, " C1 ": 

Pí/ara (wallastonita) 
Cobachi (barita) 
Los Aleones (perlita) 
Piedra Laja Sahuaripa yArivechi 
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Todo este auge minero ha producido a su vez, la llegada 
de varios cientos de profesionistas de las geociencias e 
incrementado la matrícula en las carreras y posgrados de 
geología, minería y metalurgia del estado. 

COMENTARIO FINAL 

Sonora es reconocido por sus diversas actividades 
económicas: agricultura, ganadería, industria, pesca y 
minería, siendo esta última proveedora de materia prima 
para la elaboración de productos industriales. Los minerales 
que se producen en Sonora y los extraídos en otros lugares 
del país a través de la historia, le han permitido al hombre 
tener un mejor nivel de vida, ya que sus usos tan variados, 
le han servido para comer, adornarse, sobrevivir, luchar, y 
en general, progresar en muchos aspectos vitales. 

Por todo ello, se considera importante no sólo conocer 
la ocurrencia y utilidad de los minerales producidos en 
Sonora, sino también cuantificar dichos recursos. 
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ANIVERSARIOS ASTRONÓMICOS 
EN CANANEA, SONORA 



RAÚL PÉREZ -ENRÍQUEZ* 



El 21 de marzo de 1913, en la ciudad de México, nació Guillermo Haro Barraza 
y el 27 de abril de 1988, se inauguró el Observatorio Astrofísico Guillermo Haro, 
en Cananea, Sonora. Estos dos muy importantes aniversarios astronómicos se 
conmemoraron el pasado 27 de abril. Aquí se da cuenta de la relevancia de ambos 
festejos, de la fructuosa vida del homenajeado y de las actividades desarrolladas 
ese día allá donde está instalada "Una mina de estrellas". 
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INTRODUCCION 

Fue una fresca mañana la que nos recibió en Cananea, 
la cual no podía ser más propicia para conmemorar dos 
importantes aniversarios en el campo de la Astronomía 
de nue se celebraron. Por un lado, las brillantes mañanas 
de abril, cuando la duración del día comienza a crecer, 
tienen aires de frescura que aún se dejan sentir; por 
el otro, los jardines de la Tasa Green" muy cercana a la 
Sierra de la Mariquita, donde se asienta el telescopio de 
2.1 m, en el centro de la entrañablemente sonorense 
ciudad de Cananea, eran el ambiente y marco de los 
festejos. Ahí, desde donde se abre una ventana hacia el 
resto del universo observable, nos recibieron los directivos 
encabezados por el director, investigadores y trabajadores 
del Instituto Nacional de Astrofísica, Óptica y Electrónica 
(INAOE) para homenajear 
a Guillermo Haro y uno 
de los frutos de su obra, el 
Observatorio Astrofísico 
que lleva su nombre. 

En efecto, los dos 
aniversarios a los que 
hacemos referencia son: 
El centenario del natalicio 
de Guillermo Haro 
Barraza y, los 25 años del 
Observatorio Astrofísico 
Guillermo Haro (OAGH). El 
primero, recordado por el 
Ing. César Arteaga Magaña 
quien acompañara al Dr. 
Haro en las etapas de la 
construcción; el segundo, 
entrelazado con el primero 
haciendo énfasis en la 

entrega de reconocimientos a quienes han vivido ya más 
de 20 años bajo el techo del OAGH. 

EL ASTRÓNOMO GUILLERMO HARO 

Hace cien años, el 21 de marzo, en un barrio la Ciudad de 
México nació quien llegaría a ser, si no el más importante 
astrónomo de México, si uno de sus máximos exponentes. 
En el seno y con el apoyo de la familia Haro-Barraza, 
Guillermo avanzaría en sus estudios de leyes en la UNAM, 
pero alcanzaría un doctorado en la que significaría su pasión 
al descubrirlo, la Astronomía. Después de pasar una buena 
parte de su vida en el Observatorio Astrofísico Nacional 
de Tonantzintla (OANTON), viajó a muchos observatorios, 
llegando incluso al de la Universidad de Harvard. Agudo 
investigador se le reconoce mundialmente por múltiples 
trabajos pero muy especialmente por el descubrimiento 
de los objetos, que en su honor y de un colega, llevan 
el nombre de objetos Haro-Herbig (o HH). Sería largo 
enumerar otras contribuciones del Dr. Haro al desarrollo 
de la ciencia en México; sin embargo, es importante 
destacar que se inició en la astronomía como ayudante 




en el OANTON que vendría a convertirse en el Instituto 
Nacional de Astrofísica, Óptica y Electrónica (INAOE), 
referente nacional de cara a la Astrofísica Internacional. 
Cabe recordar que al fundador del dicho observatorio Luis 
Enrique Erro, le siguió el Dr. Guillermo Haro quien junto al 
Prof. Luis Rivera Terrazas, el Dr. Luis Munch y el astrónomo 
Enrique Chavira, se encontraría entre los científicos más 
destacados de su época. Hacia fines del año 1954, este 
último y el homenajeado descubrieron un cometa que se 
conoció como Cometa Haro-Chavira. En 1953, Guillermo 
Haro recibió la Medalla de Oro "Luis G. León" de la 
Sociedad Astronómica Mexicana; y, en 1962, la Medalla 
Honorífica de la Academia de Ciencias de Armenia; todo 
ello sólo para recordar algunas de sus muchas distinciones. 
Tan solo un año después de recibir la medalla, el Dr. Haro 

se hizo merecedor del Premio 
Nacional de Ciencias. Falleció 
un 27 de abril, pero de 1988 
en la misma ciudad que lo vio 
nacer (1). 

Quisiéramos describir los 
rasgos de su compromiso con 
el conocimientoy con México, 
basta para ello, recordar 
algunas de las expresiones 
que por los años 50's 
expresaba el Dr. Haro cuando 
recibía reconocimiento 
como miembro del 
"Colegio Nacional"; y que 
muestra su espíritu de 
búsqueda, sus virtudes, que 
lo harían emprender los 
proyectos como el que este 
observatorio representa; son 
frases tomadas de su discurso inaugural denominado "En 
el cielo y en la Tierra""(2): 

"...Se dice que las ciencias florecen sólo en el seno 
de los pueblos industrializados y ricos, y que nosotros 
somos pobres y aún usamos el arado egipcio. Mas los 
pueblos ahora prósperos fueron ayer atrasados y pobres 
y, sin embrago, en ellos brotó un poderoso movimiento 
científico y tecnológico creador de riqueza, de bienestar, 
de dinamismo y poderío, y también - por qué no decirlo- 
de abuso y catástrofe". 

"Al dirigirnos al cielo, a ese cielo formado de la 
misma materia que el astrónomo, pensamos en la unidad 
de la física de nuestro universo y en la necesidad de 
entendimiento, de colaboración y de paz que deben reinar 
entre los hombres". 

"La ciencia, desde luego, no es panacea que resuelva 
todos los problemas ni conteste todas las preguntas... 
El hombre de ciencia tiene que ser consciente de las 
preguntas que no debe hacer, y sobre todo, de las 
respuestas que no puede esperar...". 

No cabe duda de que su entrega era a la ciencia pero 
sobre todo a la humanidad. 
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EL OBSERVATORIO ASTROFÍSICO 
GUILLERMO HARO 

En las palabras que acabamos de reproducir se 
ven reflejados los rasgos propios de quien habría de 
apoyar decididamente la construcción del Observatorio 
Astrofísico de Cananea, conocido ahora como Observatorio 
Astrofísico Guillermo Haro (Figura 1). Con la instalación, 
en la cima de la Sierra "La Mariquita", de un telescopio de 
2.1 2 m de longitud, a una altitud 2,480 m sobre el nivel del 
mar, nuestro país se incorporó de lleno a la investigación 
de frontera de la Astronomía y Astrofísica de aquel tiempo, 
1 988. Este telescopio que ha visto pasar por su plataforma 
levadiza primero y por su centro de control después, a 
innumerables e importantes astrónomos y astrofísicos de 
México y del resto del mundo, ha permitido la observación 
de fenómenos estelares relevantes. 

La operación del telescopio ha venido mejorando 
paulatinamente. Con la instalación del espectrógrafo 
alemán LFOSC, en 1991, empezó la etapa de observación 
de objetos ROSAT, fuentes de Rayos X. Posteriormente, 
se instaló el espectrógrafo Boller & Chivens y, más 
recientemente, en 2002, se puso en operación la cámara 
CANICA basada en imagen infrarroja de área ancha con 
HgCdTe de Rockwell (HAWAII por sus siglas en inglés). Así, 
el observatorio ha ¡do evolucionando, de manera que se 
mantiene en la frontera. Hoy, el telescopio está incorporado 
a la familia de instrumentos con visión mejorada del 
cercano infrarrojo. De hace un par de años a la fecha, la 
presencia del espectrógrafo de fibras ópticas Eber-Fastie 
ha llevado a lograr imágenes de ultra alta resolución con lo 
que la vida útil del telescopio se ve alargada (3). 




figura 1 . Observatorio Astrofísico Guillermo Haro, en 
Cananea, Sonora. 

Debemos destacar el papel relevante y simbólico de 
la antigua «Casa Greene», en Avenida Sinaloa No. 25 de 
Cananea, Sonora, sede administrativa del OAGH, en donde 
se hicieron los festejos que aquí se relatan; la casa ha 



tenido un papel muy importante en la vida de este centro 
de investigación. Esta edificación es la que guarda entre 
sus tablones y barandas, techos y ventanas, la historia y 
datos del observatorio y ha dado hospedaje a los visitantes 
y astrónomos todos estos años (4). 

LOS FESTEJOS 

En los jardines de esta simbólica casa, a manera de área 
de exposiciones, se instalaron demostraciones de los 
productos que ofrece el INAOE, telescopios por los cuales 
los asistentes pudieron ver el Sol, durante el día; y la mina 
de la que hablaba el Dr. Haro por la noche. También, se 
pudieron observar experimentos de ciencias realizados por 
estudiantes del Instituto Minerva; y la muestra de carteles y 
dibujos elaborados por estudiantes de las categorías Nova 
y Supernova que respondieron de manera entusiasta a la 
convocatoria "Cananea, una mina de estrellas". 

La UNISON, también, tuvo presencia con exposiciones 
y talleres. Por cierto, se montaron mesas de trabajo en las 
cuales niños, jóvenes y adultos podían armar modelos de 
la Luna, o del Gran Telescopio Milimétrico, maqueta esta 
última del gran proyecto que en unos meses estará en 
pleno funcionamiento en la cima del volcán Sierra Negra, 
muy cerca del Pico de Orizaba, en Puebla. 

Para completar y cerrar con broche de oro la ceremonia 
conmemorativa de los aniversarios astronómicos, 
se presentaron dos cantantes sonorenses quienes 
interpretaron piezas de Giuseppe Verdi por celebrar su 
bicentenario de nacimiento, así como piezas de autores 
mexicanos. Las voces de la soprano Nadia Lamadrid y del 
tenor Adrián Lozano emocionaron a los asistentes al seguir 
los acordes del pianista Jaime Verdugo que los acompañó. 

Viandas y aguas frescas cerraron los festejos de medio 
día, pero la celebración, con la observación nocturna y 
convivencia, habrían de rematar los aniversarios. 

La sociedad sonorense no puede dejar de conocer 
sobre estos aniversarios pues representan una de las 
muchas vetas que dan riqueza a la vida de nuestro estado. 
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Ya estaba por llegar al observatorio. Su rostro denotaba 
cansancio pero sus ojos proyectaban una rara mirada pues 
casi no había dormido por la excitación. La noche anterior 
estando en la cúpula de observación astronómica había 
avistado una nueva luz en el cielo, el cual conocía como 
la palma de su mano. Tantos años de contemplación le 
daban el conocimiento de la colocación en el firmamento 
de las miríadas de lucecitas estelares, según el paso de las 
estaciones del año. 

Había soñado repetidas veces que en la negrura 
aterciopelada de la noche una nueva luz se prendía en la 
bóveda celeste, había sido sorprendente para él visualizar 
en la límpida oscuridad de la noche anterior, totalmente 
despejada de nubes, la luz oteada en sus sueños: era como 
surgida de una intensa hoguera. Sabiendo hacia dónde 
mirar, la podía avistar casi sin esfuerzo, capacidad que le 
era muy útil para su ocupación de astrónomo. Caminando 
apresuradamente al observatorio pasó por la hilera de 
estatuas yacentes y tuvo la impresión de que en lugar de 
estar mirando hacia el horizonte lejano, todas miraran 
hacia el infinito en dirección a la novel llama. 

Al llegar a las empinadas escaleras de piedra que 
llevaban a la entrada del observatorio, subió presuroso 
para ingresar a la edificación. Ascendiendo las escaleras 
de caracol llegó a su amado lugar de trabajo, la cúpula del 
observatorio. 

Ansiosamente esperó Ya casi iba ser el momento en 
que en la ventana de observación quedara encuadrado el 



nuevo fulgor que viera la noche anterior. Seria la prueba 
concluyente que su visión no había sido una continuación 
de sus sueños. Detectarla en el momento esperado 
significaba que realmente estaba colocado en un lugar 
del espacio en que anteriormente solo había tinieblas, 
entonces podría dar la buena nueva de que el nuevo 
resplandor prometía beneficiosos augurios para la ciudad. 

Aguzó la vista y se le aceleró el corazón cuando el 
resplandor esperado apareció, ya tenía una mayor certeza 
de que la nueva estrella iba permanecer para siempre e 
imaginó a su joven ayudante observándola, cuando fuera 
viejo como él ahora y aun después cuando ya hubiera 
tomado el camino a Xibalba bé. Incluso varias generaciones 
después la seguirían observando y recordarían a quien 
tuvo la fortuna de avistarla cuando esta nueva hoguera 
celeste o Tzab-ek se había encendido. Lo recordarían a él. 

Pasó mucho tiempo en el observatorio dónde fue 
percibida una nueva Tzab-ek y aquella ciudad, a la que 
nuevos hombres llamaron Chuchen Itzá, quedó olvidada, 
también el astrónomo fue olvidado como individuo. Sin 
embargo los nuevos astrónomos reconocieron a sus 
análogos como representantes de una esplendorosa 
civilización que escogió escurrirse silenciosamente por 
las rendijas de la historia, dejando muchos enigmas para 
la posteridad: En la época presente se puede suponer que 
el astrónomo Maya pudo haber observado la extinción de 
una estrella no visible hoy pero que al morir en aquellos 
tiempos, habría refulgido con gran luminosidad. 
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Universidad de Sonora 

UNIDAD REGIONAL CENTRO 
DIVISIÓN DE INGENIERÍA 



Posgrados incluidos en el Padrón Nacional 
de Posgrados de Calidad de CONACYT 



Programa de Doctorado y Maestría 
en Ciencia de Materiales 

Objetivos del Programa 

El objetivo general del programa es formar personal de alto nivel 
académico y con autonomía para llevar a cabo en forma eficiente las 
actividades de generación, aplicación y difusión de conocimiento de 
frontera en el área de Ciencia de Materiales. 

Los objetivos específicos son contribuir a la formación de los 
recursos humanos especializados que el país requiere en Ciencia de 
Materiales. Desarrollar investigación de alto nivel comprobada a partir 
de publicaciones científicas de calidad. Participar de manera efectiva 
en la resolución de problemas tecnológicos que se presenten en el 
medio social y productivo de la región y el país, en el ámbito de las 
líneas de investigación del programa. 

Programa de Doctorado y Maestría en 
Ciencias de la Ingeniería: Ingeniería Química 

Objetivos del Programa 

Formar académicos y profesionales del más alto nivel en ingeniería 
y sus campos afines, útiles a la sociedad. Promover la práctica 
profesional de calidad en ingeniería y campos afines. Contribuir a la 
solución de problemas nacionales. Realizar investigación para generar 
nuevos conocimientos, métodos y criterios en ingeniería y sus campos 
afines. Desarrollar tecnología. 

Programa de Posgrado en Ingeniería Industrial 

Objetivos del Programa 

El objetivo general es formar recurso humano orientado a entender, 
analizar y resolver problemas complejos dentro de organizaciones 
dedicadas a la producción de bienes o servicios. 

Los objetivos específicos son profundizar en el estudio de áreas 
de interés de la Ingeniería Industrial, acorde a las necesidades 
que se deseen resolver por parte de la institución y de la sociedad. 
Hacer énfasis en la aplicación del conocimiento y la tecnología para 
resolver problemas que afectan la calidad y productividad de las 
organizaciones que producen bienes o prestan servicios. Impulsar 
la mejora de la eficiencia y eficacia de las organizaciones sobre todo 
aquellas consideradas como pequeñas y medianas (PYMES) que 
permita impulsar el desarrollo del estado y del país. 



Programa de Posgrado en Ciencias 
de la Ingeniería: Ingeniería Civil 

Objetivos del Programa 

El objetivo general del programa es formar personal del más 
alto nivel académico y profesional en el campo de la Construcción 
y la Valuación, capaz de realizar investigación para generar nuevos 
conocimientos o desarrollos que contribuyan a la solución de 
problemas de interés regional y nacional. 

Programa de Especialidad 
en Desarrollo Sustentable 

Objetivos del Programa 

Formar profesionales con una perspectiva integral de las áreas 
técnica, administrativa, económica, social, ambiental y ética de 
acuerdo a los modelos modernos del desarrollo 

Programa de Maestría en Sustentabilidad 

Objetivos del Programa 

El objetivo general del programa de maestría es coadyuvar con 
la sociedad en la transición hacia estilos de vida más sustentables, 
mediante la formación de profesionistas con una perspectiva holística 
y ética; con los conocimientos y habilidades para planear, implantar y 
evaluar intervenciones que conlleven a prevenir, eliminar y/o reducir 
los riesgos e impactos ambientales y ocupacionales en el área de su 
profesión y/o desarrollo en beneficio de su comunidad. 
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PROGRAMA DE MAESTRIA Y DOCTORADO 
EN CIENCIAS EN MATEMÁTICAS 

Objetivos del Programa 

Formar personal con alto nivel académico, para el estudio autónomo 
y comprensión de los resultados y avances de la matemática moderna y 
sus aplicaciones, con capacidad para realizar labores de asesoría y apoyo 
matemático en proyectos de investigación y desarrollo en disciplinas de 
carácter científico y tecnológico. A nivel de doctorado, además de los 
anteriores, el objetivo principal es capacitar para la realización de investigación 
original e independiente en el campo de las matemáticas y sus aplicaciones. 

PROGRAMA DE MAESTRÍA EN CIENCIAS 

CON ESPECIALIDAD EN MATEMÁTICA EDUCATIVA 

Objetivos del Programa 

El propósito de este programa académico es formar Especialistas en 
Matemática Educativa, cuya práctica profesional como docentes sintetice una 
serie de competencias éticas, ideológicas, de conocimiento teórico y sentido 
práctico, para atender los problemas de la educación matemática escolar en 
los ámbitos de la investigación, el desarrollo docente y la proyección social 
de sus egresados. 

PROGRAMA DE MAESTRÍA EN CIENCIAS EN FÍSICA 

Objetivos del Programa 

Preparar personal con un profundo conocimiento de las leyes 
fundamentales de la Física, para ejercer labores docentes a nivel licenciatura, 
posgrado y especialidad. Capacitar personal para ejercer actividades 
asociadas a la generación de conocimiento en diferentes áreas de la Física. 
Formar personal con una alta capacidad de trabajo e innovación lo cual le 
permita incidir sobre la solución de aquellos problemas de su competencia 
relacionados con el sector productivo de la región y del país. Capacitar 
personal en la utilización eficiente de diverso equipo con un alto desarrollo 
científico y tecnológico. Desarrollar programas multidisciplinarios de 
investigación científica y tecnológica. 

PROGRAMA DE MAESTRÍA EN CIENCIAS GEOLOGÍA 

Objetivos del Programa 

El objetivo general es el de formar Maestros en Ciencias de alto nivel 
de calidad y excelencia, capaces de participar en el análisis y solución de 
problemas geológicos tradicionales y de frontera, utilizando el método 
científico, y respondiendo con oportunidad a las necesidades del desarrollo 
científico y tecnológico del país. Se ofrece un espacio de formación de 
alto nivel académico en torno a la problemática derivada de tres grandes 
áreas de aplicación de las Ciencias de la Tierra como son: los yacimientos 
minerales, la geología regional, y la hidrogeología-geología ambiental. Con 
ello se pretende contribuir al fortalecimiento de la capacidad profesional 
y/o investigativa de los profesionales que se desempeñan en estas áreas del 
conocimiento. 




http://www.investigacion.uson. 
mx/programas-posgrado.htm 
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